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Surface  recombination  and  channel  have  dominated  the  electrical  characteristics, 

performance  and  reliability  of  p/n  junction  diodes  and  transistors.  This  dissertation  uses 

a  sensitive  direct-current  current  voltage  (DCIV)  method  to  measure  base  terminal 

currents  (IB)  modulated  by  the  gate  bias  (VGB)  and  forward  p/n  junction  bias  (VPN)  in  a 

MOS  transistor  (MOST).    Base  terminal  currents  originate  from  electron-hole 

recombination  at  Si/Si02  interface  traps.  Fundamental  theories  which  relate  DCIV 

characteristics  to  device  and  material  parameters  are  presented.  Three  theory-based 

applications  are  demonstrated  on  both  the  unstressed  as  well  as  hot-carrier-stressed 

MOSTs:  (1)  determination  of  interface  trap  density  and  energy  levels,  (2)  spatial  profile 

of  interface  traps  in  the  drain/base  junction-space-charge  region  and  in  the  channel 


region,  and  (3)  determination  of  gate  oxide  thickness  and  impurity  doping 
concentrations.  The  results  show  that  interface  trap  energy  levels  are  discrete,  which  is 
consistent  with  those  from  silicon  dangling  bonds;  in  unstressed  MOS  transistors 
interface  trap  density  in  the  channel  region  rises  sharply  towards  source  and  drain,  and 
after  channel-hot-carrier  stress,  interface  trap  density  increases  mostly  in  the  junction 
space-charge  region. 

As  the  gate  oxide  thins  below  3nm,  the  gate  oxide  leakage  current  via  quantum 
mechanical  tunneling  becomes  significant.  A  gate  oxide  tunneling  theory  which  refined 
the  traditional  WKB  tunneling  probability  is  developed  for  modeling  tunneling  currents 
at  low  electric  fields  through  a  trapezoidal  Si02  barrier.  Correlation  with  experimental 
data  on  thin  oxide  MOSTs  reveals  two  new  results:  (1)  hole  tunneling  dominates  over 
electron  tunneling  in  p+gate  p-channel  MOSTs,  and  (2)  the  small  gate/drain  overlap 
region  passes  higher  tunneling  currents  than  the  channel  region  under  depletion  to 
flatband  gate  voltages.  The  good  theory-experimental  correlation  enables  the  extraction 
of  impurity  doping  concentrations,  which  complements  the  DCIV  method. 

Two  fundamental  theories  of  interband  tunneling  are  developed  to  correlate  with 
the  VGB  dependence  of  drain/base  p/n  junction  currents:  (1)  direct  tunneling  at  the 
drain/base  junction  perimeter  with  and  without  the  quantization  effects  in  the  base 
surface  accumulation  layer,  and  (2)  interface  trap  assisted  tunneling  in  the  gate/drain 
overlap  region.  The  second  theory  gives  better  correlation,  which  is  further  supported 
by  the  DCIV  peaks  originated  from  interface  traps  in  the  gate/drain  overlap  region. 


VI 


CHAPTER  1 
INTRODUCTION 

One  of  the  most  significant  advancements  in  silicon  transistor  technology 
development  was  the  passivation  of  silicon  surface  by  a  thermally-grown  oxide  layer 
some  forty  years  ago  [1,2].  Today,  the  metal-oxide-semiconductor  (MOS)  transistor  has 
become  the  most  important  building  block  of  ultra-large-scale-integrated  (ULSI) 
circuits.  In  MOS  transistors,  electrical  currents  conduct  in  silicon  surface  channels 
directly  underneath  the  silicon/silicon  dioxide  interface.  The  mismatch  of  amorphous 
thermal  Si02  and  crystalline  Si  at  the  interface  gives  rise  to  point  defects  such  as  the 
unpaired  silicon  dangling  bond  which  acts  as  an  electronic  trap.  The  intrinsic  interface 
trap  density  after  the  thermal  oxidation  step  can  be  as  high  as  5xl012cm-2, 
corresponding  to  a  fractional  occupancy  of  0.6%  of  a  (100)  silicon  plane,  most  of  which 
are  passivated  by  hydrogenation  in  a  forming  gas  annealing  step,  to  give  a  residual 
active  interface  trap  density  on  the  order  of  1010cm-2  in  as-manufactured  MOS 
transistors.  However,  during  transistor  operation,  passivated  interface  traps  can  be 
reactivated  and  new  interface  traps  can  be  generated  which  degrade  MOS  transistor 
performance  by  lowering  its  subthreshold  slope  and  increasing  its  leakage  current. 

Traditional  small-signal  measurement  techniques  such  as  the  MOS  capacitance- 
voltage  method  can  only  resolve  interface  trap  density  higher  than  about  10ncm-2  and 
can  not  detect  manufacturing  residual  interface  traps  on  state-of-the-art  MOS  transistors. 


The  seminal  work  on  surface  recombination  and  channel  by  Sah  [3,4]  was  recently 
extended  by  Neugroschel  and  Sah  [5]  in  1995  to  develop  a  sensitive  technique  known  as 
the  DCIV  methodology  that  measures  the  DC  recombination  current  of  injected 
minority  carriers  at  interface  traps.  The  sensitivity  is  achieved  by  forward-biasing  a  p/n 
junction  in  a  MOS  transistor  structure  to  exponentially  raise  the  injected  minority  carrier 
concentration.  In  the  past  five  years,  many  DCIV  application  papers  and  reports  have 
been  published  which  include  delineation  of  interface  trap  generation/annealing  kinetics 
on  electrically-stressed  transistors  [6,7]  and  diagnosis  of  transistor  design  and 
manufacturing  processes  on  pre-stress  transistors  [8-10].  This  dissertation  analyzes  the 
theoretical  basis  of  the  DCIV  method  and  develops  theory-based  new  applications  for 
characterization  of  scaled  MOS  transistors  with  very  thin  gate  oxide. 

In  chapter  2,  the  gate-voltage  controlled  interface  trap  recombination  current 
from  the  base-terminal  is  analyzed  using  the  Shockley-Read-Hall  steady-state 
recombination  kinetics  to  provide  analytical  expressions  for  DCIV  lineshape,  linewidth, 
peak  voltage  and  peak  amplitude.  It  is  predicted  by  the  theory  that  DCIV  peaks  in  the 
intrinsic  to  flat  band  gate  voltage  range  originate  from  interface  traps  located  in  the 
channel  area,  while  additional  peaks  in  the  surface  accumulation  range  originate  from 
interface  traps  covering  the  gated  p/n  junction  space-charge  region.  A  new  application 
to  determine  interface  trap  density  and  energy  level  from  the  forward  p/n  junction  bias 
dependence  of  the  DCIV  peak  current  is  proposed  and  demonstrated  on  both  unstressed 
and  channel-hot-carrier  stressed  MOS  transistors.  The  distortion  of  DCIV  lineshape 
from  minority  carrier  injection  level  and  diffusion  is  also  analyzed. 


In  chapter  3,  spatial  distribution  of  interface  traps  in  the  gated  p/n  junction  space 
charge  region  (JSCR)  is  shown  by  theory  to  affect  the  DCIV  lineshape  in  the  gate 
voltage  range  between  the  weak  inversion  of  the  channel  area  and  the  weak  inversion  of 
the  gate-overlapped  drain  area.  In  this  gate  voltage  range,  the  recombination  rate  per 
unit  interface  trap  is  sharply  peaked  in  the  JSCR  due  to  rapid  lateral  variation  of 
minority  carrier  concentrations.  Based  on  this  discovery,  an  algorithm  is  developed  to 
determine  spatial  density  profile  of  interface  traps  which  has  a  resolution  of  ~3nm.  For 
interface  traps  located  outside  the  JSCR  in  either  the  channel  area  or  the  drain  area,  the 
gate  voltage  that  gives  the  peak  recombination  current  depends  only  on  transistor 
parameters  such  as  the  channel  and  drain  doping  concentrations.  Thus  these  peak 
voltages  provide  a  means  to  probe  the  MOS  transistor  design  parameters.  An  analytical 
three-parameter  formula  is  given  to  fit  the  forward  junction  bias  dependence  of  DCIV 
peak  voltage  for  the  extraction  of  gate  oxide  thickness,  impurity  doping  concentration 
and  flatband  voltage. 

MOS  transistor  has  been  aggressively  downscaled  in  the  past  25  years  to 
increase  its  speed  while  reducing  its  manufacturing  cost  at  the  same  time.  This  trend 
will  continue  in  the  near  future  [1 1-12].  In  each  technology  generation,  gate  oxide 
thickness  is  scaled  almost  in  proportion  to  the  channel  length,  while  the  Si/Si02 
interfacial  layer  can  not  be  scaled  and  becomes  a  more  significant  part  of  the  total  Si02 
thickness  in  thin  oxide  MOS  transistors.  The  DCIV  method  will  continue  to  be  a 
powerful  tool  for  evaluating  interface  reliability  as  well  as  for  characterizing  transistor 
design.   On  thin  oxide  MOS  transistors,  DCIV  measurements  show  base-terminal 


current  that  rises  with  gate  voltage,  due  to  (1)  increased  gate  oxide  tunneling  current  at 
both  inversion  and  accumulation  gate  voltages,  and  (2)  increased  current  flow  through 
forward-biased  drain/base  p/n  junction  at  accumulation  gate  voltages.  Quantitative 
understanding  of  these  two  phenomena  will  not  only  provide  accurate  baselines  for 
DCIV  peak  analysis,  but  also  enable  new  techniques  for  the  characterization  of  scaled 
MOS  transistors. 

Recent  experimental  data  on  ultrathin  oxide  (1.0nm-3.0nm)  MOS  transistors  [13] 
prompted  us  to  develop  a  gate  oxide  tunneling  theory  that  takes  into  account  all  three 
tunneling  particle  species  (electrons,  holes  and  valence  electrons),  which  is  presented  in 
chapter  4.  Theory-experimental  correlation  will  prove  two  results:  (1)  hole  is  the 
dominant  tunneling  carrier  in  P+  gate  p-channel  MOS  transistors,  which  was  not 
recognized  previously,  and  (2)  the  much  smaller  gate/drain  overlap  area  passes  higher 
current  than  the  channel  area  in  the  gate  voltage  range  between  flatband  and  strong 
inversion  [10,13].  Experimental  tunneling  currents  from  gate,  source/drain,  well  and 
substrate  terminals  will  be  matched  with  the  tunneling  theory  from  which  transistor 
design  parameters  such  as  oxide  thickness,  impurity  doping  concentrations  in  the 
channel  and  drain  areas  are  determined. 

In  chapter  5,  physical  origins  for  the  increased  DCIV  base-terminal  current  at 
accumulation  gate  voltages  are  investigated.  Two  theories  of  interband  tunneling  are 
considered:  (1)  direct  tunneling  at  the  drain/substrate  junction  perimeter  including  the 
effects  of  quantized  surface  accumulation  layer,  and  (2)  interface  trap  assisted  tunneling 
in  the  gate-overlapped  drain  area.  The  second  theory  has  a  lower  threshold  gate  voltage 


and  predicts  a  Shockley-Read-Hall  recombination  current  peak  from  interface  traps 
located  in  the  overlap  area,  both  of  which  are  supported  by  experimental  data. 
We  will  summarize  and  conclude  the  dissertation  in  chapter  6. 


CHAPTER  2 

INTERFACIAL  ELECTRONIC  TRAPS 

IN  SURFACE  CONTROLLED  TRANSISTORS 

2.1  Introduction 

Surface  recombination  and  channel  have  dominated  the  electrical  characteristics, 
performance,  and  reliability  of  p/n  junction  diodes  and  transistors.  Due  to  the 
technological  importance,  extensive  research  efforts  have  been  undertaken  to  study 
interfacial  electronic  traps  at  the  Si/Si02  interface  and  to  delineate  their  microscopic 
(atomic)  origins  [14].  Two  kinds  of  density-of-states  (DOS)  spectra  were  repeatedly 
observed  in  the  energy  gap  of  silicon:  two  DOS  peaks  near  the  midgap  and  a  U-shaped 
distribution  that  rises  towards  the  two  band  edges.  The  discrete  energy  levels  are 
characteristic  of  point  defects  at  the  Si/Si02  interface,  while  the  U-shaped  band-tail 
states  are  from  random  variations  of  Si-O  bond  length  and  bond  angles  [15].  These 
conclusions  were  supported  by  early  works  of  molecular-orbital  calculations  of  Si/Si02 
interfacial  electronic  structures  [16,17]. 

The  peaked  DOS  were  correlated  to  the  trivalent  silicon  dangling  bond  (Pb 
center)  on  oxidized  (111)  silicon  [18-20]  and  its  two  variations  (Pb0  and  Pbl  centers)  on 
oxidized  (100)  silicon  by  the  electron  spin  resonance  (ESR)  experiments  [19,21]. 
Recent  studies  using  ESR  on  oxidized  silicon  and  spin-dependent  recombination  (SDR) 


on  MOS  transistors  reveal  that  the  two  Pb  centers  on  (100)  Si  have  different  generation 
and  annealing  kinetics.  Electric  field  stress  including  the  non-uniform  channel  hot 
carrier  stress  [22]  and  the  uniform  Fowler-Nordheim  stress  [23]  creates  mostly  Pb0 
centers,  while  both  centers  can  be  either  passivated  or  generated  by  atomic  hydrogen 
[24,25].  The  structural  difference  among  the  three  Pb  centers  can  be  grossly  described 
by  the  direction  of  the  single  unpaired  Si  sp3  hybrid  orbital.  It  points  to  a  normal  <1 1 1> 
direction  for  both  the  Pb  center  on  (1 1 1)  Si  and  the  Pb0  center  on  (100)  Si,  while  it  points 
nearly  along  <21 1>  for  the  Pbl  on  (100)  Si  [26].  Multiple  DOS  peaks  of  yet  another 
origin  were  observed  by  10.2eV  light-induced  hole  injection  in  MOS  capacitors 
annealed  in  oxygen  ambient  which  could  be  from  oxygen  dangling  bonds  at  the  Si/Si02 
interface  [27]. 

While  electron  spin  resonance  has  been  used  to  identify  the  atomic  structures  of 
interface  traps  on  large  area  (~lcm2)  oxidized  silicon  wafers,  other  techniques  are  more 
traditionally  used  to  monitor  interface  traps  in  MOS  transistors  and  MOS  capacitors  in 
manufacturing  facilities.  This  includes  MOS  capacitance  [28,29]  and  conductance  [30] 
methods  and  more  recently  the  charge-pump  method  [31]  and  DCIV  method  [5]. 
Routine  manufacturing  processes  have  reduced  interface  trap  areal  density  to  below 
1010cm-2  by  slow  cooling  after  high-temperature  oxidation  steps  and  by  post-oxidation 
annealing  in  hydrogen.  At  this  low  density,  the  DCIV  technique  has  its  unique 
advantage  in  sensitivity  over  the  other  small-signal  and  transient  techniques. 

The  principle  of  DCIV  is  the  use  of  a  surface  potential  controlling  gate  terminal 
voltage,  VGB,  to  modulate  the  base-terminal  DC  current,  IB,  from  electron-hole 
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recombination  at  the  Si02/Si  interface  traps  [3,4].  The  DCIV  methodology,  which  gives 
a  family  of  IB-VGB  curves  from  MOSTs  and  gated  p/n  junctions,  was  recently  developed 
by  Neugroschel  and  Sah  [5]  to  monitor  the  density  of  oxide  and  interface  traps  generated 
by  hot  carriers  (HC)  in  bipolar-junction  and  metal-oxide-field-effect  transistors  (BJTs 
and  MOSTs)  in  order  to  determine  the  transistor  failure  rate  and  10-year  operation  time- 
to-failure  (TTF0p)  voltage  [32].  Its  application  was  extended  to  monitor  the  degradation 
rate  of  transistors  stressed  at  high  current-densities  and  low  voltages  (HJ),  and  to 
delineate  the  origin  of  hydrogen-related  interconnect  degradations  [6,7].  Recent 
demonstrations  showed  that  the  DCIV  method  could  serve  as  a  highly  sensitive 
diagnostic  monitor  for  transistor  designs  [8]  and  manufacturing  processes  [8,33].  The 
DCIV  lineshape  (i.e.  IB  vs  VGB  around  the  Ie-peak)  was  a^s0  use^  t0  obtain  the  electronic 
and  quantum-mechanical  properties  of  the  residual  Si02/Si  interface  traps  on  state-of- 
the-art  thermally  grown  thin  gate  oxides  [34]. 

This  chapter  reports  the  elementary  theory  of  the  IB-VGB  characteristics  to  help 
quantify  the  applications  of  the  DCIV  methodology  for  the  extraction  of  fundamental 
and  application-specific  properties  of  transistors  and  interconnects  and  their  materials, 
such  as  the  physical  (spatial  location  and  density)  and  electronic  (quantum  density  of 
states)  properties  of  the  residual  and  stress-generated  interface  and  oxide  traps,  the 
dopant  impurity  concentration  profiles,  and  the  hydrogen  sources.  Analytical  solutions 
and  their  physical  models  are  presented  first  to  illustrate  the  effects  of  material  and 
interface  trap  parameters  on  the  IB-VGB  lineshape  and  on  the  VGB  ak  value  at  the 
IB.peak.   Application  examples  are  then  given  for  the  process-residue  and  hot-carrier- 


generated  interface  traps  on  thin  thermally  oxidized  gate-oxide  of  production  MOS 
transistors. 

2. 1  Carrier  Recombination  in  MOS  Transistor  Channel  Region 

In  DCIV  measurements,  excess  minority  carriers  are  injected  by  a  forward- 
biased  p/n  junction  into  the  MOS-gated  Si02/Si  interface  which  covers  MOST's  channel 
region  (denoted  by  MCR  for  mid-channel  region)  and  drain-source  junction  space- 
charge  region  (JSCR).  The  p/n  junction  can  be  located  either  away  from  the  gated 
interface,  such  as  the  well/substrate  junction  (bottom-emitter,  BE)  and  a  remote  p/n 
junction  (remote-emitter,  RE)  or  under  the  gated  interface,  such  as  the  source  or  drain 
p/n  junction  (source-emitter,  SE,  and  drain-emitter,  DE,  or  top-emitter  with  source-drain 
tied  to  the  same  forward  bias  voltage).  The  Ib_^gb  theory  is  presented  for  these 
injection  bias  configurations.  The  Ib.jscr  current  is  especially  important  for  factory 
applications  to  monitor  channel-hot-carrier-generated  interface  traps,  while  the  Ib_mcr 
current  is  also  useful  for  determining  the  electronic  properties  of  the  process-residual 
and  stress-generated  interface  traps.  Modifications  of  the  theory  by  DCIV  bias 
configurations,  diffusion,  dopant-impurity  profiles  and  2-dimensional  geometries  are 
then  discussed. 

Figure  2.1  is  the  energy  band  diagram  of  the  metal-oxide-silicon  structure  with  a 
forward-biased  p/n+Si  bottom  emitter  (BE)  junction.  It  is  labeled  in  detail  to  help 
describe  the  approximations  of  the  analyses  as  follows.  All  voltages  are  normalized  by 
kT/q.  The  BE  is  forward-biased  at  UBE  while  the  DC  voltage  applied  to  the  gate 
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Fig.  2.1  Energy  band  diagram  and  cross-sectional  view  of  a  gated  Si02/p-Si/n+Si 
structure  with  the  p/n+ junction  forward  biased.  UN  and  Up  are  respectively 
the  electron  and  hole  quasi-Fermi  potentials  normalized  to  thermal  voltage 
(kT/q).  Labeled  on  the  figures  are  the  terminologies  (gate,  base,  and 
emitter),  dimensions  [XB0  and  XB(VGB,VEB)],  and  the  normalized  (to  kT/q) 
voltages  and  quasi-Fermi-potentials  (UGB,  UFB=Flat-Band,  and  UP,  UN). 


11 

relative  to  the  p-Si  base  is  UGB.  UN  and  Up  are  electron  and  hole  quasi-Fermi  potentials 
in  p-Si  and  their  difference,  known  as  the  quasi-Fermi-potential  split,  is  UPN=UP-UN 
which  is  slightly  smaller  than  the  applied  BE  voltage,  UBE.  A  formula  for  AUPN  will  be 
derived  later.  In  the  quasi-neutral  region  of  p-Si  base  (BQNR),  the  quasi-Fermi- 
potentials  are: 

Up   =    +    loge{  [  (NjA+4n?expUPN)1/2+NAA]/2ni>  (2.1) 

UN  =    -    loge{  [  (NjA+4n?expUPN)1/2-NAA]/2ni}  (2.2) 

The  electron  and  hole  concentrations  at  the  Si02/Si  interface  (Ns  and  Ps)  are  modulated 

by  the  gate  voltage  via  bending  the  Si  energy  band.  The  total  energy  band  bending  is 

denoted  by  the  surface  potential  Us=U(x=0,y). 

There  are  four  fundamental  transition  processes  between  the  continuous  band 
states  of  silicon  crystal  and  the  localized  trap  states  with  an  energy  level  ET  in  the  Si 
energy  gap,  as  shown  in  the  transition  energy  band  diagram  of  figure  2.2.  The  rate 
(event/second)  of  the  four  processes  can  be  conveniently  described  by:  (a)  electron 
capture  from  the  conduction  band  at  cnn(NTT-nT),  (b)  electron  emission  to  the 
conduction  band  at  ennT,  (c)  hole  capture  from  the  valence  band  at  c  pnT,  and  (d)  hole 
emission  to  the  valence  band  at  e  (NTT-nT).  Here  n  and  p  are  electron  and  hole 
concentrations  in  the  conduction  band  and  valence  band  respectively,  NTT  is  the  total 
density  and  nT  is  the  electron-occupied  density  of  trap  states,  and  e's  and  c's  are 
emission  and  capture  rate  coefficients  of  the  four  processes  which  depend  on  the  energy 
levels  of  both  the  trap  state  and  the  band  state.  These  transition  processes  are  mostly 
thermal  involving  emission  and  capture  of  phonons  (hco).  Since  the  indirect  energy  gap 
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Fig.  2.2  A  transition  energy  band  diagram  showing  four  fundamental  transition 
processes  between  a  band  state  and  a  gap  state  in  silicon:  capture  of  a 
condution  band  electron  by  the  trap  (a),  emission  of  an  electron  to  the 
conduction  band  from  the  trap  (b),  capture  of  a  valence  band  hole  by  the 
trap  (c)  and  emission  of  a  hole  to  the  valence  band  from  the  trap  (d).  The 
volume  density  of  band  electrons,  band  holes,  electron-occupied  traps  and 
total  traps  are  n,  p,  nT  and  N^  respectively.  The  rates  of  the  four  process 
are  shown  in  terms  of  e's  and  c's.  Purely  thermal  emission  and  capture 
processes  involve  multiple  phonons. 
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of  Si  is  1.12eV  and  the  maximum  optical  phonon  energy  in  Si  is  only  about  60meV,  the 
thermal  capture  and  emission  processes  could  involve  about  ten  phonons  for  mid-gap 
trap  levels. 

Theoretical  calculation  of  a  multi-phonon  process  was  first  done  by  Huang  and 
Rhys  [35]  and  later  refined  by  Huang  [36].  In  contrast  to  the  lattice  scattering  process 
which  can  be  treated  by  taking  the  linear  coupling  of  electron  Coulomb  potential  and 
lattice  vibration  as  the  perturbation  to  calculate  the  scattering  rate,  the  multi-phonon 
process  can  not  be  treated  in  such  a  way  since  the  linear  coupling  term  could  only  result 
in  the  emission  and  capture  of  a  single  phonon.  Instead,  the  breakdown  of  the  adiabatic 
approximation  [37]  was  used  as  a  perturbation  to  compute  the  transition  rate  between 
the  band  state  and  the  trap  state  while  the  linear  coupling  term  was  used  in  the  time- 
independent  perturbation  theory  to  compute  the  lattice  distortion/relaxation  [38]  in  the 
initial  and  final  states.  The  first-principles  calculation  of  the  capture  cross-sections  in  a 
multi-phonon  process  is  rather  laborious.  For  the  purpose  of  developing  a  theory  for  the 
DCIV  methodology,  it  would  suffice  to  use  the  phenomelogical  kinetic  theory 
developed  by  Shockley  and  Read,  and  by  Hall  [39]  which  treats  the  fundamental  capture 
and  emission  rates  as  constants  independent  of  kinetic  energies  of  band  electrons  and 
holes. 

The  steady-state  areal  rate,  Rss,  of  electron-hole  recombination  at  a  discrete- 
energy-level  interface  trap  (the  equilateral  triangle  in  Fig.  2.1),  located  at  UTI  from  the 
intrinsic  Fermi  potential,  with  a  density  of  NIT  (#/cm2),  is  given  by  the  Shockley-Read- 
Hall  formula: 
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(CnsCpsNSPS     "     enseps> 
Rss    =    .    NIT    =    RSS1-NIT  (2.3) 

CnsNS     +     ens     +     CpsPS     +     eps 

cns,  cps,  ens  and  eps  are  the  electron-hole  capture-emission  rate  coefficients  at  the 
interface  traps.  From  detailed  balance  near  thermal  equilibrium,  c's  and  e's  are  related 
by  ens  =  cnsnjexp(UTI)  and  e  s  =  c    njexp(-UTI).   Using  Boltzmann  representation, 
Ns=njexp(Us-UN)  and  Ps=njexp(UP-Us),  then  (1)  becomes 
(1/2)  (cnscps)1/2ni[exp(UPN)-l] 


-ps* 

exp(UPN/2)  cosh(Us*)    +    cosh(UTI*) 


Rss    -    •    NIT   -    RSsi*nit 


(2.4) 
where  Us*(x=0)  =  Us+loge(cns/cps)1/2-(UP+UN)/2,  and  UTI*  =  UTI  -  loge(cns/cps)1/2.  The 
expression  (2.4)  is  exact  with  no  approximations  other  than  the  thermal  Boltzmann 
distribution  with  lattice  temperature  T.  It  immediately  shows  the  presence  of  a  peak  at 
Us*=0.  The  peak  magnitude  and  the  surface  potential  at  the  peak  are 
(1/2)  (cnscps)1/2ni[exp(UPN)-l] 

RSS-peak    =  *     NIT 

exp(UPN/2)    +    cosh(UTI*) 

=    RSSl-peak'NIT  (2-5) 

Us-peak     =    (UP+UN)/2    -    loge(cns/cps)1/2  (2.6) 

The  peak  formula  (2.6),  was  derived  by  Sah  et  al.  in  1957  [40].  The  half  width  at  half 
maximum  (HWHM)  of  the  Rss-Us  line  is: 

AUS    =    cosh"1  [2  +  exp(-UPN/2)  cosh  (UTI* )  ]  (2.7) 

«    cosh-1  (2),  for   UPN   >    2(|UTI*|     +    4) 

=    UTI*    -     (UPN/2)  for   UPN    <    2(|UTI*|     -    4) 
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The  HWHM  drops  to  the  narrowest,  (cosh_12)*(kT/q)=1.32kT/q=34mV,  at  high  VPN 
and  widens  to  UTI*  at  low  VPN.  Results  (2.4)-(2.9)  are  valid  at  all  injection  levels  in 
both  p-  and  n-Si.  For  single-discrete  level  interface  trap,  they  predict  a  peak  body- 
terminal  recombination  current,  given  by  Ie-peak  =  Q^ssi-peaiJ  NITdydz  that  is 
proportional  to  exp(UPN/n)  with  a  transition  of  n  from  1  to  2  over  a  small  range  of 
forward  bias,  about  AUPN=4  or  AVPN  =  100mV.  This  result  gives  two  important 
consequences.  First,  for  a  single  species  of  interface  traps  with  a  spatially  varying 
concentration,  the  DCIV  lineshape  will  not  change  due  to  the  y-distribution  of  NIT. 
Second,  if  there  are  many  interface  trap  species,  each  with  a  different  energy  level  E^*, 
IB  is  then  just  the  superposition  of  the  bell-shaped  Rss  curves  one  from  each  of  the 
energy  levels  or  species  with  different  peak  locations  in  VGB,  resulting  in  a  broadening 
and  generally  bell-shaped  multi-peak  DCIV  curve  [34]  and  an  n-slope  between  1  and  2 
over  a  larger  range  of  VPN. 

For  low  injection  levels,  traditionally  defined  as  N  <  NAA/10  in  p-Si,  we  have 
Up=UF>0  for  p-Si  (UN=UF<0  for  n-Si),  where  UF  is  the  majority  carrier  Fermi  potential. 
This  is  the  common  application  range  of  the  DCIV  methodology.  At  low  levels, 
equation  (4)  becomes: 

Us-peak   =      UF    -    UPN/2    -    loge(cns/cps)1/2  (p-Si)  (2.8) 

=      UF   +   UPN/2    -    loge(cns/cps)1/2  (n-Si).  (2.9) 

Thus,  the  peak  is  in  the  flat-band  to  intrinsic  gate-voltage  range  (0<Us<UF,  for  p-Si).  In 

this  gate-voltage  range,  Sah  [41]  had  given  a  formula  for  VS(VGB)  listed  in  (2.10)  below 

which  we  shall  use  to  obtain  the  analytical  solutions  for  the  IB-VGB  linewidth. 
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vs  =  vGB  -  vFB  - 

2.sign(Vs)  .V^t  [l+(VGB-VFB-kTAD/q) /V^]172  -  1}  (2.10) 
VFB  is  the  flatband  voltage,  VAA=esqNAA/2CoX  where  C0x=e0x^ox  *s  tne  oxide 
capacitance  per  unit  area,  and  AD=1  in  this  range  [41].  The  last  term  in  (2.10)  is  the 
voltage  drop  across  the  oxide  layer,  whose  error  is  negligible  when  using  AD=1  until  Vs 
is  about  kT/q  from  the  flatband  condition.  The  VGB  at  the  peak  can  be  solved  from 
(2.10)  in  terms  of  V^^  given  by  (2.8)  and  (2.9): 

^GB-peak     ~     ^FB     =    ^S-peak     + 

2.sign(Vs_peak).(VAA.|Vs_peak|)1/2  (2.11) 

Thus,  the  DCIV  peak  location  is  determined,  via  Vs  ak,  by  the  substrate  impurity 
concentration,  trap  parameters  in  terms  of  loge(cns/c  )  and  oxide  thickness  (from  VAA). 
The  half-width  at  half  maximum  (HWHM)  of  the  DCIV  peak  is: 


AVGB    =    AVS    +    2J^    [    ^l|Vs_peak|  -      ^l|Vs_peak|-AVs    ] 

(flat-band   side)  (2.12) 


=   AVS   +    2fv~    [    ^|Vs_peak|+AVs    -      fj 


^O-npaV 


S-peak i         VS  'I     S-peak 

(intrinsic    side)  (2.13) 

Since  Rss  is  symmetric  around  the  peak  in  Vs,  (2.12)  and  (2.13)  show  that  the  DCIV  IB- 
VGB  lineshape  is  asymmetric  and  slightly  wider  on  the  flat-band  side  of  the  peak  than 
on  the  intrinsic  side.  The  difference  is  on  the  order  of  0.5(VAAVS  ak)1/2(AVs/Vs  ak)2, 
which  is  more  pronounced  in  transistors  with  thick  oxide  and  high  surface  impurity 
concentration. 
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At  high  injection  levels  with  N  >  10NAA,  we  have  UP=-UN  or  the  electron  and 
hole  concentrations  are  nearly  equal  in  the  BQNR,  and  the  maximum  surface 
recombination  rate  is  near  the  flat-band.  The  exact  result  is  Us_peak  =  loge(c  /cns)1/2 
derived  from  (2.6).  So  the  peak  could  not  move  deeply  into  accumulation  range  even  at 
extremely  high  forward  bias  since  c  /cns  is  not  likely  to  be  outside  of  the  range  of  0.01 
to  100  for  any  physically  realistic  bound  state.  The  linewidth  at  the  high  injection  level 
limit  is  obtained  from  the  general  MOS  voltage  equation  in  p-Si,  VGB(VS): 
VGB-VFB-VS    =    sign(Vs)  •  (kT/q)  • 

{2UAA[  (eUs-l-Us)  (a-l)  +  (e-Us-l+Us)  (a+1)  ]  }1/2  (2.14) 
where  a  =  [1  +  4n?exp  (UPN) /NAA]  1/2  .  At  high  injection  levels,  a  = 
2exp[(UPN/2)-UF]  »  1.  Only  one  term  in  the  {  }  of  (2.14)  will  dominate,  depending 
on  the  sign  of  Us: 

vgb-vfb-vs   *    2  (kT/q)  (UAA)1/2sign(Us)  . 

exp[  (  |Us|+UPN/2-UF)/2]  (2.15) 

The  half-width  of  the  DCIV  line  is  then 

AVGB   "   Avs   = 

2  (kT/q)  (e+AUs/2-l)  [XJjJ  cps/cnsexp  (UPN/2-UF)  ] 1/2 

(Ss>cns'Vs>0)  (2.16) 

2  (kT/q)  (l-e-AUs/2)  [U,J  cns/cpsexp  (UPN/2-UF)  ] 1/2 

(Ss<cns'Vs<0)  (2.17) 

These  results  show  that  the  Full  Width  Half  Maximum  (FWHM)  has  the  exp(UPN/4) 

dependence  on  forward  bias  at  high  injection  levels.  Thus,  at  low  injection  levels,  the 

IB-VGB  linewidth  can  be  large  which  is  determined  by  the  trap  level,  ETI*.   The 
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linewidth  then  decreases  with  increasing  VPN  (to  about  FWHM=70mV  for  very  thin 
oxide)  until  the  onset  of  high  injection  level  condition,  beyond  which  it  increases 
exponentially  with  VPN. 

Figure  2.3  shows  the  theoretical  and  normalized  DCIV  lineshapes  or  the  unit 
steady-state  recombination  rate,  RSS1  vs  VGB.  The  shape  changes  with  increasing 
forward  biases,  VPN,  applied  to  a  n+/p  junction  in  p-Si  with  an  acceptor  impurity 
concentration  of  NAA=1017cm~3,  an  oxide  thickness  of  5.0nm,  and  a  discrete  interface 
trap  level  at  ETI*=200meV.  For  VPN<ETI*/q,  there  is  a  flat-top  region  around  the 
maximum  RSS1.  The  linewidth  gradually  decreases  with  increasing  VPN  and  reaches  the 
high  injection  level  asymptote  as  indicated  by  the  dashed  line  at  VPN=900mV. 

2.3  Carrier  Recombination  in  MPS  Transistor  Junction  Space-Charge  Region 

The  preceding  IB-VGB  theory  applies  to  interface  traps  over  MOST's  mid- 
channel  region.  This  IB  component  will  be  denoted  as  the  Ib_mcr-  The  MCR  formulas 
assumed  that  there  is  no  lateral  variation  of  the  carrier  concentrations,  thus,  the  surface 
energy  band  bending  Us  is  constant,  independent  of  y,  at  a  given  gate  voltage.  For 
interface  traps  over  MOST's  source  and  drain  JSCR,  the  foregoing  analysis  needs  to  be 
extended  to  take  into  account  of  Us(y)  in  the  JSCR.  These  interface  traps  contribute  a 
second  component  to  IB,  denoted  as  IB-jscR'  wnen  source,  drain  or  both  of  them  are 
forward  biased.  Figure  2.4  illustrates  the  cross-sectional  view  and  the  energy  band 
diagram  at  the  Si02/Si  interface,  x=0,  of  a  forward-biased  abrupt  p/n+  junction.  In  Fig. 
2.4,  the  build-in  potential  of  the  p/n+  junction  is  given  by 
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Fig.  2.3  Theoretical  normalized  DCIV  lineshapes  as  a  function  of  minority  carrier 
injection  levels,  VPN.  Parameters  used  are  as  follows.  VPN=0-900mV  at 
lOOmV  intervals  with  OmV  and  900mV  the  asymptotes.  p-Si  with  NAA  = 
1017cm"3.  n^lO^cm-3.  Oxide  thickness  XOx=5.0nm.  Discrete  interface 
trap  at  ETI  =  200meV  from  the  intrinsic  Fermi  level  (~mid-gap)  with 
cns=c=l(T8cm3/s,  and  NIT=low— >0  but  arbitrary  since  normalized. 
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Fig.  2.4  The  x-cross  section  view  of  the  two-dimensional  energy  band  diagram  of  a 
p/n+ junction  on  the  Si02/Si  interfacial  plane,  x=0.  The  one-sided  junction 
space  charge  region  extends  from  y=0  to  y=YEB  in  the  p-Si.  VBI  =  built-in 
potential  of  the  p/n+ junction,  VB  =  p/n+  barrier  height  =  VBI  -  Vs  -  VPN, 
Vs  =  surface  band  bending  outside  the  JSCR,  and  VPN  =  Vp  -  VN  =  forward 
bias  voltage  applied  to  the  p/n-i- junction  disregarding  voltage  or  potential 
drops.  All  labeled  potentials  are  normalized  to  (kT/q)  with  subscripted  (##) 
symbol  U##=qV##/kT. 
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VBI=(kT/q)loge(NAANDD+/nf ).  The  total  normalized  barrier  height  at  the  surface  of  the 
p/n+ junction  is: 

UB=UBI-UPN-USX+USE  (2.18) 

where  Usx  and  USE  are  the  surface  band  bendings  outside  of  the  JSCR  in  the  MCR  and 
in  the  emitter  (n+drain)  region  respectively.  We  use  U(y)  to  represent  the  normalized 
surface  potential  in  the  JSCR  along  the  interface  from  y=0  at  the  junction  boundary  to 
y=YEB  at  the  p-edge  of  the  JSCR,  with  U(YEB)=USX  given  by  boundary  condition. 

For  the  one-sided  p/n+  junction,  the  n+  side  JSCR  is  much  smaller  than  the  p- 
side  JSCR  and  the  lateral  (y-direction)  voltage  drop  across  the  n+  side  JSCR  can  be 
neglected.  However,  the  vertical  (x-direction)  surface  band  bending  in  the  n+drain  is 
significant  at  large  negative  gate  biases  even  for  highly  degenerate  doping 
concentrations.  Figure  2.5  shows  surface  band  bendings  Vsx  in  part  (a)  and  VSE  in  part 
(b)  for  the  gated  p/n+ junction  as  a  function  of  gate  to  p-substrate  voltage  VGB  with  the 
drain  to  p-substrate  forward  bias  VDB=-VPN  as  a  parameter.  The  results  were  obtained 
by  solving  one-dimensional  MOS  capacitor  equations  taking  into  account  of  non- 
equilibrium  minority  carrier  concentrations.  At  VGB<-1  V  when  the  p-substrate  surface 
is  strongly  accumulated,  surface  band  bending  in  n-i-  drain  with  NDD=1019cm-3  becomes 
more  significant  than  that  in  the  p-substrate  with  NAA=1017cm-3.  VSE  reaches  -1.0V  at 
VGB=-3V  under  a  lOOmV  forward  bias,  while  Vsx  is  at  -0.25V  at  the  same  gate  voltage 
and  is  independent  of  forward  bias.  It  shows  that  under  the  surface  accumulation  VGB 
range,  Vsx  is  the  dominant  term  in  eqn.  (2.18)  for  UB  and  therefore  indispensable  in 
evaluating  the  surface  potential  U(y)  in  the  JSCR. 
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Fig.  2.5  Surface  band  bending  in  (a)  p-substrate,  Vsx  and  (b)  n+  drain,  VSE  as  a 
function  of  gate  voltage,  with  the  forward  bias  on  the  drain/substrate 
junction  as  a  parameter  that  varies  from  OV  to  800mV  with  a  step  size  of 
lOOmV.  Transistor  parameters  include  doping  concentrations: 
NAA=1017cm-3  in  p-substrate,  NDD=1019cm-3  in  n+  drain  and  5xl019cm-3 
in  n+  polysilicon  gate,  and  gate  oxide  thickness  of  3.0nm. 
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The  carrier  concentration  at  0<y<YEB  in  the  JSCR  are: 
Ns(y)    =   niexp[U(y)-UN]    =   n^xp  [Usx+U(y) -U(YEB) -UN]  (2.19) 

Ps(y)    =   niexp[Up-U(y)  ]    =   r^exp  [Up-Usx-U  (y) +U  (YEB)  ]  (2.20) 

IB.JSCR  is  then  obtained  by  integrating  Rss,  given  by  (2.3)  over  the  JSCR,  which  gives 

YEB  NIT(y)dyZ 


■^B-JSCR    ~     -'"BJO 


(2.21) 


0  exp  (UPN/2)  cosh  (Us*j)    +   cosh(UTI*) 

where  for  the  p/n+ junction 

IBJ0    =    q(cnscps)1/2(ni/2)  [exp(UPN)-l]        (Ampere)  (2.22) 

Us*j    =    Usx    +    U(y)     -    U(YEB)     -     (Up+UN)/2    +    loge  (cns/cps)  1/2 

(2.23) 

Two  NIT  spatial  distributions  will  be  considered:  (1)  Localized  NIT(y)=NIT5(y- 
YT)  and  (2)  uniform  NIT(y)=  NIT  =  constant.  The  localized  interface  trap  distribution  is 
expected  for  channel-hot-carrier  generated  traps  in  MOSTs  from  Sah's  Si:H  bond 
breaking  model  [42]  in  which  the  generated  silicon  dangling  bonds  are  concentrated  at 
the  location  near  the  drain  junction  towards  the  source  where  the  secondary  hot  carrier 
kinetic  energy  reaches  the  Si:H  bond  energy  (3.1eV). 

For  the  localized  spatial  distribution,  Ib_jscr  nas  a  peak  when  Us*j=0  in  (2.23), 

which  gives: 

Usx-JSCR-pk    =     (UP+UN)/2    -    loge(cns/cps)1/2    -     [U(YT)     -    U(YEB)] 

=    Usx_MCR_pk  -     [U(YT)     -    U(YEB)] 

(2.24) 

This  shows  a  general  result:  the  IB-jscr  Pea^  *s  always  at  the  accumulation  gate-voltage 
side  of  the  Ib-mcr  Pea^-  F°r  interface  traps  located  closer  to  the  p/n+drain  junction 
boundary,  y=0,  [U(YT— >0)-U(YEB)]  in  (2.24)  is  larger  and  the  peak  is  further  into  the 
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accumulation  VGB  range.  Quantitative  results  can  be  obtained  from  analytical  solutions 
when  \dEx/dx\«\dE Jdy\  inside  JSCR  which  we  shall  call  the  abrupt  model,  i.e.  VGB 
does  not  modulate  the  potential  variation  inside  JSCR,  only  outside  in  the  MCR.  Then, 
using  the  depletion  approximation,  U(y)  assumes  the  quadratic  form  illustrated  in  Fig. 3 
given  by 

U(y)     -    U(YEB)     =    (q/kT)  (qNAA/2es)  .  (YEB   -    y)2 

=     (YEB-y)2/2L2)  0<y<YEB  (2.25) 

where  LD  =  [(kT/q)(es/qNAA)]1/2  is  the  Debye  screen  length  in  p-substrate,  U(YEB)=USX 
and  the  JSCR  barrier  height  is  given  by  UB  =  U(0)  -  U(YEB)  =  Y|B/2L§.  The  thickness 
of  the  JSCR,  YEB,  is  determined  from  the  barrier  height  UB.  Using  (2.18)  and  (2.25), 
we  derive  the  following  expression  for  the  surface  potential  at  the  trap  location  y=YT 
inside  the  JSCR  in  terms  of  surface  band  bendings  outside  of  the  JSCR,  Usx  and  USE: 

U(YT)     -    UBI    -    UpN    +    USE    - 

2[UYT(UBI-USX-UPN+USE)  ]1/2    +    UYT  (2.26) 

where  UYT  =  Y|/2Lq.  At  IB.JSCR  peak,  U(YT)  satisfies  equation  (2.6),  which  together 
with  (2.26)  provide  the  following  equation  to  relate  the  two  surface  band  bendings  at  the 
peak: 

usx-JscR-pk   =   UBI-UPN+USE_JSCR_pk-  (1/4UYT)  x 

[UBi-UPN/2-UF+USE_JSCR_pk+UYT+loge(cns/cps)1/2]2  (2.27) 
The  surface  band  bendings  outside  the  JSCR  in  p-substrate  and  n+  drain  are  both 
determined  by  the  gate  voltage  through  the  one-dimensional  MOS  capacitor  equation. 
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Thus  eqn.  (2.27)  provides  a  means  to  numerically  solve  the  two  surface  band  bendings 

at  Ib-jscr  Peak- 

The  result  of  IB-VS  HWHM  (denoted  by  AVS)  in  equation  (2.7)  applies  to  all  NIT 
locations,  while  the  result  of  IB-VGB  HWHM  in  (2.10)  and  (2.1 1)  applies  to  NIT  in  JSCR 
only  if  U(YT)-U(YEB)  <  Usx.MCR.pk  so  that  Usx.JSCR.pk  in  (2.24)  will  not  change  sign, 
which  corresponds  to  NIT  located  near  the  p-edge  of  the  JSCR.  If  the  condition  is  not 
satisfied,  then  AVS  translates  differently  to  AVGB  because  the  peak  is  now  in  the 
accumulation  range  while  (2.10)  of  VGB(VS)  is  valid  only  for  the  depletion  range.  To 
give  an  analytical  solution  in  the  accumulation  range,  we  use  the  asymptotic  VS(VGB) 
equation  given  by  [41]: 

Us    =    -loge[  (U^-UpB-Ug)2/^^   +    AA]  (p-Si)        (2.28) 

with  AA=  1,  then  the  HWHM  of  IB-VGB  is  given  by 

AVGB   =   AVS    +    2(kT/q)VUAAexp(|Usx_JSCR_pk|/2)  [exp(AUs/2)-l] 

(accum.    side)     (2.29) 

AVGB   =   AVS    +    2(kT/q)VUAAexp(|Usx_JSCR_pk|/2)  [1-exp  ( -AUs/2 )  ] 

(FB   side)  (2.30) 

To  compare  the  linewidth  of  JSCR  peak  to  that  of  the  MCR  peak,  we  examine 

the  change  of  surface  potential  U(YT)  with  respect  to  an  infinitesimal  change  of  VGB  by 

taking  the  derivative  of  eqn.  (2.26),  which  gives: 


AU(YT)     =     [1-(UYT/UB)1/2]AUSE    +     (UYT/UB)1/2AU 


sx 


=     [1-(YT/YEB)1/2]AUSE    +     (YT/YEB)1/2AUSX  (2.31) 

In  the  limit  that  the  interface  traps  are  located  near  the  MCR  (YT— >YEB),  the  JSCR  peak 
is  close  to  the  MCR  peak  with  gate  voltage  corresponding  to  the  weak  inversion  of  p- 
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substrate  surface  where  |AUSE|«|AUSX|.   In  the  other  limit  that  the  interface  traps  are 
near  the  drain  region  (YT— >0),  the  JSCR  peak  is  close  to  the  drain  peak  where 
|AUSE|»|AUSX|,  refer  to  Fig.  2.3.  Therefore  the  two  asymptotic  limits  of  (2.31)  are: 
AU(YT)    «     (YT/YEB)1/2    AUSX  YT->YEB,    near   MCR  (2.32) 

AU(YT)    *     [1-  (YT/YEB)1/2]    AUSE  YT->0,    near   drain  (2.33) 

The  above  equations  clearly  show  the  DCIV  linewidth  dependence  on  the  spatial 
location  of  the  interface  trap  YT.  If  the  interface  traps  are  near  the  two  edges  of  JSCR, 
the  DCIV  linewidths  approach  those  of  the  MCR  peak  and  the  drain  peak  respectively, 
as  AU(YT)=AUSX  at  YT=YEB  and  AU(YT)=AUSE  at  YT=0.  The  DCIV  linewidth 
increases  as  a  function  of  (YEB/YT)1/2  as  the  location  of  NIT  moves  away  from  the  MCR 
edge  into  the  JSCR,  it  increases  as  a  function  of  l/[l-(YT/YEB)l/2]  as  the  NIT  location 
moves  away  from  the  drain  edge. 

In  case  of  a  spatially  constant  NIT,  we  shall  evaluate  the  Ib-jscr  integral,  (2.21), 
analytically  to  describe  the  dependence  of  this  current  component  as  a  function  of  VPN. 
Since  the  integrand  is  expected  to  be  a  bell-shaped  function  with  a  maximum  surface 
recombination  rate  at  y=YM  that  satisfies  Us+J  =0  in  (2.23) ,  we  can  use  the  Taylor  series 
to  expand  the  integrand  about  y=YM  and  keep  terms  up  to  the  second  order.  The  result 
is  then 


2f2q^lcnscpsniNITZ 

• 

3(q|EM|/kT) 

[exp(UpN)     -    1] 

exp(UPN/4)  [exp(UPN/2)+cosh(UTIJ  ] 1/2 


Ib= 


(2.34) 
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EM  is  the  electric  field  at  the  position  YM(VGB).  EM  varies  with  a  low-power-exponent 
of  VPN  so  that  the  IB  dependence  on  VPN  is  dominated  by  the  exponential  dependence  of 
VPN  shown  in  numerator  and  denominator  of  (20),  which  predicts  a  n=1.33  to  n=2 
transition  from  low  forward  biases  to  high  forward  biases  with  a  transition  range  of 
about  UPN~8=200mV. 

2.4  Diffusion  Limitation  and  Other  Non-Ideal  Factors 

DCIV  measurement  on  MOS  transistors  can  use  four  different  bias 
configurations  to  inject  minority  carriers  to  the  Si02/Si  interface  [8],  as  illustrated  in 
figure  2.6.  They  are  the  drain-emitter  (DE),  source-emitter  (SE),  top-emitter  (TE)  with 
both  source  and  drain  junctions  tied  and  forward-biased  together,  and  the  bottom-emitter 
(BE)  in  transistor  structures  having  a  well/substrate  p/n  junction  which  can  be  forward 
biased.  In  each  of  the  bias  configurations,  the  p/n  junctions  not  forward  biased  are  zero- 
biased  (they  can  also  be  reverse-biased  or  even  forward-biased  at  a  lower  voltage).  In 
the  first  three  gated-emitter  configurations  (DE,  SE,  TE),  the  body-terminal  current  IB 
has  two  surface  recombination  components:  from  the  MCR,  Ib_mcR'  anc^  fr°m  me  one  or 
two  emitter-base  JSCRs,  Ib.jcr>  while  in  the  ungated  BE  configuration,  IB  contains  only 
IB.MCR  because  the  drain  and  source  p/n  junction  collectors  are  grounded  so  their 
IB_JSCR=0.  However,  the  zero-biased  source  and/or  drain  p/n  JSCR's  also  reduce  the 
interface  minority  carrier  concentrations  in  the  MCR.  This  lowers  Ib_mcr  an<^  changes 
the  IB-VGB  lineshape  which  becomes  appreciable  in  short  channels.  A  more  prominent 
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Fig.  2.6        Four  DCIV  bias  configurations:  (a)  Drain-emitter,  (b)  Source-emitter,  (c) 
Top-emitter,  and  (d)  Bottom-emitter. 
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lowering  and  shape-change  of  Ib_Mcr  arises  from  the  2-dimensional  geometrical  edge 
effect  due  to  the  proximity  of  the  grounded  collector  and  forward-biased  emitters  in  the 
thin-base  whose  base  thickness  under  the  drain  and  source  junction,  is  comparable  with 
the  metallurgical  channel  length.  Another  source  of  deviation  from  the  simple  ideal 
theory  comes  from  the  2D  impurity  doping  profiles  as  a  result  of  channel  and  junction 
engineering.  Minority  carriers  injected  from  the  body/well  junction  in  BE-DCIV  will 
have  a  drift  component  due  to  the  built-in  field  Ex(x)  =  (kT/q)d(logeNAA)/dx  from  the 
vertical  impurity  profile.  The  vertical  profile  also  gives  a  variation  of  the  injected 
minority  carrier  concentration  near  the  source  and  drain  junction  perimeters.  These 
would  give  a  difference  between  the  TE  and  BE  lineshapes.  The  channel  regions  that 
have  different  NAA(x=0,y)  and  different  oxide  charge  density,  Q0t(v)>  wiU  smft  VGB  of 
the  DCIV  peaks  which  are  combined  into  a  broader  peak. 

Diffusion  of  the  injected  carriers  in  the  BQNR  (Base  Quasi  Neutral  Region)  can 
affect  and  even  limit  the  terminal  IB,  especially  the  IB-VGB  lineshape.  For  example,  the 
interface  recombination  rate  could  be  so  high  that  the  minority  carrier  concentration  at 
the  Si02/Si  interface  is  decreased  below  the  injected  value  at  the  injecting  p/n  junction 
boundary,  Ns  <  N0  =  (n^/NAA)exp(UBE).  In  the  limit  of  NIT^°°,  Ns=(nf/NAA)  or 
UPN(x=0)=0.  The  interface  recombination  and  base  diffusion  mechanisms  are  in  series. 
In  contrast,  bulk-JSCR  recombination  (the  Sah-Noyce-Shockley  or  SNS  current 
component)  and  base  diffusion  mechanisms  are  in  parallel  [40].  Thus,  the  gate- 
controlled  terminal  IB  in  the  DCIV  methodology  is  limited  by  the  smaller  of  the  two 
series  current  mechanisms,  base  diffusion  current  at  low  VPN  and  interface 
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recombination  current  at  higher  VPN.  A  simple  analytical  solution  to  include  diffusion 
limitation  is  obtained  using  VBE  =  VPN+AVpN  and  assuming  that  AVPN  is  mainly  in  the 
BQNR  since  the  surface-space-charge  layer  is  much  thinner  than  the  BQNR.  Then, 
AVPN(VBE)  at  the  Ie-peak  can  De  obtained  as  follows  by  equating  the  interface 
recombination  current  to  the  diffusion  current: 

JB-peak    =      ^ifCnsSs^^IT/2^ 

[exp(UPN)-l]  /  [exp(UPN/2)+cosh(UTIJ  ]  (2.35) 

=       JD     =     (qDB/XB)  (nf/N^expCJpN)  [exp  (AUPN) -1  ]  (2.36) 

which  gives 

exp(AUPN)     -1 

=     (N^Xg/riiDg)  [  (cnscps)1/2NIT/2]  • 

[l-exp(-UPN)  ]/[exp(UPN/2)+cosh(UTI*)  ]  (2.37) 

=    1000.  [l-exp(-UPN)  ]/[exp(UPN/2)+cosh(UTIJ  ]  (2.38) 

=    1000./cosh(UTIJ  , 

for    0<UPN<2  [logecosh(UTI*)-2.3]=2  (  |UTI*|  -3.0)      (2.39) 
=    1000.  exp  (-UPN/2)  , 

for  UPN>2[logecosh(UTI*)+2.3]=2(  |UTI*|+1.6)  (2.40) 
DB  is  the  minority  carrier  diffusivity  in  the  base  well  of  XB  thick.  The  value  of  1000  is 
obtained  from  using  the  following  estimated  typical  values:  cns=c  =10_8cm3/s  for  a 
0.3nm  diameter  neutral  potential  well  of  the  interface  traps  and  NIT=2xl09cm-2  (which 
give  a  surface  recombination  velocity  of  S=(cnscps)1/2NIT=20cm/s);  NAA=1017cm-3  at 
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x  =  XB,  nj  =  1 0 10cm~3,  DB=10cm2/s,   and  XB  =  10-4cm  which  gives 
(NAAXB/2DBnj)=50s/cm. 

It  is  evident  from  the  approximate  results  given  in  (2.39)  and  (2.40)  that  there 
could  be  a  significant  AUPN  only  at  low  UPN's  and  only  if  the  interface  trap  energy  level 
is  near  the  mid-gap.  The  DCIV  lineshape  in  this  case  will  be  distorted  and  have  a 
diffusion-limited  flat  top.  An  interface-recombination-limited  threshold  VPN  can  be 
defined  by  equating  the  asymptotic  current  equations  limited  by  recombination  (2.36) 
and  diffusion  (2.37).  Using  the  numerical  values  just  listed,  this  gives 

uPN-th   =   2  [6.908    -    cosh(UTI*)]  (2.41) 

The  diffusion  current  JD  given  by  (2.36)  for  constant  NAA  can  be  generalized  to 
graded  and  retrograded  profiles,  NAA(x),  by  replacing  the  pre-exponential  factor  with 
the  exact  integral  given  by  j  [NAA(x)/DB(x)n^(x)]dx  which  is  integrated  over  the  quasi- 
neutral  base,  x=0  to  x=XB.  This  is  the  well-known  Gummel  Number. 

2.5  Application  Examples 

DCIV  measurements  have  been  performed  extensively  since  its  revival  in  1995 
[5]  on  silicon  micron  and  submicron  MOSTs  manufactured  by  1988-1998  technologies 
with  channel  aspect  ratios  of  W/L=100um/100n,m  to  10(im/0.1|0,m  and  gate  oxide 
thicknesses  from  30nm  to  lnm.  This  section  will  present  the  DCIV  curves  obtained  on 
both  unstressed  and  hot-carrier-stressed  MOSTs  to  illustrate  the  theoretical  analyses  of 
the  experiments.  The  data  were  from  the  large  pMOSTs  to  minimize  the  2-D  effects. 
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The  pMOSTs  have  W/L=50nm/50^m,  10.5nm  gate  oxide,  NDD(interface)  =  2xl0,7cm~3 
and  a  graded  n-well  with  p+buried  bottom  emitter. 

The  first  examples  are  from  unstressed  pMOSTs  where  the  Ib_^gb  Pea^ 
originates  from  the  residue  interface  traps  of  the  manufacturing  processing  steps.  Figure 
2.7(a)  shows  a  family  of  normalized  experimental  BE-DCIV  curves  in  solid  lines  and 
theoretical  curves  in  dashed  lines  based  on  ETI*=228meV  and  NIT=7.8xl09cm-2.  These 
numerical  values  were  obtained  from  a  2-parameter  least-squares-fit  [33]  of  the  Ig.pe^ 
vs  VPN  shown  in  Fig.  2.7(b)  which  shall  be  coined  as  the  DCIV  Peak-Plot  (P-Plot).  The 
P-plot  fits  very  well  only  the  discrete-trap-level  theory  while  the  data  deviate 
substantially  from  the  theory  computed  from  constant  or  U-shaped  interface  trap  energy 
level  distributions.  Diffusion  current  was  not  included  in  this  fit.  Fig. 7(a)  indicates  a 
good  theoretical  account  of  the  lineshape  at  the  strong  inversion  side  while  the 
experimental  halfwidth  is  wider  on  the  flatband  side.  Several  sources  contribute  to  the 
difference  including:  diffusion  limitation  at  these  low  VPN  values  which  would  broaden 
the  linewidth,  larger  surface  impurity  concentration  than  assumed  from  C-V 
measurements  and  graded  vertical  doping  profiles. 

Figure  2.8  shows  the  differences  among  the  DE-,  SE-,  TE-  and  BE-DCIV  curves 
measured  at  VPN=400mV  on  the  same  unstressed  pMOST.  The  three  gated-emitter 
configurations  (DE,  SE  and  TE)  all  show  one  narrow  peak  indicating  low  NIT 
concentrations  inside  the  gated  JSCR.  Overlap  of  SE  and  DE  curves  shows  good 
source/drain  symmetry.  The  different  VGB  ^  locations  between  TE  and  DE  or  TE  and 
SE  come  from  the  lateral  variation  of  the  minority  carrier  (holes)  concentration  at  the 
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Fig.  2.7  Experimental  BE-DCIV's  and  Peak-Plots  on  unstressed  50/50  pMOST 
compared  with  theories.  Parameter  values  are:  ETI*=228meV  and 
NIT=7.8xl09cm~2  from  least-squares-fit  to  the  discrete-level  theory  using 
cns=c  =10_8cm3/s  (assumed),  n^lO^cm-3,  and  from  Cg-VGB  data, 
NAA(x=0)=2xl017cnT3,  and  Xox=10.5nm.  (a)  Peak  Plots  (IB.Peak  vs  vpn)  of 
experimental  data  (points)  and  theory  from  three  assumed  DOS's  of 
interface  traps:  a:  Discrete-level  given  above,  b:  constant  DOS  over  EG.Si 
with  DIT  =  7.  1  x  1  09/cm2-e  V,  and  c:  U-shaped  DIT  = 
2.5xl09cosh(ETI/200meV)/cm2-eV.  (b)  Lineshape  comparison  with  one- 
discrete-level  theory. 
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interface,  P$(y),  from  Ps(0)=0  at  the  shorted  source  junction  to  Ps(L)=Pss[exp(UPN)-l] 
at  the  forward-biased  drain  junction  in  the  DE-DCIV  bias  configurations  or  vise  versa  in 
the  SE-DCIV  bias  configuration.  Ps(y)  is  further  reduced  due  to  the  finite  diffusion 
length  towards  the  MCR  from  the  injection  edge  of  the  drain  or/and  source  junction. 
These  give  a  smaller  effective  channel  interface  recombination  area  and  also  a  smaller 
effective  VPN  which  would  shift  the  peak  location  towards  the  -VGB  side  as  indicated  by 
(2.9).  Figure  2.8  shows  that  the  Ib"Vgb  linewidth  of  the  ungated  emitter  (BE)  is  almost 
doubled  in  comparison  with  those  of  the  gated  emitter  (SE,  DE  and  TE),  which  is  due  to 
carrier  injection  against  a  diffusion  barrier  from  graded  BQNR.  These  effects  on  the 
lineshape  from  junction  and  geometrical  edge  effects  and  impurity  concentration  profile, 
indeed,  provide  flexible  variables  that  can  be  leveraged  in  the  monitor  application  of  the 
DCIV  methodology  for  optimization  of  transistor  designs  and  fabrication  processes. 

For  the  second  sets  of  application  examples,  the  DCIV  theory  is  used  to  evaluate 
the  interface  traps  generated  during  the  substrate-hot-carrier  (SHC)  and  channel-hot- 
carrier  (CHC)  stresses  of  the  same  50/50  pMOSTs  used  in  the  preceding  unstressed 
examples.  During  the  SHC  stress,  the  p+buried/n-well  junction  is  forward  biased  to 
inject  a  high  concentration  of  holes  into  the  strongly  inverted  hole  surface  channel 
which  is  reverse  biased  by  a  high  reverse  voltage  applied  simultaneously  to  the  drain 
and  source  p+/n-well  junctions.  This  high  reverse  voltage  accelerates  the  injected  holes 
to  a  kinetic  energy  greater  than  the  Si:H  interfacial  bond  energy  (3.1eV)  in  order  to 
release  the  hydrogen  and  generate  the  silicon  dangle  bonds,  Si»,  which  are  the  interface 
traps.  Reasonably  good  lateral  uniformity  of  interface  trap  generation  was  expected 
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from  using  low-enough  VBE  and  strong-enough  inversion  VGB  to  minimize  the  lateral 
drop  of  the  vertical  acceleration  potential  due  to  channel  ohmic  resistance.  Figure  2.9(a) 
shows  the  IB-VGB  DCIV  lineshapes  which  indicates  a  larger  experimental  linewidth 
suggesting  diffusion  limitation.  Figure  2.9(b)  shows  the  DCIV  Peak-Plot,  le-peak'^PN' 
which  includes  also  the  unstressed  data  that  coincide  the  stressed  data  for  VPN<400mV 
further  indicating  diffusion  limitation.  The  stress-generated  interface  traps  are  still  quite 
discrete  as  indicated  by  the  good  fit  to  the  discrete-energy-level  model  at  high  VpN's. 
The  least-squares-fit  (LSF)  stressed-generated  trap  energy  level,  E^*  =  280meV,  should 
shift  towards  midgap  when  diffusion  is  taken  into  account. 

A  second  example  of  stress-generated  interface  traps  was  obtained  using  the 
channel-hot-carrier  stress  on  the  same  50/50  pMOST  whose  DCIV  at  VPN=350mV  and 
P-Plot  are  shown  in  Fig.  2.10.  The  unstressed  peak  near  the  intrinsic  gate  voltage 
(VBG=+0.1V)  does  not  grow  at  all  which  is  consistent  with  its  MCR  location  where 
there  are  few  if  any  hot  carriers.  The  growing  post-stress  peak  is  broad  (FWHM  -  2V) 
and  in  the  accumulation  range  (VGB=+1.8V).  It  has  an  n— >2  over  many  decades  in  the 
P-Plot  (inset  of  Fig. 8)  which  suggests  that  the  interface  traps  are  mainly  generated  in  the 
drain  JSCR,  as  anticipated  by  the  JSCR-DCIV  theory  just  described,  and  the  CHC  stress 
theory  [42]  where  the  channel  hot  holes  and  the  secondary  hot  electrons  are  generated  in 
the  drain  junction  space  charge  region  in  the  pMOST  during  stress. 
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2.6  Summary 

Theory  and  experiments  on  gate-voltage  modulation  of  the  recombination  rate  at 
the  interfacial  electronic  traps  are  quantitatively  correlated  in  prestress  and  two 
poststress  (channel  and  substrate  hot  carrier  stresses)  DCIV  applications.  The  wider 
experimental  linewidths  could  be  an  indication  of  diffusion  limitation  in  the  low  forward 
bias  range,  and  from  2-dimensional  effects  at  high  forward  bias  range  due  to  lateral 
variations  of  the  interface  minority  carrier  concentration  along  the  gated  conduction 
channel.  The  examples  illustrate  the  identification  of  discrete  interface  trap  energy 
levels  and  the  near  midgap  energy  level  positions. 


CHAPTER  3 

SPATIAL  DENSITY  PROFILE  OF  INTERFACE  TRAPS 

FROM  DCIV  MEASUREMENTS 

3.1  Introduction 

The  DC  basewell  terminal  (or  pad)  current,  IB,  modulated  by  the  applied 
gate/base  DC  voltage,  VGB,  in  MOS  transistors  was  a  method  recently  reactivated  to 
monitor  electric-field-stress  generated  interface  traps  as  a  transistor  reliability  monitor 
[5]  and  to  serve  as  a  pre-stress  diagnostic  monitor  for  transistor  design  and  processing 
[8].  The  VGB-modulated  IB  arises  from  recombination  of  the  injected  minority  carriers 
by  a  forward-biased  p/n  junction  (Drain/Base,  Source/Base,  or  Substrate/Basewell)  with 
the  majority  carriers  in  the  basewell,  at  the  Si02/Si  interface  traps  under  the  gate  oxide. 
Electron-hole  recombination  can  occur  at  interface  traps  under  the  gate  oxide  in  three 
regions  [3,4]:  (1)  the  mid-channel  region  MCR  (2)  the  drain  and  source  junction  space- 
charge  regions  D-JSCR  and  S-JSCR,  and  (3)  drain  and  source  lowly-doped  extension 
regions  D-LDER  and  S-LDER.  A  detailed  theoretical  analysis  was  presented  in  chapter 
2  for  the  MCR  and  JSCRs.  Some  of  our  recent  application  demonstrations  [8,33,34] 
focused  on  interface  traps  NIT  in  the  MCR.  However,  recombination  at  interface  traps 
in  the  other  two  regions,  JSCR  and  LDERs,  becomes  increasingly  important  in 
unstressed  transistors  as  the  channel  length  is  scaled  down  and  it  is  well-known  that 
recombination  in  JSCR  dominates  in  stressed  transistors  [8,43,44] 
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regardless  of  the  channel  length  [8,33].  It  was  anticipated  [8]  and  analytically 
demonstrated  in  the  previous  chapter  that  the  Ib_Vgb  Hneshape  and  its  VPN 
dependencies  are  affected  by  two  lateral  variations:  the  interface  traps,  NIT(y),  [8]  and 
the  interface  concentration  of  the  injected  minority  carriers  N(x=0,y)  [8]  in  n-channel 
MOSTs  or  P(x=0,y)  in  p-channel  MOSTs.  This  spatial  dependency  was  analytically 
shown  in  chapter  3  to  be  substantial  in  the  JSCR  due  to  the  rapid  spatial  variation  of  the 
carrier  concentrations  in  the  JSCR.  We  will  now  describe  the  procedures  for  obtaining 
the  NIT(y)  distribution  and  some  transistor  design  parameters  (gate  oxide  thicknesses 
and  dopant  surface  concentration).  The  theoretical  analyses  and  computed  illustrations 
are  given  in  section  3.2.  Analyses  of  the  experimental  data  to  give  transistor  design 
parameters  in  the  MCR  and  NIT(y)  profiles  in  both  the  JSCR  and  MCR  are  presented  in 
sections  3.3,  3.4  and  3.5  respectively.  A  summary  is  given  in  section  3.6. 

3.2  Theoretical  Analysis 

The  analytical  formulas  were  derived  and  described  in  chapter  3  for  the 
Shockley-Read-Hall  (SRH)  steady-state  recombination  rate  Rss  at  interface  traps  in  the 
MCR  and  JSCRs.  There  is  an  additional  VGB-dependence  of  Rss  from  NIT  in  the 
LDERs,  which  will  be  included  in  this  chapter.  The  present  analyses  for  the  LDER 
follow  the  steps  for  the  MCR  with  two  noticeable  differences:  (i)  the  nearly  zero 
flatband  voltage  in  the  LDERs  due  to  the  same  dopant  type  and  nearly  the  same  high 
dopant  impurity  concentrations  in  the  Si-gate  and  drain  and  source  LDERs  according  to 
the  traditional  designs  of  digital  CMOS  transistors  and  (ii)  the  gate-to-drain  voltage  is 
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now  VGD  =  VGB-VDB  =  VGB±VPN  (+  for  nMOST  and  -  for  pMOST)  which  gives  IB  the 
additional  VPN  dependency,  both  are  leverageable  for  NIT  characterization  described  as 
follows. 

As  illustrated  by  the  computed  potential  variation  curves  in  Fig.l,  the  (VPN) 
forward-bias-injected  minority  carriers  reduce  the  interfacial  potential  barrier  height  on 
both  the  MCR-side  (lower  curves)  and  the  LDER-side  (upper  curves)  of  the  JSCR 
shown  in  Fig.  3.1(a),  and  the  total  potential  barrier  height  through  the  JSCR  shown  in 
Fig.  3.1(b).  This  VPN  dependency  significantly  affects  the  DCIV  or  IB-VGB  lineshape 
from  NIT's  located  in  the  JSCR.  [The  interfacial  electric  potential,  V(x,y)  at  the  Si02/Si 
interface  x=0,  has  been  denoted  by  the  two  traditional  and  two  recent  engineering 
notations:  \j/(x=0,y)  =  \|/s(y)  =  Vj(x=0,y)  =  VIS(y).  The  subscript  I  for  intrinsic  Fermi 
potential  will  be  dropped  for  arbitrariness  of  the  reference  electric  potential  and  for 
notation  compactness  using  kT/q  normalization:  Us(y)  =  (q/kT)VIS(y).]  The  curves  in 
Figs.  3.1(a)  and  (b)  were  computed  using  the  carrier  depletion  as  the  zeroth 
approximation  in  the  JSCR  [45]  but  they  do  include,  in  the  MCR  and  the  LDERs,  the 
injected  carrier  concentration  from  the  forward  bias  (VPN)  just  described,  and  also  the 
Fermi-Dirac  carrier  distribution  and  impurity  deionization  for  high  impurity  surface  or 
interface  concentrations  [46].  The  upper  part  of  Fig.  3.1(a)  gives  the  total  lateral  drop  of 
the  interfacial  potential  in  the  LDER  with  10~19cm  3.  The  potential-drop  is  about  100 
times  smaller  than  that  in  the  MCR  with  1017cm-3  shown  in  the  lower  part  of  Fig.  3.1(a) 
as  expected  from  the  concentration  ratio.  Figure  3.1(b)  shows  that  at  VGB=-TV  (about 
flatband  in  the  MCR),  the  p/n  junction  barrier  height,  VB,  at  the  interface  is  reduced 
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Fig.  3.1  Gate-Basewell  voltage  (VGB)  dependence  of  the  potential  variation  through 
the  space-charge  layer  at  the  Si02/Si  interface  of  the  drain  or  source  n+/p- 
Base  junction  covered  with  a  n+Si  gate  in  the  carrier-depletion  and  thick- 
oxide  limits  with  spatially  constant  dopant  impurity  concentrations  of 
Nnn=lxl019cm~3,  and  PRR=lxl017cm~3,  and  oxide 
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thickness  of  Xox=3.0nm.  (a)  The  potential  drop  in  the  drain  and  basewell 
sides  of  the  junction  space-charge  region,  VD.JSCR  and  VB.JSCR  for  forward 
applied  dc  biases  of  VpN=0mV  to  800mV.  (b)  The  potential  variation  along 
the  SiQ2/Si  interface  for  VPN=0mV  and  800mV. 
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Fig.  3.2  The  length  of  junction  space-charge  region  as  a  function  of  the  applied  dc 
gate  voltage,  VGB,  for  applied  junction  dc  biases  VPN=0mV  to  800mV.  (a) 
Total  length,  YJSCR.  (b)  Length  on  the  basewell  and  drain  (or  source)  sides, 
YjscR-Base  anc*  YjSRC-Drairr  Transistor  parameters  are  the  same  as  those 
used  in  Fig.  3.1. 
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nearly  600mV  by  the  injected  carriers,  from  VB=720mV  at  VPN=OmV  to  VB=120mV  at 
VPN=800mV.  The  geometrical  effect,  shortening  of  the  length  of  the  JSCR  by  VPN  in 
Fig.  3.1(b),  should  also  be  noted  which  is  leveraged  in  NIT  profiling.  This  is  further 
illustrated  in  Fig.  3.2(a)  for  the  total  length,  YJSCR,  and  in  Fig.  3.2(b)  respectively  for 
YJSCR's  component  on  the  Channel-side  (left  Y-axis)  and  the  LDER-  (or  Drain-  and 
Source-)  sides  (right  Y-axis)  whose  curves  are  nearly  the  same  with  ratio  100  from  the 
dopant  impurity  concentration  ratio  at  the  interface,  1019/1017=100. 

The  effects  of  this  lateral  interfacial  potential  variation  with  forward  bias  are 
shown  in  Figs.  3.3(a)  and  3.3(b)  in  which  the  normalized  SRH  recombination  rate, 
Rss/Rss  k,  or  normalized  basewell-pad  current,  Ib/IbdIc  at  a  midgap  interface  trap  in  the 
JSCR  are  computed  as  a  function  of  VGB.  As  described  in  the  following  two 
paragraphs,  the  Ib"Vgb  lineshape  provides  the  basis  for  calculating  the  spatial 
dependence  of  the  interface  trap  density,  NIT(y)  from  experimental  DCIV  Ib_Vgb 
curves. 

In  Fig.  3.3(a),  the  effects  of  spatial  variation  of  NIT(y)  in  the  JSCR  on  the  Ib-Vgb 
characteristics  (at  VPN=200mV)  are  illustrated  for  the  seven  NIT(y)  distributions: 
NIT(y)=NIT0=constant  in  the  (i)  MCR,  (ii)  Drain  (or  Source)  LDER,  and  (iii)  Drain  (or 
Source)  JSCR;  and  the  four  discrete  N1T(y)=NIT18(y-YT)  distributions  located  at  (iv)- 
(vii)  y=YT=10,  20,  30  and  40nm  from  the  metallurgical  p/n  junction  boundary.  There 
are  several  obvious  features  in  Fig.  3.3(a)  useful  for  NIT(y)  characterization  from 
experimental  data.  First  consider  the  six  distributions  except  the  constant-NIT  in  the 
Drain-  or  Source -JSCR,  (iii).  (1)  The  VGB  k  moves  from  inversion  to  depletion- 
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Fig. 3  Effects  of  spatial  variation  of  interface  trap  concentration  (NIT)  and  forward 

bias  voltage  (VPN)  on  the  normalized  DCIV  lineshape  (le-Normaiize  or 
^ss-Normalize  vs  ^gb)  from  an  Si-midgap  (ETI*=0mV)  interface  trap  located 
in  the  drain  and  source  lowly-doped-extension  regions  (D-LDER  and  S- 
LDER),  the  junction  space-charge-region  (JSCR),  and  the  mid-channel- 
region  (MCR).  (a)  Three  spatially  constant  NIT's:  curvel=MCR, 
curve2=D-LDER  and  S-LDER,  and  curve3=JSCR,  and  four  discrete  (delta- 
function)  NIT's  located  at  YT=10,  20,  30,  and  40nm  (curve4,5,6,7)  from  the 
p/n  impurity  concentration  boundary,  at  a  forward-bias  of  VPN=200mV.  (b) 
Forward-bias  dependence  (VPN=0mV  to  800mV)  for  the  constant  N,T  in  the 
JSCR. 
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flatband-accumulation  as  the  discrete  NIT(y=YT)  moves  from  the  MCR  edge  towards  the 
p/n  boundary  y=0.  (2)  The  IB-VGB  lineshape  from  the  constant-NIT  in  the  MCR  and  D- 
LDER  (or  S-LDER)  and  from  the  four  spatially  discrete  NIT(y)  all  are  sharply  peaked 
Gaussian-like  with  increasing  width  as  YT  moves  from  the  MCR  towards  the  Drain  (or 
Source)  LDER.  These  two  lineshape  properties  provide  the  basis  for  its  utility  as  a 
delta-function  probe  for  the  lateral  variation  of  NIT(y). 

For  the  spatially  constant  NIT(y)=NIT0  in  the  JSCR,  (iii),  the  following  effects  on 
the  lineshape  and  position  are  noted.  From  the  VPN=200mV  curve  in  Fig.  3.3(a):  (31) 
the  VGBpk  is  near  that  of  the  NIT  located  in  the  MCR  (i),  and  (32)  the  upper  part  of  the 
lineshape  is  similar  to  that  of  MCR-NIT  (i)  also.  The  curves  in  Fig.  3.3(b)  shows  the 
effect  of  increasing  VPN  or  injected  minority  carrier  concentration:  (33)  VGB  k  shifts 
toward  flatband-accumulation  (negative  VGB  in  nMOST),  (34)  the  lineshape  broadens, 
(35)  a  substantial  plateau,  about  25%  for  VPN— >0,  diminishes,  (36)  the  upper  20%  of  the 
IB-VGB  lineshape  is  nearly  Gaussian  with  a  slightly  skew  on  the  accumulation  side,  and 
(37)  the  entire  half  of  the  IB-VGB  lineshape  on  the  inversion-side  (positive  VGB  in 
nMOST)  is  nearly  Gaussian-like. 

Figures  3.4  and  3.5  provide  theoretical  estimates  of  the  confidence  level  or  the 
accuracy  of  using  the  delta-function  approach  to  analyze  the  DCIV  IB-VGB  data  for 
extracting  slowly  and  rapidly  varying  NIT(y)  profiles  respectively.  The  assumed 
(continuous  lines)  and  the  extracted  (crosses)  profiles  are  shown  in  Figs.  3.4(a)  and 
3.5(a),  respectively  for  the  exponential  profile,  exp(-y/3nm),  and  the  Gaussian  profiles, 
exp[-(y-YT)2/Ly]  located  at  four  distances  Y-^lO^O^O^Onm  from  the  p/n  boundary. 
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Fig.  3.4  Illustration  of  the  use  of  the  DCIV  line-shape  for  profile  extrapolation  of  the 
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Fig.  3.5  Illustration  of  the  use  of  the  DCIV  line-shape  for  profile  extrapolation  of  the 
mid-gap  interface  trap  concentration,  NIT(y),  for  four  Gaussian  profiles 
located  at  YT=10,  20,  30  and  40nm  from  the  n+D/p-B  boundary  inside  of 
the  JSCR  at  VPN=200mV.  Transistor  parameters  are  the  same  as  those  of 
Fig.  3.1.  (a)  Assumed  and  extrapolated  profiles,  (b)  Fractional  errors,  (c) 
Current  integrals  of  the  profile-extrapolation  algorithm. 
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The  difference  or  fractional  error  of  the  extracted  profiles  are  shown  in  Figs.  3.4(b)  and 
3.5(b),  indicating  excellent  accuracy  (<0.07  or  7%)  for  the  slowly  varying  exponential 
N|T  in  Fig.  3.4(b),  but  rather  large  differences  for  the  rapidly  varying  Gaussian  NIT  in 
Fig.  3.5(b),  the  latter  as  expected  since  Rssi^'^gb)  is  ^ar  from  a  delta-function, 
nevertheless,  Fig.  3.5(a)  shows  that  extrapolated  profile  are  peaked  and  Gaussian-like 
with  a  peak  location  near  the  true  physical  location  and  with  a  linewidth  about  1.5-times 
wider  which  increases  ~2-times  wider  as  the  peak  moves  towards  the  base-side  or 
channel-side  of  the  JSCR.  Thus,  theoretical  corrections  could  be  made  to  correct  part  of 
such  a  systematic  difference. 

The  NIT(y)  profile  extraction  procedure,  used  to  give  Figs.  3.4(a),  3.4(b),  3.5(a) 
and  3.5(b),  starts  with  the  delta-approximation  for  evaluating  the  general  formula  (2.21) 


IR-.T.qrR(VnR)     =    qW 


YJSCR 

NIT(y)  .Rssl(y,VGB)dy  (3.1) 

0 


where  W  is  the  channel  width.  The  integrand  Rssi(v'Vgb)  *s  me  SRH  recombination 
rate  per  unit  interface  trap.  Its  lineshape  in  y  and  VGB  are  both  bell-shaped  and  rather 
sharp,  resulting  in  a  simple  but  accurate  mapping  of  the  experimental  Ib~Vgb  to  the 
N,T(y)-y  profile.  The  bellshape  characteristics  are  illustrated  by  the  following  examples. 
The  maximum  of  Rssi(y)-y  at  a  given  VGB  [See  Fig.  3.6(a).]  occurs  at  the  position 
y=^M(^GB)  which  was  given  by  Eq.  (2.4): 

us*j   =    U(YM)-(Up+UN)/2    +    loge(cns/cps)1/2    =    0  (3.2) 

Up  and  UN  are  the  hole  and  electron  quasi-Fermi  potentials,  respectively.  Using  the 
U(y)  from  the  depletion  approximation  and  the  abrupt  dopant  impurity  concentration 
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Fig.  3.6  Gate-voltage  dependence  of  the  spatial  variation  of  the  steady-state 
recombination  rate  at  VPN=200mV.  Other  parameters  are  the  same  as  those 
of  Fig.  3.1.  (a)  Gate  voltage  dependence  of  RSS1  vs  y.  (b)  Location,  YM,  of 
the  peak  RSS1  vs  gate  voltage. 
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model,  YM  has  the  following  dependence  on  VGB  where  Usx  and  USE  are  total  surface 
potential  changes  in  the  x-direction  normal  to  the  Si02/Si  interface  in  the  MCR 
(Basewell)  and  the  D-LDER  or  S-LDER  (Emitters)  respectively. 

YM(VGB)    =    {7ld    Hubi-Usx-Upn+Use    "    ^luF-UPN/2-Usx-loge(cns/cps)1/2] 

(3.3) 

LD  is  the  Debye  length,  UBI  the  built-in  potential  of  the  n+D/p-B  junction.  This 
equation  maps  a  gate  voltage  to  a  spatial  location  in  the  JSCR  in  order  to  extract  the 
NIT(y)  profile  from  the  experimental  IB-VGB.  Shown  in  Fig.  3.6(a)  for  the  nMOST  with 
midgap  trap  at  VPN=200mV,  as  VGB  is  swept  from  subthreshold-inversion  (VGB=-0.6V) 
to  flatband-accumulation  (VGB=-1.4V),  the  position  of  the  RSS1  maximum,  YM,  moves 
from  the  MCR-edge  of  the  JSCR  (YM=46nm)  towards  the  metallurgical  p/n  boundary  or 
the  LDER-edge  of  the  JSCR  (YM=9nm  -»  0).  This  shift  is  computed  from  (3.3)  and 
shown  in  Fig.  3.6(b)  indicating  almost  a  linear  dependence  of  30nm/V  in  the  flatband- 
accumulation  range  (VGB<-0.94V).  Another  feature  is  the  narrowing  of  the  Rns(y) 
linewidth  towards  flatband-accumulation  shown  in  Fig.  3.6(b)  which  sharpens  the  probe 
for  the  NIT(y)  profile.  The  sharp  bellshaped  Rssi(y)  can  De  approximated  by  a  delta- 
function  to  evaluate  the  IB  integral,  (3.1): 


ib-jscr(vgb>    -   QWNIT(YM) 


YJSCR 

Rssi<y<VGB>dy  =   NIT(YM)*IB_1(VGB)  (3.4) 

0 


where  Ib_jscr  *s  trie  experimental  base-terminal  current  or  theoretically  generated 
current  using  (3.1).  The  interface  trap  profiling  formula  in  the  JSCR  is  then  given  by 

NIT(YM)     =    IB-jscr<Vgb)/Ib-i(Vgb)  (3.5) 
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where  IB.i(VGB)  defined  by  (3.4)  is  the  recombination  current  from  a  spatially  constant 
interface  trap  concentration  inside  the  JSCR  normalized  to  NIT(0<y<YM)=l. 

To  estimate  the  error  in  the  extracted  NIT  profile  using  the  delta-function 
approximation  in  (3.4),  an  assumed  profile  of  NIT  in  the  JSCR  is  used  to  compute 
IB.j(VGB)  from  (3.1)  and  IB.i(VGB)  from  (3.4).  These  were  given  in  Fig.  3.4(c)  for  the 
exponential  profile  and  Fig.  3.5(c)  for  the  four  Gaussian  profiles.  The  ratio  of  these  two 
curves  then  give  the  extracted  NIT  profiles  in  Figs.  3.4(a)  and  3.5(a)  and  the 
corresponding  difference  errors  in  Figs.  3.4(b)  and  3.5(b).  Practical  applications  using 
this  procedure  to  extract  pre-stress  and  post-stress  NIT  profiles  are  given  in  section  3.4. 

Figures  3.7-3.9  give  additional  illustrations  of  the  dependencies  of  the  DCIV  IB- 
VGB  characteristics  on  the  forward-bias  VPN  or  injected  minority  carrier  concentration. 
These  VpN  dependencies  can  be  used  to  determine  additional  properties  of  the  interface 
traps.  The  trap-spatial-location  effects  on  the  VGB  k-VPN  is  shown  in  Fig.  3.7  for  the 
NIT  distributions  given  in  Fig.  3.3(a)  which  indicate  that  increasing  forward  bias,  VPN, 
or  minority  carrier  concentration  would  decrease  VGB  k  towards  flatband-accumulation 
if  NIT(y)  is  constant  in  the  JSCR  and  MCR,  or  localized  near  MCR  (YT>30nm)  while 
increase  VGB  k  towards  subthreshold  and  inversion  if  NIT(y)  is  localized  (delta-function) 
in  the  JSCR  (YT<30nm)  or  in  the  LDER. 

The  trap-energy-location  effects  on  the  shape  of  the  Rss  or  IB  vs  VGB  curve  is 
shown  in  Fig.  3.8(a)  for  a  spatially  constant  NIT  in  the  JSCR  and  four  trap  energy  levels, 
^,=0.0  V  (midgap),  0.1V,  0.2V  and  0.3V.  The  shallower  traps  (E^O^V,  0.3V)  have  a 
broader  Rss'^gb  or  ^b'^gb  an<^  steePer  and  higher  plateau  in  the  flatband-accumulation 
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Fig.  3.7        Forward-bias  (VPN)  dependence  of  the  VGB  at  IBpk  for  the  seven  NIT(y) 
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Fig.  3.8  Effects  of  energy  level  of  the  interface  trap  in  the  JSCR  on  the  DCIV 
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in  Fig.  3.1.  (a)  Constant  NIT  in  the  JSCR.  (b)  Discrete  NIT  at  YT=20nm 
from  the  n+D/p-B  boundary. 
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range  (VGB<-1  V  in  this  nMOST  example).   Figure  3.8(b)  is  for  localized  NIT  in  the 

JSCR  which  shows  broadened  IB-VGB  when  the  energy  level  of  the  interface  trap  moves 

from  the  midgap  EIT  =  0  towards  one  of  the  band  edge,  EIT  =  0.3eV. 

The  trap-energy-location  effects  on  the  shape  and  the  normalized  reciprocal 

slope  of  the  log(IB  k)  vs  VPN  characteristics  are  shown  in  Figs.  3.9(a)  and  3.9(b).  (We 

called  this  the  Peak  Plot  or  PP  in  contrast  to  the  Gummel  Plot  and  the  Shockley  Plots.) 

This  PP  is  probably  the  most  powerful  result  of  the  DCIV  methodology  for  rapid 

interface  trap  energy  level  determination  because  the  VPN-dependence  of  the  peak  IB 

current  is  a  function  of  only  the  effective  interface  trap  energy  level,  ETI*  =  ETI  + 

(kT)loge(cns/c    )1/2.  It  is  completely  independent  of  the  surface  energy  band  bending 

and  other  material  properties,  such  as  surface  impurity  concentrations,  oxide  thickness, 

and  carrier  capture  rates.  Although  ETI*  contains  the  carrier  capture  rate  ratio  it  would 

only  give  an  uncertainty  in  the  extracted  ETI  of  (kT/q)loge(cns/c    )1/2  ~  l.l(kT/q)  ~ 

30mV  since  the  carrier  capture  ratio  is  usually  in  the  range  of  0.1  to  10  for  neutral  trap 

potential  wells  of  the  intrinsic  and  silicon-dangling-bond  interface  traps  observed 

experimentally  thus  far.  Figure  3.9(a)  gives  the  IBpk-VPN  plots  for  four  Ey^eV)  (0,  0.1, 

0.2,  0.3)  with  the  seven  spatial  NIT  distributions  of  Fig.  3.3(a).  Figure  3.9(b)  gives  the 

reciprocal  slope  defined  by 

exp(UPN/n)-l    =    [exp(UPN)-l]/[exp(UPN/2)+cosh(UTIJ]  (3.6) 

Six  of  the  NIT  distributions  have  identical  PP's  (unbroken  lines)  in  Fig.  3.9(a)  which  are 

the  two  constant-NIT's  in  the  MCR  and  LDERs  (Drain  and  Source),  and  the  four  delta- 

function-NIT's  in  the  JSCR.  Their  n-factor,  shown  Fig.  3.9(b),  approach  the  ideal 
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Fig.  3.9  Effects  of  energy  level  of  the  interface  trap  in  the  JSCR  and  MCR  on  the 
peak  recombination  current  vs  forward  bias.  Other  transistor  parameters  are 
the  same  as  those  in  Fig.  3.1.  (a)IB  k  and  KBpk  vs  VPN.  (b)  Reciprocal 
slope  or  voltage-swing  vs  VPN. 
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Shockley  value  of  n=l  at  low  VPN  for  the  three  non-midgap  traps  (ETI*  =  O.leV,  0.2eV, 
0.3eV),  increasing  to  n=2  at  high  VPN  or  P§=NS  »  n;  in  a  small  range  of  VPN  about 
200mV  (10%  to  90%  rise).  It  degenerates  to  n=2.0  at  all  VPN  for  the  midgap  trap 
Ej^O.  The  seventh,  constant-NIT  in  the  JSCR,  has  a  slightly  higher  PP  (broken  line)  in 
the  lower  VPN  range  and  a  higher  n  than  the  other  six  distributions.  For  this  case,  Fig. 
3.9(b)  also  illustrates  the  contribution  to  IB  k  from  the  NIT  at  the  channel  or  MCR-edge 
of  the  JSCR.  This  geometric  effect  from  YJSCR(VPN)  increases  n  above  2  at  low  VPN  for 
the  midgap  trap  as  indicated  by  the  dotted  line  in  Fig.  3.9(b).  To  remove  this  geometric 
effect,  the  YJSCR(VPN)  dependence  can  be  taken  out  approximately  by  normalizing  the 
numerically  integrated  IB  k  integral  to  YJSCR(VpN)  before  using  (3.6)  to  compute  n. 
Such  a  procedure  is  desired  in  this  theoretical  illustration  since  the  reciprocal  slope  is 
defined  for  the  purpose  of  quickly  assess  the  properties  of  the  interface  traps  from  the 
experimental  VPN  dependence  of  IB  k  used  in  the  initial  analyses  of  the  gate-controlled 
base-terminal  current  [5,8]  and  the  1957  SNS  theory  for  the  recombination  at  traps  in 
the  bulk  JSCR  [40]. 

The  insensitivity  on  the  spatial  distribution  of  NIT  of  the  six  distributions  first 
considered  above  allows  an  accurate  determination  of  the  effective  energy  level  of  the 
interface  traps,  E^*,  and  the  thermal  activation  energy  of  the  IB  which  is  an  important 
component  of  the  standby  current  in  CMOS  operations.  In  practice,  the  interface  traps 
have  always  been  observed  near  the  midgap  (ETI*=0eV).  Some  extracted  shallow-trap 
energy  levels  in  previous  analyses  of  data  [8,33,34]  originated  from  emitter-injection- 
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limited  and  diffusion-limited  IB  current  through  a  thick  basewell  layer  or  long  channels 
rather  than  from  shallow  interface  traps. 

3.3  Extraction  of  MOST  Device  Parameters 

One  of  the  applications  of  the  DCIV  base-recombination-current  methodology 
was  its  use  as  the  diagnosis  monitor  for  transistor  design  and  fabrication  processes  [8]. 
In  this  section,  a  demonstration  is  given  on  the  extraction  of  the  surface  concentration  of 
the  base  dopant  impurities,  Nxx,  and  the  oxide  thickness,  Xox,  and  their  comparisons 
with  the  design  values.  A  sample  of  the  experimental  Top-Emitter  DCIV  (TE-DCIV) 
data  matrix  used  in  this  demonstration  is  shown  in  Figs.  3.10(a)  and  (b)  for  two  drawn 
channel  lengths,  W/L=10|Lim/10|lim  and  10u,m/0.25^im,  of  pre-stressed  pMOSTs  on  8- 
inch  wafers  fabricated  by  a  0.25u,m  technology.  The  other  4-sets  in  the  6-set  matrix 
were  obtained  for  L  =  5,  2,  0.8  and  0.4u.m  drawn  channel  lengths.  Similar  application 
results  have  been  obtained  for  earlier  (0.5u,m  and  0.35jxm)  and  later  (<0.18u;m) 
technologies  [8]. 

The  broke  lines  in  Figs.  3.10(a)  and  (b)  show  the  location  of  the  (iBpk'^GBpk) 
whose  shift  towards  flatband-accumulation  (positive  VGB  for  pMOSTs)  with  increasing 
drain  and  source  forward-bias,  VPN=0-700mV,  is  characteristic  of  an  interface  traps  in 
the  MCR.  The  IB-VGB  curve  around  the  IB  k  is  least-squares-fitted  (LSF)  to  a  Gaussian 
with  a  baseline, 

IB    =    Ai+   A2.exp{-[  (VGB-A3)/A4]2}  (3.7) 
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Fig.  3.10  Experimental  TE-DCIV  of  pMOSTs  from  quarter-micron  CMOS 
technology  for  VPN=50mV  to  700mV.  Channel  Width  to  Length  aspect 
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The  four  fitting  parameters  are  the  baseline  Aj=IB0,  the  peak  current  A2=IBpk,  the  peak 
gate  voltage  A3=VGB  k  and  the  linewidth  A4.  These  four  LSF  parameters  were  obtained 
for  the  84  DCIV  curves  from  pMOSTs  with  the  six  L  at  the  fourteen  VPN  values  (50, 
100,  ...,  650,  700mV).  Figure  3.11  shows  a  typical  LSF  result  for  the  L=10u,m  and 
0.25|im  pMOST  at  VPN=400mV.  The  simple  Gaussian  cannot  be  used  to  fit  the  entire 
gate  voltage  range  because  (i)  the  DCIV  lineshape  from  NIT  in  MCR  is  skewed  toward 
the  flatband-accumulation  side  of  the  peak  (towards  positive  VGB  in  pMOSTs),  with  a 
slightly  wider  half-width-half-maximum  (HWHM)  as  previously  shown  chapter  3  and 
the  skew  was  also  observed  in  the  analytical  solution  shown  in  Fig.  3.3(b)  for  a  constant 
NIT  in  the  JSCR  and  (ii)  constant  or  peaked  NIT  in  the  JSCR  contributes  to  additional  IB 
to  that  from  NIT  in  the  MCR  thereby  lineshape  widening  on  the  flatband-accumulation 
side  of  the  MCR  peak.  To  circumvent  these  two  asymmetry  possibilities,  only  data 
points  from  the  inversion  range  and  a  small  range  on  the  accumulation  side  of  the  IB  k 
are  LS  Fitted  to  (3.6).  These  are  the  open  circles  (L=10u,m)  and  squares  (L=0.25u.m) 
shown  Fig.  3.1 1.  The  difference,  Ib.data  ~  ^B-FIT'  are  snown  by  tne  tw0  lower-height 
humps  labeled  JSCR  in  Fig.  3.1 1  which  are  from  the  NIT  distributions  in  the  drain  and 
source  JSCRs  or  LDERs.  The  VGB  k  (or  A3)  from  the  LSFs  are  plotted  as  a  function  of 
VPN  in  Figs.  3.12(a)  and  3.12(b)  for  the  main  IBpkl.  The  increasing  VGB  kl  towards 
flatband-accumulation  with  increasing  VPN  (pMOST)  suggest  that  the  IBpkl  arises  from 
NIT  located  in  the  MCR.  Thus,  we  proceeded  to  fit  the  (VGBpkl,VPN)  datapoints  to  the 
analytical  formulas  developed  in  chapter  3: 

VGBpk    =    VFB    +   VSpk    +    2.sign(VSpk).[Vxx(|VSpk|-kT/q)]1/2  (3.8) 
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Fig.  3.11  Least-squares-fits  of  the  TE-DCIV  curve  of  L=10jim  and  0.25nm  pMOST  at 
VPN=400mV  shown  in  Fig.  3.10.  Decomposition  into  two  components:  one  from 
NIT  in  the  MCR  and  one  from  NIT  in  the  JSCR. 
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Fig.  3.12  Experimental  channel  length  dependence  of  VGB  k  at  IB  k  at  six  L(fxm)  from  0.25 
to  10  for  extrapolation  of  dopant  impurity  concentration  at  the  oxide/silicon 
interface,  (a)  VGB  k  vs  VPN  with  3-parameter  and  2-parameter  least-squares-fit  for 
L=2.0jim.  (b)  VGBpk  vs  L  for  VPN=50mV  to  700mV. 
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where 

VSpk    =    VF    ±    VPN/2    -    loge(cns/cps)1/2  (3.9) 

where  Vs  k  is  the  Si  substrate  surface  energy  band  bending  at  the  peak,  Vxx  = 
£sqNxx/2C0X  where  Cox  =  Eox^oX'  anc*  ^f  =  (kT/q)»loge(Nxx/nj)  is  the  majority 
carrier  quasi-Fermi  potential.  The  term  (-kT/q)  in  (8)  was  dropped  in  (2.11)  and  is 
retained  here  to  provide  better  accuracy  at  high  VPN's  where  Vs  k  approaches  zero  or 
the  flatband  condition  as  indicated  by  the  decreased  drain/base  p/n  junction  barrier 
height  at  the  interface  shown  in  Fig.  3.1.  Combining  (3.8)  and  (3.9)  gives  the  3- 
parameter  formula  to  fit  the  experimental  VGB  k-VpN  data: 

^GBpk 


VGBPk    =    A1+A2+(VPN/2)-2.{A3.  [-A2-(VPN/2)-(kT/q)  ]}1/2 


(pMOST)     (3.10) 
The  three  LSF  parameters  are:  Aj  =  VFB,  A2  =  VF'=  VF  -  (kT/q)»loge(cns/c  s)1/2  and  A3 

=  vxx. 

As  an  example  of  the  LSF  procedure  we  discuss  the  steps  to  fit  the  L=2\im  data 
to  (10)  which  gives  the  solid  line  in  Fig.  3.12(a).  Only  9  (VPN=150mV-550mV)  of  the 
14  datapoints  (50-700mV)  were  used  which  are  the  unfilled  triangles  of  the  L=2p.m  data 
in  Fig.  3.12(a),  in  order  to  limit  the  experimental  VGB  k  uncertainty  to  <0.5mV.  For 
VPN<100mV,  VGB  k  errors  are  larger  than  0.5mV  due  to  the  ~3fA  noise  in  IB.  In  the 
VPN>600mV  range,  a  more  tedious  fit  to  parametric  equations  instead  of  (3.10)  is 
needed  due  to  the  high  injection  level.  The  9-point  VGB  k  data  fits  the  3-parameter 
equation  (3.10)  very  well,  as  shown  visually  in  Fig.  3.12(a),  with  a  Chi-Square  of  only 
0.16%.   The  LSF  value  of  VFB,  VF'  and  Vxx  are  also  listed  in  Fig.  3.12(a).   The 
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VFB=0.96V  is  consistent  with  the  designed  polysilicon  p+Gate  for  pMOST.  The 
electron/hole  capture  rate  ratio  cns/c  is  on  the  order  of  unity  for  near  mid-gap  traps, 
therefore  A2  =  VF'  =  VF,  which  gives  a  surface  concentration  for  the  n-Basewell  donor 
(implanted  As  or  P)  impurity  of  Nxx=NDD=l.lxl018cm'3  which  is  consistent  with  the 
nominal  design  value  of  ~1018cm"3.  The  gate  oxide  thickness  can  then  be  determined 
from  A3=VXX  which  gives  Xox=2.4nm  compared  with  the  nominal  design  value  of 
2.5nm.  The  exponential  dependence  of  Nxx  on  the  Fermi  potential  VF  gives  a  much 
larger  uncertainty  in  Xox  and  Nxx  from  the  6%  uncertainty  in  VF.  The  error  propagates 
as  follows: 

ANdd/Ndd   =   AVF/(kT/q)  (3.11) 

^ox/Xqx    =     <AVXX/VXX    +    ANDD/NDD)/2  (3.12) 

The  data  was  taken  at  room  temperature  giving  kT/q~25mV  and  using  AVF=0.03V,  the 
percentage  error  in  Nxx  is  1 16%  and  Xox,  64%  in  this  global  3-parameter  fit.  However, 
if  one  of  the  three  parameters  (VFB,  Nxx  and  Xox)  is  known  from  another  independent 
measurement  or  is  fixed  in  a  regional  2-parameter  fit,  then  the  two  LSF  parameters  may 
be  determined  more  precisely.  For  example,  we  fixed  the  value  of  VFB  at  0.964V  which 
was  obtained  by  the  3-parameter  fit,  and  then  we  made  a  2-parameter  LSF  to  (3.10). 
The  results  are  listed  in  Fig.  3.12(a)  showing  a  large  reduction  of  uncertainties  in  VF' 
and  Vxx.  The  percentage  errors  translated  to  Nxx  and  Xox  are  only  4%  and  2% 
respectively. 

To  test  the  sensitivity  of  this  2-parameter  extraction  to  the  fixed  value  of  VFB, 
we  vary  VFB  by  ±2%  (~±20mV).   An  alternative  test  would  be  to  fix  Xox  from  an 
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independent  measurement,  such  as  the  C-V.  The  value  of  the  device  parameters  (Nxx 

and  Xox)  ar>d  their  uncertainties,  computed  from  the  A's  from  the  2-parameter  LSF  with 

VFB  fixed,  are  summarized  in  the  following  table: 

vfb<v>  Nxx(1018cirT3)  Xox(nm)  *2(%) 

0.945  0.76+0.02  2.87±0.06  0.15 

0.964  1.13±0.04  2.40±0.05  0.15 

0.984  1.68±0.06  2.00+0.05  0.15 

The  above  table  shows  that  the  DCIV  method  can  extract  Nxx  and  Xox  to  an  accuracy 
of  <10%  provided  that  the  flatband  voltage  is  known  to  better  than  about  (kT/q)/5  ~ 
5mV.  In  the  preceding  analysis,  assumption  is  made  that  cns=cps  so  VF'=VF.  In  factory 
application,  a  calibration  transistor  with  known  substrate  doping  concentration  and 
oxide  thickness  may  be  used  to  determine  the  cns/c  ratio,  from  which  VF  can  be 
extracted  and  the  gate  dopant  impurity  concentration  can  be  computed  from 
^fg=^fb+^f  wnere  Vfg  denotes  the  Fermi -potential  in  the  p+Gate.  For  |VFG|  >  0.50V, 
accurate  gate  doping  concentration  can  be  obtained  only  by  including  impurity 
deionization  and  Fermi-Dirac  statistics.  A  similar  fit  to  the  IBpk2  can  be  made  to  give 
the  Pxx  and  Xox  in  the  LDERs. 

3.4  Profiling  Interface  Traps  in  Junction  Space-Charge  Regions 

Examples  are  now  given  on  profiling  the  interface  traps  in  the  drain  and  source 
junction  space-charge  regions  before  and  after  channel-hot-carrier  (CHC)  stress.  An 
earlier  profiling  technique  employed  the  transient  charge-pump  current  using  its 
dependence  on  the  drain  reverse  biased  [47]  which  was  shown  to  correlate  with  the 
DCIV  method  [48].   A  recent  demonstration  of  the  DCIV  technique  for  profiling 
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interface  traps  near  the  drain  junction  measured  the  differential  basewell-terminal 
conductance  by  varying  the  drain  reverse  bias  [43].  The  differential  DCIV  method  had  a 
higher  sensitivity  and  resolution  and  could  profile  a  large  portion  of  the  channel. 
However,  these  methods  are  all  limited  to  the  region  outside  of  the  reverse-biased 
JSCRs.  The  new  approach,  as  described  in  section  3.2,  gives  the  spatial  NIT  variation 
inside  the  forward-biased  drain  or  source  JSCR  directly  from  the  ratio  of  the 
experimental  DCIV  curve  and  the  reference  theoretical  DCIV  curve  corresponding  to  a 
constant  NIT  distribution.  The  rather  good  spatial  resolution  and  accuracy  of  this  new 
approach  were  already  illustrated  by  using  analytical  NIT(y)  models  in  Figs.  3.4-3.5  and 
described  in  section  3.2.  Oxide  charge  built-up,  which  is  usually  associated  with 
interface  trap  generation  during  CHC  stress  under  large  drain  reverse-bias  [32,49], 
changes  the  local  flatband  voltage  which  complicates  data  analysis.  In  this 
demonstration,  the  CHE  stress  is  kept  at  low  drain  voltages  to  prevent  injection  of 
secondary  hot  holes  over  the  Si02/Si  barrier  (4.25eV  for  hot  holes  and  3.12eV  for  hot 
electrons)  [32]. 

The  example  to  be  given  is  the  extraction  of  NIT(y)  in  a  0.35u>m-nMOST  testbed 
for  reliability  assessment  of  low-K  interconnect  materials  [6].  The  specific  nMOST  in 
this  demonstration  has  a  drawn  aspect  ratio  of  W/L  =  50|im/0.4(j.m.  The  pre-stress  and 
post-stress  DCIV  data  are  plotted  in  Fig.  3.13(a)  with  the  VGB-independent  Ie-baseiine 
subtracted.  The  sharp  pre-stress  peak,  IB  k=0.6pA,  at  VGB  k=-0.5V  is  in  the 
subthreshold  range  and  comes  from  recombination  at  the  interface  traps  in  the  MCR. 
This  peak  gate  voltage  is  consistent  with  the  nominal  XOx=7-0nm  and 
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Fig.  3.13  Demonstration  for  the  extrapolation  algorithm  applied  to  channel-hot- 
electron  stress-generated  interface-trap  concentration  profile  of  an  nMOST 
from  a  0.35fj.m  low-k-interconnect  CMOS  technology,  (a)  Prestress  TE- 
DCIV  and  poststress  DE-DCIV  curves,  (b)  Prestress  and  poststress  NIT(y) 
profile  extrapolations. 


69 

Pxx(boron)=1017cm_3.  The  broad  pre-stress  structure  (~0.45pA)  shown  in  Fig.  3.13(a) 
in  the  range  VGB<-2.5V  (flatband-accumulation  range  for  this  nMOST)  is  from 
recombination  at  NIT  in  the  JSCR  close  to  the  metallurgical  n+Drain/p-Basewell 
boundary  y=0  at  the  Si02/Si  interface  x=0.  This  nMOST  was  then  CHE  stressed  for 
3000  seconds  at  VGB=+1.5V,  VDB=+5.0V  and  VSB=0V.  No  positive  oxide  charging 
was  observed  under  this  stress  condition  since  the  kinetic  energy  of  the  hot  carriers  is 
only  qVGB-q(VDB-VTH)  =  5.0eV-(1.5eV-0.5eV)  =  4.0eV,  which  is  insufficient  for  the 
secondary  hot  holes  to  surmount  the  4.25eV  hole  barrier  at  Si/Si02  interface  [49]. 
Figure  3.13(a)  shows  that  after  the  CHE  stress,  the  broad  pre-stress  structure 
(VGB<-2.5V)  grew  almost  200  times  in  amplitude  (to  85pA)  and  became  even  broader 
(VGB<-1  V).  This  indicates  a  substantial  increase  of  interface  trap  density  in  the  JSCR. 
A  small  post-stress  hump  can  be  observed  around  VGB=-0.5V  which  was  the  peak  gate 
voltage  due  to  the  pre-stress  NIT  in  the  MCR,  indicating  that  some  interface  traps  were 
also  generated  in  the  MCR  by  the  secondary  hot  holes  generated  by  the  hot  electrons  in 
the  JSCR,  but  at  much  smaller  magnitude  than  that  in  the  JSCR  due  to  attenuation  of  the 
hot  hole  density  [32]. 

The  pre-  and  post-stress  IB-VGB  data  in  Fig. 3. 13(a)  can  be  analyzed  using  (3.5) 
to  give  the  lateral  distribution  of  the  interface  trap  density.  The  results  are  shown  in  Fig. 
3.13(b).  The  pre-stress  NIT  was  about  2xl010cm~2  from  the  p/n  boundary,  y=0,  to 
y=10nm  into  the  JSCR  on  the  basewell  or  channel  side,  and  it  dropped  to  a  nearly 
constant  value  of  5xl09cm-2  beyond  y=16nm.  The  post-stress  NIT  increased  to 
8xl012cm-2  within  18nm  of  the  p/n  boundary  y=0  and  dropped  continuously  towards 
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the  MCR,  and  is  below  1012cm~2  at  y=40nm.  The  sharp  rise  near  y=0  or  the  D-LDER  in 
0<y<3nm  may  come  from  the  error  due  to  the  parabolic  potential  of  the  depletion 
approximation  and  the  resolution  limitation  of  steeply  varying  NIT(y),  as  indicated  by 
the  theoretical  models  shown  in  Figs.  3.4  and  3.5.  Thus,  this  example  shows  a  3nm 
spatial  resolution. 

3.5  Interface  Trap  Profile  in  the  Channel  Region 

A  single  DCIV  curve  would  not  provide  sufficient  information  to  profile  the 
interface  traps  in  the  MCR.  Because  the  surface  potential  is  essentially  constant  along 
the  entire  length  of  the  MCR,  every  interface  trap  along  the  MCR  gives  the  same  IB-VGB 
dependency  in  the  TE-DCIV  measurement  regardless  of  the  NIT  location.  This  is 
readily  seen  by  the  general  proof  using  (3.1)  for  the  interface  traps  in  the  MCR  which 
reduces  to  the  simple  formula: 


IB    "     IB1 


LEL 


NIT(y)Wdy  h    ib1    *    NITA(LEL)  (3.13) 

0 


IB1=qRssl  is  the  basewell  recombination  current  from  a  single  interface  trap  and  NITA  is 
the  number  of  interface  traps  in  the  MCR  with  the  electrical  channel  length  of  LEL 
between  the  edge  of  source  and  drain  JSCRs.  Thus,  the  spatial  variation  of  NIT(y)  is 
given  by 

NIT(LEL)*W=     (d/dLEL)NITA(LEL)     =     (d/dLEL)  (IB/IB1)  (3.14) 

Since  LEL  is  modulated  by  VPN,  the  forward  bias  applied  to  the  Source/Basewell  and 
Drain/Basewell  junctions,  the  VpN-dependent  IB  k  measured  in  MOSTs  covering  a  range 
of  drawn  channel  length,  L,  can  give  the  NIT(y)  profile  near  the  basewell  edge  of  the  two 
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JSCR's.  In  addition  to  TE-DCIV  measurement,  the  SE-DCIV  and  DE-DCIV 
measurements  would  provide  further  information  regarding  to  the  symmetry  of  the  NIT 
distribution  by  utilizing  the  additional  variable,  i.e.,  the  nearly  linear  variation  of  the 
injected  minority  carrier  interfacial  concentration  from  the  emitter,  such  as  the  DE,  to 
the  collector,  such  as  the  zero-biased  Source  Collector  or  SC  in  the  DE-DCIV  bias 
configuration.  In  DE-DCIV,  interface  traps  near  the  drain-emitter  (DE)  dominate  the  IB 
while  in  SE-DCIV,  those  near  the  SE  dominate.  Thus,  large  asymmetry  or  differences 
between  DE-DCIV  and  SE-DCIV  would  be  expected  in  the  stressed  transistors  and  in 
pre-stressed  transistors  with  designed  drain-source  asymmetries  such  as  the  FLASH 
memory  transistor. 

The  above  differential  method  will  be  applied  to  the  0.25u>m  family  of  pMOSTs 
with  6  drawn  channel  lengths  whose  typical  TE-DCIV  data  at  15-VPN's  (OmV  to  700mV 
with  50mV  steps)  were  shown  in  Fig.  3.10(a)  for  L=10|xm  and  Fig.  10(b)  for  L=0.25^im. 
The  I^'s  were  obtained  from  the  LSF  to  the  Gaussian  as  described  by  eq.  (3.7)  with 
uncertainties  of  4-8%.  The  ratio  IB  k(L)/IB  k(L=10jum)  versus  VPN  for  the  five  drawn 
L's  are  plotted  in  Fig.  3.14(a)  and  the  4%-8%  uncertainties  of  the  datapoints  are  not 
perceptible.  The  nearly  VPN-independent  value  of  this  ratio  is  plotted  as  a  function  of  L 
in  Fig.  3.14(b)  for  three  VPN's  (300,  400  and  500mV)  showing  two  features:  (i)  a  rise 
proportional  to  L,  suggesting  an  increasing  dominance  of  the  NIT  in  MCR  due  to 
increasing  length  or  recombination  area,  and  (ii)  an  asymptotic  constant  value  of  0.25 
with  decreasing  L,  suggesting  a  large  contribution  to  IB  from  NIT  inside  or  near  the  edge 
of  the  drain  and  source  JSCRs  due  to  its  L-independence.  The  data  points  can  be  fitted 
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Fig.  3.14  NIT  profile  extrapolation  from  channel  length  dependence  of  DCIV  peak 
current,  (a)  Normalized  IBpk  vs  VPN  for  L(|j.m)=0.25  to  10.  (b)  Normalized 
IBpk  vs  L  at  VPN=300,  400,  and  500mV.  Two  components  are  evident,  one 
each  from  the  MCR  and  JSCR. 
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to  the  following  geometric  equation  where  L  is  Lithographic  or  drawn  channel  length 
which  will  be  taken  as  the  distance  between  the  metallurgical  drain  and  source  p/n 
junction,  LMCR  is  the  length  of  the  mid-channel  region,  YD(=YDJSCR)  and  YS(=YSJSCR) 
are  respectively  the  length  of  the  drain  and  source  junction-space-charge  region  at  the 
Si02/Si  interface  which  are  equal,  YD+Ys=2YD=Yj  for  conventional  symmetrical 
cMOS  transistor  design,  and  YMCR  =  L  -  Yj,  and  where  the  J's  are  the  areal  density  of 
the  recombination  current  in  the  MCR  and  the  drain  and  source  JSCR  as  indicated  by 
the  subscripts  which  are  to  be  shortened  to  two  characters,  MC,  DJ,  SJ. 

^pk^)  JMCr(L  Yj)      +     JDJSCR*YD     +     JSJSCR*YS 


IBpk(1°)  JMCR<10"     Yj)      +     JDJSCR*YD     +     JSJSCR*YS 

=  y  =     (x    +    A^/dO    +   Ax)  (3.15) 

where  x=L,  Al  =  (JDJ*YD  +  Jsj*Ys)/Jmc  -  Yj  which  reduces  to  Al  =  [(JDJ/JMc)  ~  !]*Yj 
for  identical  drain  and  source.  The  data  in  Fig.  3.14(b)  gives  Aj  =  3.25±0.13|im.  For  an 
abrupt  source  and  drain  n+/p-B  junction  of  ~1018cirT3,  YD=Ys=Yj/2=0.04|im  giving 
(JJSCR/JMCR)=NIT(SCR)/NIT(MCR)  ~  40.  This  suggests  that  the  NIT(y)  profile  rises 
sharply  in  a  short-distance  (<0.04u,m  or  <40nm)  in  the  MCR  near  or  inside  the  JSCR, 
which  was  depicted  by  the  theoretical  example  of  a  sharper  exponential  profile, 
NIT(y)/NIT(0)  =  exp(-y/3nm)  given  in  Fig.  3.4(a). 

The  geometric  modeling  from  the  L-dependence  just  presented  that  showed  a 
MCR  trap  and  a  high-NIT  JSCR  trap  can  be  analyzed  more  thoroughly  using  the  IB 
integral  given  by  (3.1)  to  cover  the  entire  channel.  In  view  of  the  good  accuracy  and 
resolution  demonstrated  by  the  analytical  model  of  an  exponential  NIT  profile  in  the 
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JSCR  shown  in  Fig.4(a)-(c),  we  shall  assume  the  exponential  NIT  profile  again  to  model 

the  L-dependence  that  showed  a  high-NIT  near  the  edge  or  inside  of  the  JSCR  over  the 

entire  VPN  range  (50mV  to  700m V)  which  sweeps  this  edge  over  a  region  about  20nm 

to  40nm  on  the  MCR  side  of  the  p+D/n-B  boundary.   Thus,  the  NIT  profile  can  be 

modeled  by 

NIT(y)     =    A1    +   A2{exp(-y/A3)     +    exp  [  (y-YMCR) /A3  ]  }  (3.16) 

Aj  is  the  constant  concentration  in  the  middle  of  the  channel  in  asymptotic  long  channel, 

A2  is  the  peak  concentration  near  the  edges  of  the  two  JSCR's  and  A3  =  LIT  is  the  decay 

length  of  the  NIT  distribution. 

With  the  TE-DCIV  (both  Drain  and  Source  are  forward-biased  at  the  same  VPN) 

the  minority  carrier  interfacial  concentration,  P(x=0,y),  in  the  n-Basewell  of  a  pMOST  is 

spatially  constant  in  the  MCR  from  y=0  to  y=YMCR  (y-axis  shifted).  In  the  DE-DCIV  or 

SE-DCIV  bias  configuration,  it  is  linearly  decreasing  from  the  forward-biased  drain  (or 

source)  emitter  to  the  zero-biased  source  (or  drain)  collector.  These  are  given  below 

with  the  origin  of  the  y-axis  shifted  to  the  collector-edge  of  the  MCR. 

P(x=0,y)    =   Ps(y)    =   PS(YMCR)  TE-DCIV 

(3.17) 
=    ps(YMCR)*(y/YMCR)       DE-DCIV,     SE-DCIV 

These  three  minority  carrier  concentration  profiles  can  be  used  with  the  symmetrical 
exponential  NIT  profile  of  (3.16)  to  evaluate  the  IB  integral  such  as  (3.1)  with  the 
integration  limit  extended  over  the  entire  channel.  To  provide  the  analytical  solutions 
for  IB  k,  we  shall  use  the  peak  value  of  RSS1  in  (3.1),  which  occurs  at  c  Ps  =  cnsNs.  For 
En  ~  0  or  cpsPs=cnsNs  »  eps  and  ens,  it  is  given  by 
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Rssi-pk   =     d/2)  (cpsPscnsNs)1/2  (3.18) 

=     (1/2)  (cpscns)1/2[niexp(qVPN/2)  ]  TE-DCIV 

=     (1/2)  (cpscns)1/2[niexp(qVPN/2)  ]  (y/YMCR)1/2      DE/SE-DCIV 

Then,  from  (3.1)  using  (3.16)  and  (3.18),  and  L  instead  YMCR  for  the  integration  limit, 

and  noting  that  YMCR  =  L  -  2YJSCR(VPN),  we  have 

IBPk    =     f  A1L+2A2A3.  [l-exp(-L/A3)  ]  }IB1W       TE-DCIV  (3.19) 

=     {  (2/3)A1L+aA2A3.  [l-exp(-L/A3)  ]  }IB1W       DE/SE-DCIV    (3.20) 

Where  a- 1  if  A3  «  L,  or  the  exponential  NIT(y)  is  localized  near  the  Source  and  Drain. 

In  the  other  limit  that  NIT(y)  varies  slowly  across  the  channel,  A3  >  L,  cc=4/3.  These 

show  that  IBpk  of  TE-DCIV  is  1.5  or  2.0  times  larger  than  that  of  DE-DCIV  and  SE- 

DCIV  for  very  long  and  short  channels  respectively. 

The  TE-DCIV,  DE-DCIV  and  SE-DCIV  data  at  VPN=400mV  are  fitted  to  (3.19) 

giving  the  following  results  which  are  used  for  the  LSF  lines  in  Fig.  3.15. 

A1IB1W(pA/Um)  A2IB1W(pA/um)  A3  (run) 

TE  4.3±0.2  90±40  80±40 

DE  4.1±0.2  90±50  80+50 

SE         4.8±0.2 

These  show  consistency  among  the  three  bias  configurations  and  a  NIT(y)  that  rises  from 
the  MCR  by  a  factor  of  20  to  the  edge  MCR  or  JSCR  (90/4.2=21).  The  two  LSF 
parameters,  A!  and  A2,  can  be  converted  to  NIT.  Using  W=10u,m,  nj=1010cm-3,  and 
assuming  cns=c  =10_8cm3/sec,  we  obtained  Aj=2.2xl09cm-2  and  A2=4.68xl010cm"2 
for  the  interface  trap  densities  in  the  middle  of  the  channel  and  near  the  drain  and  source 
junctions. 


76 


10-9E 


< 

Q_10-io 


XL 

Q. 


Q010        = 


10 


-12 


=        I      I    I  I  I  INI          I      I    I  I  I  III 

I     DCIV  Data    LSF  Eq(17) 
TE       o        — 

DE         a 

SE      o 

I           I        I      I     I    III- 

pMOST  = 

W=10|im     : 
VPN=400mV  - 

ii |       i    i  i  m 

_3 * 

— 

i     i    i 

I     I    I    III                    I           I        I      I     I    I    III 

I 

I     MINI 

11  — 


10 


-1 


1  10 

LEL/(1|am) 


10s 


Fig.  3.15  Forward  bias  configuration  dependence  of  the  Iepk'L  data  from  TE-DCIV, 
DE-DCIV  and  SE-DCIV  measurements  at  VPN=400mV,  showing  the  0.5 
ratio  at  short  L  from  contributions  of  NIT  located  inside  and  near  the  edge  of 
the  ~30nm  JSCR  and  the  0.66  ratio  at  long  L  from  NIT  located  in  the  MCR 
which  are  illustrated  more  explicitly  in  Fig.  3.16. 
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The  experimental  data  for  the  DE/TE  and  SE/TE  ratios  of  the  IB  k  are  shown  in 
Fig.  3.16(a)  and  (b)  as  a  function  of  forward  bias,  VPN,  showing  the  predicted  variation 
with  channel  L  by  (3.19)  from  0.5  to  0.667.  Finally,  the  LSF  profile  are  plotted  in  Fig. 
3.17  at  VPN=400mV  with  an  estimated  YJSCR=40nm  for  the  shorter  channels, 
L=0.25u,m,  0.4u,m,  0.8|im  and  2.0fim. 

Outside  of  the  MCR  or  in  the  JSCRs,  NIT  seemed  to  drop  as  shown  by  the  two 
broken  lines  in  Fig.  3.17,  as  suggested  by  three  experimental  evidences:  (1)  The  VPN- 
independence  of  the  IB  k  in  Figs.  3.14(a),  3.16(a)  and  3.16(b);  and  channel-length 
dependence  of  the  (2)  lineshape,  and  (3)  VGB  k  as  indicated  in  Fig.  3.11.  The  0.25u.m 
pMOST  has  a  higher  shoulder  and  a  larger  HWHM  at  the  base-accumulation  side  of  the 
peak  which  is  consistent  with  contributions  from  NIT  in  the  JSCR  anticipated  by  the 
analytical  theory  discussed  in  section  3.3.  There  is  also  a  20mV  shift  of  peak  gate 
voltages  between  the  L=10u,m  and  0.25u,m  curves  which  could  be  from  two 
possibilities:  (1)  NIT  in  the  JSCR  near  MCR  edge  would  shift  VGBpk(0.25u.m)  pMOST 
towards  flatband-accumulation  (positive  VGB)  with  respect  to  VGB  k(10fxm)  and  (2) 
injection-diffusion-limited  minority  carrier  reduction  in  the  10u.m  pMOST  would  shift 
VGB  k(10|im)  towards  subthreshold-inversion  (negative  VGB  direction)  with  respect  to 
VGB  k(0.25u,m).  The  first  explanation  is  favored  because  NIT  inside  the  JSCR  and  near 
the  MCR/JSCR  boundary  can  account  for  the  lineshape  change  as  well  as  peak  gate 
voltage  shift  as  indicated  in  Fig.  3.1 1.  Furthermore,  injection-diffusion-limitation, 
proportional  to  L,  is  expected  to  dominate  at  low  VPN  due  to  its  lower  value  from  the 
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Fig.  3.16  Channel  length  dependence  of  the  Iepk'^PN  anc*  ^ts  theoretical  limits  for 
short  and  long  channels,  (a)  DE-IB  k/TE-IB  k  and  (b)  SE-IB  k/TE-IB  k, 
showing  the  theoretically  anticipated  0.5  ratio  at  short  L  from  NIT  located 
inside  and  near  the  edge  of  ~30nm  thick  JSCR  and  the  0.66  ratio  at  long  L 


from  N,T  located  in  the  MCR. 
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Fig.  3.17      NIT(y)  distribution  for  L(u.m)=0.25  to  2.0  of  pMOSTs  extrapolated  from 
experimental  data  using  the  profiling  algorithm  described  in  the  text. 
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exp(qVPN/kT),  which  is  inconsistent  with  the  nearly  L-independence  of  the  experimental 

^GBpk  at  a^  ^pn's  shown  in  Fig-  3.12(b). 

The  presence  of  the  interface  traps  inside  JSCR  should  not  affect  the  validity  of 
MOST  parameter  extraction  using  the  VPN  dependence  of  VGB  k  as  described  in  section 
3.4.  The  parallel  shift  in  Fig.  3.12(b)  indicates  that  the  accuracy  of  oxide  thickness  and 
substrate  doping  concentration  extraction  would  not  be  affected  by  the  channel  length  of 
the  MOST.  However,  there  could  be  a  small  offset  (<20mV)  in  the  extracted  flatband 
voltage. 

3.6  Summary 

A  self-consistent  picture  on  the  effects  of  spatial  distributions  of  interface  traps 
on  DCIV  characteristics  has  been  obtained  and  presented.  Three  applications  of  the 
DCIV  method  have  been  demonstrated  with  experimental  data:  (1)  Extraction  of 
transistor  design  parameters  using  the  weak  inversion  gate  voltage  at  the  peak  base 
terminal  current  from  carrier  recombination  at  interface  traps  in  the  channel  region;  (2) 
Profiling  the  concentration  of  process-residual  and  channel  hot  carrier  stress  generated 
interface  traps  in  the  forward-biased  drain  junction  space  charge  region  using  a  portion 
of  DCIV  curve  between  the  weak  inversion  gate  voltage  and  strong  accumulation  gate 
voltage;  and  (3)  Profiling  the  concentration  of  process  residual  interface  traps  in  the 
channel  region  using  the  channel  length  dependence  of  base-current  peak.  These 
demonstrations  provide  further  application  capabilities  for  industrial  applications  of  the 
DCIV  method  [50,51]. 


CHAPTER  4 

GATE  TUNNELING  CURRENTS  IN 

ULTRATfflN  OXIDE  MOS  TRANSISTORS 

4.1  Historical  Survey 

Electron  tunneling  through  a  classically-forbidden  energy  barrier  was  first 
employed  by  Fowler  and  Nordheim  [52]  (F-N)  in  1928  to  successfully  explain  the 
emission  of  electrons  from  the  surface  of  a  tungsten  electrode  into  vacuum  under  the 
influence  of  an  external  electric  field  of  some  106  V/cm.  Schottky  attempted  to  explain 
the  phenomena  by  thermionic  emission  through  an  image  force  lowered  surface 
potential  barrier  [53],  but  a  complete  reduction  of  the  4.5eV  barrier  in  tungsten  requires 
a  field  of  at  least  108V/cm.  In  the  F-N  theory  of  field  emission  [54],  the  electrons  inside 
the  tungsten  with  kinetic  energy  lower  than  4.5eV  have  an  energy  dependent  probability 
of  tunneling  through  the  field-induced  triangular  potential  barrier  as  determined  by 
Schrodinger's  wave  mechanics,  while  the  energy  distribution  of  electrons  was  obtained 
following  Sommerfeld's  electron  theory  of  metal. 

With  the  advancement  of  quantum  theory  and  materials  technology, 
investigation  of  tunneling  phenomena  was  extended  to  (1)  extrinsic  semiconductors  first 
by  Zener  [55]  which  eventually  led  to  the  invention  of  p-n  tunnel  diode  by  Esaki  [56]  in 
1957,  and  (2)  metal-oxide-metal  (MOM)  and  metal-oxide-semiconductor  (MOS) 
structures  which  were  made  possible  by  the  thermal  oxidation  technology  to  grow 
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surface  oxide.  In  terms  of  band  theories  of  semiconductors,  Kane  [57]  interpreted  the 
interband  tunneling  in  p-n  diode  as  'electron  penetrating  the  forbidden  gap  along  the 
imaginary  k  axis'  and  under  the  constant-field  approximation  derived  an  electron 
tunneling  theory  that  took  into  account  of  perpendicular  momentum  conservation.  The 
p-n  tunnel  diode  showed  negative  resistance  under  small  forward  biases  due  to  the 
narrowing  of  energy  range  for  direct  interband  tunneling  with  increasing  forward-bias 
but  its  current  never  went  to  zero  even  when  the  constant  energy  tunneling  path  ceased 
to  exist.  This  'excess'  current  was  from  many  trap-assisted  inelastic  tunneling 
pathways,  such  as  those  studied  by  Sah  in  gold-doped  narrow  p-n  diode  [58].  Bardeen 
considered  tunneling  in  MOM  structure  as  transferring  of  a  quasi-particle  across  the 
oxide  in  a  many-particle  system  [59].  He  calculated  the  transition  probability  using  the 
time-dependent  perturbation  approach  and  showed  that  the  matrix  element  for  the 
transition  is  that  of  the  current  density  operator  in  the  barrier  region.  This  was  followed 
by  Harrison  [60]  to  treat  tunneling  through  a  forbidden  region  with  spatially  varying 
band  structures  from  an  independent-particle  point  of  view.  Harrison  showed  that  the 
traditional  Wentzel-Kramers-Brillouin  (WKB)  approximation  can  be  readily  applied  to 
band  structures  which  vary  slowly  over  a  wavelength  of  the  particle,  but  at  the  sharp 
boundaries  the  wavefunctions  must  be  matched  to  assure  the  conservation  of  electron 
flux.  The  latter  resulted  in  an  energy-dependent  pre-exponential  factor  in  the  expression 
of  the  transition  probability  which  was  used  by  Shewchun  and  Temple  [61]  in  their 
computations  of  tunneling  current  through  semiconductor-insulator-semiconductor 
strictures.  Simmons  [62]  introduced  the  average-barrier  approximation  to  evaluate  the 
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WKB  tunneling  probability  for  barriers  of  arbitrary  shape  and  derived  an  analytical 
tunneling  current  formula  for  MOM  structures  at  zero  temperature.  Simmon's  formula 
was  used  by  Chang,  Stiles  and  Esaki  [63]  to  correlate  their  data  on  electron  tunneling 
through  a  trapezoidal  barrier  in  a  Al/Al203/SnTe  MOS  structure  while  neglecting  band 
bendings  in  the  highly-degenerate  p-type  SnTe.  By  making  the  WKB  approximation, 
detailed  tunneling  current  structures  due  to  the  sharp  potential  boundaries  are  obscured. 
Gundlach  [64]  showed  by  solving  the  Schrodinger  equation  exactly  that  tunneling 
current  through  a  triangular  barrier  can  be  an  oscillatory  function  of  applied  voltage  due 
to  partial  reflections  and  interferences  of  electron  waves  in  the  conduction  band  of  the 
insulator  at  the  sharp  potential  boundary. 

Since  the  advent  of  silicon  transistor  technologies  in  the  late  1950's  [1], 
tunneling  through  metal/Si02/Si  structures  has  received  much  attention.  Lenzlinger  and 
Snow  [65]  performed  the  original  experiment  of  electron  tunneling  from  silicon  into 
thermally  grown  Si02  and  found  that  the  tunneling  current  density  as  a  function  of 
electric  field  J(E)  followed  the  classical  F-N  theory.  They  derived  a  constant  effective 
mass  of  0.42m0  in  the  forbidden  energy  gap  of  Si02  and  showed  that  the  simple 
parabolic  dispersion  relationship  is  equivalent  to  that  of  the  empirical  Franz-type  which 
reduces  to  k=0  at  both  conduction  band  edge  E=EC  and  valence  band  edge  E=EV.  We 
note  that  this  equivalence  is  only  valid  in  the  triangular  barrier  range  (or  F-N  range) 
while  some  recent  authors  applied  it  to  trapezoidal  barriers.  The  band  bandings  and 
carrier  quantization  effects  in  silicon  inversion  layer  were  considered  by  Weinberg 
[66,67]  and  those  of  silicon  accumulation  layer  considered  by  Krieger  and  Swanson  [68] 
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who  also  included  transverse  momentum  conservation  and  image  force  correction  in 
their  theory.  In  view  of  the  accuracy  of  Weinberg's  tunneling  J(E)  data  [66]  on  thick 
-lOOnrn  Si02  films,  we  will  use  them  as  a  reference  to  calibrate  our  analytical  tunneling 
theory  in  the  F-N  range.  On  thin  ~5.0nm  oxide  MOS  structures,  Lewicki  and  Maserjian 
[69]  observed  oscillatory  tunneling  J(V)  characteristics  in  the  F-N  range  as  predicted  by 
the  Gundlach  theory.  From  the  decay  of  the  oscillation  amplitude  with  the  tunneling 
path  length  in  Si02  conduction  band,  they  obtained  a  value  of  scattering  mean-free-path 
(MFP)  in  Si02  conduction  band,  which  was  later  re-interpreted  to  be  0.65nm  by 
Fischetti  et  al.  [70]  who  correlated  this  MFP  to  the  electron-LO  phonon  scattering. 
Recently,  the  oxide  fields  at  which  oscillatory  extrema  occur  were  used  by  Zafar  et  al. 
[71]  to  extract  oxide  thickness  in  the  range  of  4.0nm  to  6.0nm. 

In  MOS  structures,  tunneling  carrier  species  includes  electrons  and  holes,  which 
were  first  separated  by  Weinberg  et  al.  who  diffused  a  shallow  p-n  junction  in  the  silicon 
substrate  under  the  insulator  [72,73].  In  a  MOS  transistor  or  a  MOS  capacitor 
surrounded  by  a  source  region  (sourced  MOS  capacitor,  or  SMOSC),  the  inversion 
channel  and  the  well/substrate  p-n  junction  (if  available)  provide  a  means  for  charge- 
carrier  separation.  However,  this  separation  is  not  unambiguous  because  the  carrier  type 
detected  are  those  in  the  terminals  of  the  p-n  junction  or  the  MOST  rather  than  those  in 
the  oxide.  The  following  three  processes  lead  to  the  detection  of  hole  currents  as  a 
result  of  electron  tunneling:  (1)  valence  electron  tunneling  from  Si  leaves  a  hole  behind 
(2)  electron  tunneling  from  the  gate  electrode  can  impact  generate  a  hole  or  (3) 
recombine  with  a  hole  in  the  silicon  substrate.  Careful  analysis  of  charge-separation 
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data  were  demonstrated  by  Yau  [74]  in  his  identification  of  electron  as  the  major 
conduction  carrier  in  Si3N4.  Eitan  and  Kolodny  [75]  observed  a  hole  current  in  the  p- 
substrate  terminal  of  a  7.8nm  oxide  n+sourced  MOS  capacitor  under  positive  gate  bias 
and  attributed  it  to  valence  electron  tunneling  from  the  p-substrate.  This  was  argued 
against  by  Weinberg  and  Fischetti  [76]  because  their  theoretical  estimates  for  valence 
electron  tunneling  was  about  103  times  lower  than  the  experimental  data.  Recently, 
Rasras  et  al.  [77]  used  a  remote  collector  to  show  that  this  hole  current  under  large 
positive  gate  voltage  was  due  to  electron-hole  pair  generation  in  the  p-substrate  by 
photons  travelled  from  the  gate  electrode.  For  thin  oxide  (<3.0nm)  under  low  (<3V) 
gate  voltages,  Shanware  et  al.  [78]  showed  that  their  valence  band  tunneling  theory  with 
a  tunneling  effective  mass  of  0.33m0  correlated  with  the  hole  substrate  currents  from 
which  oxide  thickness  can  be  extracted. 

At  the  Si/Si02  interface,  the  valence  band  offset  is  4.25eV  while  the  conduction 
band  offset  is  3.13eV  [79].  The  1.12eV  difference  in  tunneling  barrier  height  makes  the 
hole  tunneling  much  smaller  in  magnitude  and  elusive  for  experimental  detection. 
However,  secondary  hot  hole  generation  [80]  and  injection  [81]  from  the  anode  initiated 
by  F-N  tunneling  of  primary  electrons  from  the  cathode  was  identified  to  cause  positive 
oxide  charge  and  dielectric  breakdown.  The  threshold  Si02  voltage  to  give  secondary 
hole  over-the-barrier  injection  was  found  by  Lu  and  Sah  [81]  to  be  (1)  4.25V  by  Auger 
generation  and  (2)  5V  by  impact  generation  of  energetic  holes.  Primary  hole  tunneling 
was  not  reported  until  very  recently  in  p+gate  pMOST  [82-84]  at  the  low  inversion  gate 
voltage  range,  but  no  quantitative  correlation  with  hole  tunneling  theory  was  given. 
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4.2  Recent  Technology  Motivations 

Literatures  on  fundamental  theories  and  experiments  of  tunneling  are  abundant 
as  reviewed  in  the  first  section.  Most  of  them  focused  on  tunneling  through  a  triangular 
barrier.  Over  the  past  three  decades,  semiconductor  industry  has  sustained  a  growth  rate 
of  15%  per  year  which  is  fueled  by  relentless  down  scaling  of  silicon  transistors.  The 
channel  length  of  a  MOST  is  scaled  from  25u,m  in  1962  [1]  to  about  0.13u,m  in  the 
current  production  technology.  Following  Dennard's  constant-field  scaling  law  [85], 
gate  oxide  thickness  and  power  supply  voltage  need  to  be  scaled  accordingly.  Oxide 
thickness  is  scaled  down  to  about  3.0nm  in  production  technologies  while  more 
aggressive  scaling  to  below  1.5nm  has  been  pursued  in  research  and  development 
establishments  [86-87]  which  is  limited  by  oxide  tunneling  currents.  Power  supply 
voltage  is  limited  by  the  MOST  threshold  voltage  to  about  1.0V.  At  such  values  of 
oxide  thickness  and  power  supply  voltage,  the  shape  of  the  potential  barrier  for 
tunneling  is  trapezoidal  rather  than  triangle,  and  gate  oxide  tunneling  could  become  a 
limiting  factor  in  MOST  off-state  power  dissipation.  Also,  tunneling  current  hinders 
transistor  parameter  extraction  by  traditional  measurement  techniques  such  as 
capacitance-voltage  method  [10,13,88].  Accurate  modeling  of  trapezoidal  tunneling 
through  ultrathin  (<3.0nm)  oxide  is  needed  in  new  transistor  technology  development. 

Tunneling  through  a  trapezoidal  Si02  barrier  has  been  modeled  for  conduction 
band  electron  tunneling  from  nMOST  inversion  layers  [89-92]  and  from  pMOST 
accumulation  layers  [93,94].  When  the  silicon  surface  is  strongly  inverted  or  strongly 


87 

accumulated,  electrons  at  the  Si/Si02  interface  are  confined  by  the  surface  electric  field 
on  the  silicon  side  and  the  3.13eV  potential  barrier  on  the  Si02  side.  Due  to  the  finite 
probability  of  tunneling  through  the  Si02  barrier,  electrons  are  in  quasi-bound  states. 
Full  quantum  mechanical  models  compute  the  tunneling  current  by  solving  the  lifetime 
of  the  quasi-bound  state  [89,91,93],  which  is  computationally  demanding  and  not 
suitable  for  compact  modeling.  A  comparative  study  of  quantum  mechanical  and 
classical  modeling,  however,  reveals  similar  results  in  tunneling  current-voltage 
characteristics  [93].  This  maybe  due  to  the  two  quantum  mechanical  effects  that  act 
against  each  other  as  far  as  tunneling  current  is  concerned:  (1)  energy  quantization  gives 
a  ground  state  energy  above  the  silicon  conduction  band  minimum  which  reduces  the 
effective  barrier  height,  and  (2)  the  centroid  of  electron  distribution  in  the  surface  layer 
moves  away  from  the  Si/Si02  interface  which  increases  'tunneling  distance'. 
Analytical  approximations  to  evaluate  lifetimes  [92,94],  however,  render  the  tunneling 
current  as  a  function  of  oxide  electric  field  and  oxide  thickness  only.  Those  results  give 
zero  tunneling  current  at  zero  oxide  electric  field  rather  than  at  the  zero  gate  voltage, 
which  is  inaccurate  for  any  MOST  structures  that  has  work  function  difference  between 
gate  and  substrate  materials. 

We  will  develop  a  compact  tunneling  current  model,  which  is  applicable  to  the 
full  gate  voltage  range  for  both  nMOSTs  and  pMOSTs.  The  model  will  treat  conduction 
band  electrons,  valence  band  electrons  and  valence  band  holes  on  equal  basis.  The  latter 
is  indispensable  to  describe  the  tunneling  currents  in  the  inversion  range  of  a  p+gate 
pMOST.   Carrier  concentrations  in  silicon  layers  will  be  taken  into  account  both 
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explicitly  and  implicitly,  the  latter  through  the  gate-voltage  vs.  oxide  field  relationship, 
and  tunneling  in  both  directions  will  be  included.  This  will  allow  an  accurate 
comparison  of  relative  importance  of  various  energetic  tunneling  pathways  and 
geometric  tunneling  pathways  in  MOSTs  at  low  gate  voltages.  We  will  correlate  the 
model  to  experimental  data  and  explore  its  utility  for  transistor  parameter  extraction. 

4.3  Theory  of  Tunneling  in  MPS  Structures 

In  deriving  our  tunneling  current  formula,  we  follow  Bardeen's  transition 
probability  approach  [59]  and  Harrison's  independent-particle  tunneling  model  [60]. 
The  probability  per  unit  time  of  transition  for  an  electron  in  state  a  on  one  side  of  the 
tunneling  barrier  to  state  (3  on  the  other  side  is  described  by  Fermi's  golden  rule: 

Pa/3    =     (2JT/B)  |Ma/3|2    ppfa{l-fp)  (4.1) 

Map  is  the  matrix  element  to  be  evaluated  between  the  two  states  a  and  (3  in  the 
tunneling  region  where  M  is  proportional  to  the  current  density  operator  as  shown  by 
Bardeen  [59].  pp  is  the  density  of  states  at  (3,  and  fa  and  f«  are  occupation  probabilities 
of  the  two  states  respectively. 

The  net  tunneling  current  density  is  derived  by  summing  over  all  allowed  states 
a  and  subtracting  the  reverse  tunneling  current  under  the  constraint  of  transverse 
momentum  (kt)  conservation  [60]: 


4;rq 

J   =   X 

h        K 


lMa/?l       PaPfS^a-£^     dE  <4-2> 

0 


The  additional  factor  2  accounts  for  spin  degeneracy,  q  is  the  charge  of  an  electron  and 
zero  to  Em  represents  allowed  energy  range  for  tunneling.   The  summation  over 
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transverse  momentum  can  be  converted  to  integrals  if  the  electrons  are  free  to  move  in 
the  transverse  directions:  Z  — >  l/(27t)2jdkyjdkz.  The  double  integration  over  transverse 
momentum  can  be  carried  out  explicitly  for  the  Fermi  functions  (fa-fp).  The  matrix 
element  can  be  evaluated  if  the  wavefunctions  inside  the  tunneling  barrier  for  state  oc 
and  |3  are  written  in  terms  of  WKB  approximations.  The  result  is  the  familiar 
exponential  WKB  tunneling  probability.  Figure  4.1  shows  an  energy  band  diagram  of  a 
Si/Si02/Si  trapezoidal  barrier.  The  Si/Si02  conduction  band  and  valence  band  offsets 
are  denoted  by  <J)Bn  and  <j)B  respectively.  The  initial  state  on  the  right  side  of  the  barrier 
with  kinetic  energy  Ej  and  wavefunction  X¥SI  tunnels  to  the  final  state  on  the  left  side 
with  kinetic  energy  E2  and  wavefunction  Y^.  At  the  two  abrupt  boundaries,  x=0  and 
x=Xox,  both  *F  and  (l/m^dW/dX  must  be  continuous  to  conserve  electron  flux. 
Wavefunction  matching  of  *PS1  at  x=0,  and  W^  at  x=X0x  results  in  an  electron  kinetic 
energy  dependent  pre-exponential  factor. 

Applying  the  above  for  electron  tunneling  through  Si/Si02/Si  structure  (see 
Appendix  for  detailed  analysis),  we  derive  the  following  tunneling  current  formula  for  a 
single  ellipsoid  conduction  band  sub-valley  in  the  silicon: 
4  7TqmdkT 


JN  - 


h3 


xlog. 


T0(Ex)exp[-2u(Ex) ] 
0 


l+exp[ (EF-Ec-Ex)/kT] 


dEx  (4.3) 


l+exp[ (EF-Ec-Ex-qVGB)  /kT] 

md  is  the  density-of-states  effective  mass  of  the  sub-valley,  EX=(E-Ec)-Et  is  the 
electron  kinetic  energy  in  the  tunneling  direction  x,  EF  is  the  Fermi  energy,  Ec  is  the 
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Fig.  4.1  Trapezoidal  barrier  for  electron  tunneling  in  a  Si/Si02/Si  structure. 
Electron  wavefunctions  of  the  initial  and  final  states  are  denoted  by  Wsl  and 
4^2  respectively.  Boundary  conditions  at  x=0  and  x=Xox  must  be  satisfied 
to  conserve  electron  flux  across  the  two  abrupt  boundaries. 
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conduction  band  minimum  in  the  semiconductor  near  the  tunneling  boundary.  T0(EX) 
and  i)(Ex)  are  the  pre-exponential  factor  and  the  exponent  of  the  WKB  tunneling 
probability  respectively  to  be  described  as  follows  in  conjunction  with  figure  4.1. 


16 


T0(EX)     = 


moxklx(°)  ^l^ox^H 

+    


msl^ox(°>  moxklx(°)J      Lms2Kox<XOX>        moxk2x<XOX> 


lox^2xv^OX'         ms2Nox  v^OXJ 

+ 


(4.4) 


and 


V(Ey)    = 


"•ox 


K(x)dx 


.0  (4.5) 

msi  (i=l,  2)  are  the  normal  effective  mass  in  the  two  silicon  regions,  and  mox  is  the 
tunneling  effective  mass  in  the  Si02  layer  with  two  boundaries  at  x=0  and  x=Xox 
respectively.  klx(0)  and  k2x(Xox)  are  rea'  wave  numbers  at  the  silicon  side  of  the 
boundaries,  while  k(0)  and  k(Xox)  are  imaginary  wave  numbers  at  the  Si02  side  of  the 
boundaries.  Assuming  parabolic  conduction  band  of  Si02  and  extending  the  dispersion 
relationship  analytically  into  the  energy  gap,  we  have  the  traditional  expression  for  k, 
which  is: 

K(x)     =     [2mox(0Bn-qEoxx-Ex)]1/2/ri  (4.6) 

Eox  is  the  electric  field  inside  Si02  and  is  constant  for  a  trapezoidal  barrier.  We  will 
include  some  discussions  to  the  above  formulas  (4.3)  to  (4.6)  in  the  following 
paragraphs. 

The  WKB  approximation  is  valid  when  the  de  Broglie  wavelength  (X)  of  the 
electron  (or  its  imaginary  analogue  inside  the  barrier)  is  much  smaller  than  the 
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characteristic  length  over  which  the  potential  varies  appreciably  [95].  For  the  Si/Si02 
electron  barrier,  X  can  be  estimated  as  follows: 

1  h 

1    =    =    =    1.5A 

K  (2mox<V1/2  (4.7) 

This  result  shows  that  the  WKB  approximation  should  be  valid  for  oxide  thickness  as 
thin  as  10-15A  which  is  already  6-10  times  larger  than  the  (imaginary)  electron 
wavelength. 

The  pre-exponential  factor  T0(EX)  of  (4.4)  is  valid  for  trapezoidal  barrier  only  in 
which  both  boundaries  are  abrupt.  This  kinetic-energy  dependent  factor  plays  an 
important  role  in  explaining  the  gate  voltage  dependence  of  thin  oxide  tunnel  data  as 
will  be  shown  later.  It  gives  the  largest  correction  near  the  flatband  voltage  as  compared 
to  the  traditional  unity  pre-exponent.  Equation  (4.4)  is  symmetric  to  the  two  silicon 
layers,  which  is  consistent  with  the  notion  that  the  tunneling  probability  of  a  trapezoidal 
barrier  is  independent  of  the  tunneling  direction.  In  the  triangular  barrier  range,  only  an 
approximate  form  of  pre-exponential  factor  can  be  derived  (see  Appendix).  A  similar 
factor  was  used  by  Krieger  and  Swanson  [68].  As  it  does  not  affect  the  gate  voltage 
dependence  of  tunneling  current,  we  simply  use  T0=l  in  the  triangular  barrier  range. 

The  band  structures  of  Si02  was  studied  by  Li  and  Ching  [96].  However,  the 
exact  dispersion  relationship  in  the  wide  energy  gap  of  Si02  is  not  known.  We  use  the 
simplest  parabolic  band  approximation  in  (6).  The  empirical  Franz-type  dispersion 
relationship  has  been  used  by  several  authors  [92,94],  which  treated  the  conduction  band 
barrier  and  valence  band  barrier  on  the  same  footing.   However,  it  assumed  equal 
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effective  masses  near  the  valence  band  edge  and  near  the  conduction  band  edge,  i.e., 
mcb=mvb.  This  is  not  supported  by  theoretical  band  structure  calculations  [96]  that 
showed  a  much  smaller  curvature  near  the  valence  band  edge  (mvb>mcb).  Using  the 
form  suggested  by  Freeman  and  Dahlke  [97],  a  more  appropriate  empirical  E(k) 
relationship  in  the  energy  gap  is: 

=    h2/[2mcb(0B-qEoxx-Ex)  ]+n2/{2mvb[Eg-(0B-qEoxx-Ex) ]  }  (4.8) 

where  E  =8.5eV  is  the  energy  gap  of  Si02.  Equation  (4.8),  however,  does  not  yield  an 
analytical  solution  for  the  tunneling  probability.  Since  mvb  is  much  larger  than  mcb,  it 
would  be  a  better  approximation  by  simply  dropping  the  second  term  on  the  R.H.S.  of 
(4.8)  (the  parabolic  band  approximation)  than  retaining  both  terms  and  making  mvb=mcb 
(Frantz-type  dispersion  relationship). 

It  would  be  straightforward  to  prove  that  the  equation  (3.3)  reduces  to  the 
classical  Fowler-Nordheim  formula  under  these  conditions:  (1)  low  temperature:  T— >0, 
(2)  highly  degenerate  n-type  surface  layer:  EF-Ec>kT  (3)  triangular  barrier  and  (4)  unity 
pre-exponential  factor:  T0(EX)=1.  The  classical  F-N  formula  predicts  that  a  semilog  plot 
of  JFN/E0x  vs-  1/Eox  's  a  straight  line: 

JFN   =   aE02    exp(-0/Eox)  (4.9) 

a   =   q3md/ (167T2hmox</>B) 

p   =    4(2mox0B3)1/2/(3qh) 

The  total  conduction  band  electron  tunneling  current  is  obtained  by  summing 

over  contributions  from  all  six  conduction  band  subvalleys  of  silicon.   The  result 


94 

depends  on  silicon  orientation.  For  <100>  Si,  there  are  two  longitudinal  subvalleys  with 
normal  mass  msi=0.916  and  density-of-state  mass  md=0.190,  and  four  transverse 
subvalleys  with  msi=0.190  and  md=0.417.  For  <11 1>  Si,  there  are  six  equivalent 
subvalleys  with  msi=0.258  and  md=0.358  [98]. 

Valence  electron  tunneling  and  hole  tunneling  formulas  are  derived  in  a  similar 
way.  Silicon  valence  band  is  consisted  of  light-hole,  heavy-hole  and  split-off  bands. 
For  simplicity,  we  neglect  the  contribution  from  the  split-off  bands.  The  minimum 
barrier  height  for  valence  electron  tunneling  is  <J)Bn+E  (Si)=3.13eV+1.12eV=4.25eV, 
while  the  maximum  barrier  height  for  valence  hole  tunneling  equals  the  valence  band 
offset  (J)B  =4.25eV.  From  each  of  these  valence  bands,  the  tunneling  current  formula  for 
valence  electrons  (JVe)  ana"  holes  (Jp)  are  listed  as  follows: 


47Tqmd 
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T0(Ex)exp[-2D(Ex)  ] 
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l+exp[ (EV-EF-EX)  /kT] 


dEx  (4.11) 


ll+exp[  (Ev-EF-Ex-qVGB)  /kT] 

In  both  formulas,  Ex  is  referenced  to  the  valence  band  maximum  Ev.  The  valence 
electron  formula  has  the  natural  logarithmic  term  replaced  by  qV0x-Eg-Ex,  because  of 
the  nearly  unity  occupation  factor  of  the  valence  electron  state.  The  upper  limit  of  the 
integration  qV0x_Eg  is  determined  by  constant  energy  tunneling  pathways.  The 
crossing  of  the  conduction  band  and  valence  band  opens  the  constant  energy  tunneling 
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pathways,  which  determines  a  threshold  oxide  voltage  for  the  onset  of  valence  electron 
tunneling:  V0X=EJq.  The  hole  tunneling  formula  is  almost  the  same  as  the  conduction 
band  electron  tunneling  formula  except  for  the  band  edge  reference  and  the  barrier 
height.  The  values  for  the  normal  mass  and  density-of-state  mass  are  derived  from  an 
expansion  of  the  warped  valence  band  E(k)  relationship  around  the  k=0  [99].  They  are: 
msi=0.28  and  md=0.45  for  the  heavy-hole  band  and  msi=0.20  and  md=0.16  for  the  light- 
hole  band. 

The  only  remaining  parameter  in  the  three  tunneling  current  formulas  (4.3), 
(4.10)  and  (4.1 1)  is  the  oxide  tunneling  effective  mass,  mox,  which  can  be  determined 
by  least-squares-fit  (LSF)  of  the  conduction  band  electron  tunneling  formula  (4.3)  to 
Weinberg's  data  of  electron  tunneling  from  <100>  Si  through  a  triangular  barrier  of 
Si02  [66]  shown  as  a  solid  line  in  Figure  4.2.  Weinberg's  F-N  tunneling  current  data 
was  obtained  on  thick  (~100nm)  Si02  samples  under  high  voltage  and  therefore  has  less 
ambiguity  in  the  values  of  oxide  electric  field  than  the  more  recent  thin  oxide  data.  We 
fit  the  theoretical  tunneling  currents  through  a  MOS  structure  with  thin  (5nm)  oxide 
computed  from  (4.3)  to  the  experimental  data  in  the  triangular  barrier  range  at 
Eox>6MV/cm.  The  three  theoretical  curves  with  n-type  impurity  doping  concentrations 
from  1016cm~3  to  1020  cm-3  converge  in  this  high  field  range.  The  tunneling  effective 
mass  thus  determined  is  mox  =  0.437±0.005m0.  This  value  is  used  to  compute  both 
electron  and  hole  tunneling  currents  through  2.0nm  and  5.0nm  oxides.  Notice  that  the 
ratio  of  hole  currents  and  electron  currents  (IN/IP)  increases  as  the  oxide  is  thinned  or  the 
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Fig.  4.2  Calibration  of  tunneling  theory  to  the  experimental  data  (solid  line)  of 
electron  tunneling  from  <100>  Si  through  a  triangular  barrier  of  Si02  in 
reference  [66].  This  calibration  determines  the  value  of  tunneling  effective 
mass  which  is  used  to  compute  the  theoretical  curves  of  electron  and  hole 
tunneling  for  two  gate  oxide  thickness  (2.0nm  and  5.0nm)  and  three 


substrate  doping  concentrations  (1020  cm  3, 


1018cnT3andl016cm-3). 
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oxide  field  is  decreased,  pointing  to  the  importance  of  hole  currents  in  thin  oxide  MOST 
at  low  voltages  which  has  not  been  studied  in  detail. 

4.4  Theoretical  Tunneling  Current- Voltage  Characteristics 

We  will  next  compute  electron  and  hole  tunneling  currents  as  a  function  of  gate 
voltage  in  the  trapezoidal  barrier  range.  At  each  applied  terminal  voltage,  the  voltage 
drops  in  the  gate,  oxide  and  substrate  layers  are  obtained  by  numerically  solving  the 
one-dimensional  MOS  equation  taking  into  account  of  Fermi-Dirac  distribution  and 
impurity  deionization  effects.  The  voltage  drops  in  the  gate  and  substrate  layers 
determine  the  density-of-states  factor  at  the  two  interfaces  while  the  oxide  voltage 
together  with  the  oxide  thickness  determines  the  barrier  transmission  probability. 

Figure  4.3  shows  all  three  constant  energy  tunneling  current  components  in  two 
geometrical  areas  of  a  p+gate  pMOST  with  2.0nm  gate  oxide:  over  the  n-substrate 
channel  area  in  part  (a)  and  over  the  gate-overlapped  p+source  and  p+drain  areas  in  part 
(b).  The  impurity  doping  concentrations  are  those  of  a  typical  pMOST:  5x10  cm-3  in 
p+gate,  1018cm-3  in  n-substrate  and  1019cmf 3  in  p+S  and  p+D.  The  three  current 
components  in  each  part  are:  (1)  conduction  band  electron  currents  JN,  (2)  valence 
electron  currents  JVE  and  (3)  hole  currents  Jp.  There  are  two  distinctive  features  at  low 
voltages.  First,  hole  tunneling  current  Jp  is  dominating  in  both  geometric  areas  which  is 
overtaken  by  JVE  at  large  negative  voltages  and  by  JN  at  large  positive  voltages. 
Secondly,  the  overlap  area  passes  much  higher  hole  current  density  than  the  channel 
area.  The  dominance  of  hole  currents  over  electron  currents  is  due  to  the  low 
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Fig.  4.3  Gate  voltage  dependence  of  oxide  tunneling  current  densities  for  three  types 
of  tunneling  carriers  in  pMOST's  (a)  channel  area  and  (b)  gate-overlapped 
source  and  drain  areas.  The  three  carrier  types  are  (1)  conduction  band 
electrons,  (2)  holes  and  (3)  valence  band  electrons.  Parameters  for  the 
pMOST  are:  dopant  impurity  concentrations:  PGG=5xl019cm-3, 
NBB=1018cm-3  and  PDD=1019cm-3,  and  gate  oxide  thickness  Xox=2.0nm. 
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concentration  of  conduction  band  electrons  at  the  interface  of  p+Si/Si02.  Tunneling  of 
conduction  band  electrons  from  p+Si  could  become  significant  only  when  the  surface  of 
the  p+Si  starts  to  invert,  such  as  those  at  VGB<-2.6V.  The  large  difference  in  Jp  of  the 
two  geometric  areas  is  due  to  the  ~1V  difference  in  gate  voltages  required  for  the 
flatband  condition  (or  Eox=0)  in  the  two  areas.  Flatband  occurs  at  -IV  in  the  channel 
area  and  at  ~0V  in  the  overlap  area.  At  0<VGB<1  V,  Vox<0  in  the  channel  area  which 
impedes  hole  tunneling  from  p+gate,  while  Vox>0  in  the  overlap  area  which  aids  hole 
tunneling.  The  negative  oxide  voltage  lowers  Jp  by  raising  the  integrated  tunneling 
barrier,  or  1)  in  (4.5).  When  the  channel  is  strongly  inverted  at  VGB<-0.4V,  source  and 
drain  areas  are  electrically  connected  by  the  inversion  channel  and  the  large  difference 
in  Jp  of  the  two  geometric  areas  disappears. 

Figure  4.4  compares  the  three  energetic  tunneling  current  components  in  the  two 
geometric  areas  in  a  n+gate  nMOST  with  the  complementary  impurity  type  and  the 
same  doping  concentrations  as  those  of  the  pMOST  in  Fig.  4.3.  Here  conduction  band 
electron  current  JN  is  dominant  over  the  entire  voltage  range  and  in  both  areas  as  the 
concentration  of  conduction  band  electron  is  not  a  limiting  factor  in  n+Gate.  The 
overlap  area  still  passes  much  higher  JN  than  the  channel  area  at  low  voltages  for  the 
same  reason  described  in  the  example  of  PMOST. 

In  the  next  two  figures,  we  illustrate  the  effects  of  substrate  impurity  doping 
concentrations  on  the  tunneling  current-voltage  characteristics  of  a  MOS  capacitor 
(MOSC).  Figure  4.5  shows  JP  in  part  (a)  and  JN  in  part  (b)  for  p+gate  MOSCs,  with 
substrate  doping  concentration  as  a  parameter  which  varies  from  degenerate  p-type  to 
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Fig.  4.4  Gate  voltage  dependence  of  oxide  tunneling  current  densities  for  three  types 
of  tunneling  carriers  in  nMOST's  (a)  channel  area  and  (b)  gate-overlapped 
source  and  drain  areas.  The  three  carrier  types  are  (1)  conduction  band 
electrons,  (2)  holes  and  (3)  valence  band  electrons.   Parameters  for  the 
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nMOST  are:  dopant  impurity  concentrations:  NGG=5xl0     cm 
PBB=1018cm~3  and  NDD=1019cm-3,  and  gate  oxide  thickness  Xox=2.0nm. 
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intrinsic  to  degenerate  n-type.  The  general  lineshape  for  Jp-V  in  (a)  and  JN-V  in  (b)  is 
different,  although  the  oxide  voltage  is  the  same  for  both  electron  and  hole  tunneling. 
This  reflects  the  difference  in  conduction  band  electron  concentrations  at  the  cathode 
and  valence  band  hole  concentrations  at  the  anode.  The  most  distinctive  feature  of  Fig. 
4.5  is  the  hole  current  JP  at  0<VGB<1V  which  changes  its  magnitude  and  slope  with 
substrate  doping  concentrations.  This  feature  is  attractive  for  transistor  diagnosis  as  the 
slope  of  logJ-V  curve  is  independent  of  the  channel  area.  At  this  gate  voltage  range, 
holes  tunnel  from  the  p+gate  to  the  substrate  and  Jp  is  mainly  determined  by  the  oxide 
voltage.  For  intrinsic  substrate,  oxide  voltage  is  nearly  zero  until  the  Fermi-level  in  the 
substrate  reaches  conduction  band  edge  or  valence  band  edge,  this  gives  a  nearly  flat  Jp. 
For  the  three  p-type  substrates,  flatband  voltages  are  close  to  zero.  Higher  p-substrate 
doping  concentration  results  in  less  voltage  drop  in  the  substrate  and  more  voltage  drop 
across  the  oxide  layer  which  increases  the  magnitude  and  slope  of  the  tunneling 
currents.  For  the  three  n-type  substrate,  however,  the  flatband  voltage  VFB  is  close  to 
IV  and  VOx<0  at  VGB<VFB.  Therefore  larger  |Vox|  from  higher  n-substrate  doping 
concentration  would  reduce  the  magnitude  of  tunneling  current  at  VGB<VFB.  At  large 
positive  voltages,  when  the  surface  of  n-substrate  is  strongly  accumulated  or  the  surface 
of  p-substrate  is  strongly  inverted,  the  Jp  curves  would  converge  as  the  oxide  voltage  no 
longer  depends  on  substrate  doping  concentrations.  This  can  be  understood  from  the 
following  approximated  expressions  of  oxide  voltage  drop  in  the  surface  inversion  range 
and  surface  accumulation  range  of  either  p+gate/n-substrate  MOSC  (pMOSC)  or 
n+gate/p-substrate  MOSC  (nMOSC). 
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VOX    ~    VGB    VFB    VSX    VSG 


~   VGB-VSG  strong    inversion   range  (4.12) 

~  VGB-E_/q  strong  accumulation   range  (4.13) 

Vsx  and  VSG  are  voltage  drops  in  the  substrate  and  gate  respectively.  The  knee  points 

in  Jp  curves  at  VGB<0  and  in  JN  curves  at  VGB>0  reflect  the  condition  when  the  Fermi- 
level  in  the  substrate  reaches  the  conduction  band  edge  (JN)  or  valence  band  edge  (Jp). 

Figure  4.6  shows  JN  and  Jp  of  n+gate  MOSCs  for  the  same  variation  of  substrate 
dopant  type  and  concentrations.  The  lineshapes  of  JN  and  Jp  are  complementary  to  Jp 
and  JN  of  p+gate  MOSCs  respectively:  JN  in  Fig.  4.6(a)  is  nearly  a  mirror  image  of  Jp  in 
Fig.  4.5(a)  while  Jp  in  Fig.  4.6(b)  a  mirror  image  of  JN  in  Fig.  4.5(b).  This  feature  is 
expected  because  both  the  carrier  concentration  and  oxide  voltage  for  JN  (Jp)  at  VGB  in 
n+gate  MOSC  are  the  same  as  those  for  Jp  (JN)  at  -VGB  in  p+gate  MOSC.  The 
magnitude  difference  is  from  the  ~leV  difference  in  the  barrier  height  between  CB 
electron  tunneling  and  VB  hole  tunneling. 

Figure  4.7  shows  the  gate  impurity  doping  concentration  dependence  of 
tunneling  currents  in  the  two  complementary  structures:  n+gate/p-substrate  and 
p+gate/n-substrate  MOSCs.  The  lineshape  of  hole  currents  in  part  (b)  is  complementary 
to  that  of  electron  currents  in  part  (a).  As  the  gate  doping  concentration  NGG  is 
increased  from  5xl018cm~3  to  5xl020cm-3  in  the  n+gate/p-substrate  MOSC,  electron 
currents  JN  increases  but  the  responsible  mechanism  varies  with  the  gate  voltage  range. 
At  VGB<0,  it  is  due  to  increased  electron  concentration  at  the  n+gate  while  at  VGB>0,  it 
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Fig.  4.6  Substrate  dopant  impurity  concentration  dependence  of  gate  oxide 
tunneling  currents  in  n+gate  MOS  structures,  (a)  Electron  tunneling  and 
(b)  Hole  tunneling.  NGG=5xl019cnT3,  Xox=2.0nm. 


105 


10' 


CNJ 


CO 

c 

10"4 

0) 

Q 

-•— • 

10"* 

c 

CD 

i_ 

Z3 

10"8 

o 

t=l  I  I  I  I  I  I  I  I  I  I  I 
Electrons 


i  i  i  i  i 


-2-1012 

Gate  Voltage /(1V)         (a) 


i  i  i  I  i  i  i 

Xox=2.0nm 

-3 


Pvy/Nvv=1018cm-3 


Nqg/Pgg  (cm 

5xio20 

5x1 019 

5x101 


p+Gate/n-sub  /jf      i 


1/ 

>! 
>  i 
//i  i 

/  /n+Gate/p-sub- 

,'.l  /  /  /  i  1  i  i  /'  i  i  I  i  i  i 


-3-2-1  0  1  2 

Gate  Voltage  /(1V)         (b) 


Fig.  4.7  Gate  impurity  doping  concentration  dependence  of  tunneling  currents  in 
n+gate/p-sub  and  p+gate/n-sub  MOS  structures,  (a)  Electron  tunneling 
currents,  (b)  Hole  tunneling  currents.  Substrate  impurity  doping 
concentration  is  1018cm-3  and  Xox=2.0nm. 


106 

is  primarily  due  to  the  increase  of  oxide  voltage  through  the  suppression  of  poly- 
depletion  effect.  At  large  positive  voltages,  when  both  n+gate  and  p-substrate  are 
strongly  inverted  such  as  the  case  of  the  NGG=5xl01 8cm-3  curve  at  VGB>2.0V  in  part 
(a),  electron  currents  from  the  two  complementary  MOS  structures  converge.  This  is 
because  the  oxide  voltage  of  the  nSMOSC  in  strong  inversion  as  described  by  (4.12) 
approaches  that  of  the  pSMOSC  in  strong  accumulation  as  described  by  (4.13). 

Figure  4.8  shows  the  oxide  thickness  dependence  of  tunneling  currents  for  the 
two  complementary  MOSC  structures:  pMOSC  in  part  (a)  and  nMOSC  in  part  (b).  Hole 
current  is  dominant  at  low  voltages  in  pMOSCs,  consistent  with  the  result  of  Fig.  4.3. 
In  the  strong  accumulation  gate  voltage  range,  the  ratio  of  hole  current  and  electron 
current  increases  from  0.001  to  0.1  as  the  oxide  is  thinned  from  3.0nm  to  l.Onm.  To  our 
knowledge,  hole  tunneling  theories  have  not  been  reported  in  the  literatures.  The  result 
of  electron  tunneling  from  accumulation  layer  of  a  pMOSC  (Fig.  14  in  [90])  agrees  with 
those  of  Fig.  4.8(a)  at  VGB>1  V,  while  electron  tunneling  from  inversion  layer  of  a 
nMOSC  (Fig.  8  in  [92]  and  Fig.  15  in  [90])  compares  well  with  those  of  Fig.  4.8(b)  at 
VGB>0.  For  example,  the  current  result  gives  JN=0.23  A/cm2  for  a  nMOS  with  2.0nm 
oxide  at  VGB=2.0V,  which  is  between  JN=0.15A/cm2  in  [92]  and  JN=1.0  A/cm2  in  [90]. 

Metal  gate  was  used  in  the  early  stages  of  MOS  technology  development  and  yet 
it  may  be  revitalized  in  the  future  to  replace  polysilicon  gate  in  the  current  MOS 
technology  to  eliminate  gate  depletion  effect  as  well  as  boron  out-diffusion  from  the  p+ 
polysilicon  gate  [100].  In  figure  4.9,  we  compare  electron  tunneling  currents  through 
MOSCs  with  different  gate  materials  including  three  metals  (aluminum,  tungsten  and 
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Fig.  4.8  Oxide  thickness  dependence  of  electron  and  hole  tunneling  currents  in  MOS 
structures,  (a)  p+gate/Si02/n-substrate.  (b)  n+gate/Si02/p-substrate 
nMOS.  Impurity  dopant  concentrations  are  5xl019cm-3  in  the  gate  and 
1018cm~3  in  the  substrate. 
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gold)  in  addition  to  n+  silicon.  The  work  functions  of  the  three  metals  are:  (Al)4.2eV, 
(W)4.57eV  and  (Au)  5.22eV  [101].  At  VGB>0,  electrons  tunnel  from  the  intrinsic 
silicon  substrate.  The  parallel  shift  among  the  three  metal-gate  curves  reflects  the  work 
function  differences,  while  the  slower  rise  of  the  silicon-gate  curve  is  due  to  silicon  gate 
depletion  effect.  At  VGB<0,  electrons  tunnel  from  the  gate.  The  work  functions 
determine  the  tunneling  barrier  height  as  well  as  the  oxide  voltage  and  therefore  affect 
tunneling  current  density  in  a  more  complicated  manner. 

Finally  we  examine  the  theoretical  temperature  dependence  of  electron  tunneling 
currents  in  Fig.  4.10.  The  MOM  structure  with  tungsten  gate  and  tungsten  substrate 
shows  very  weak  temperature  dependence,  which  is  expected  from  our  tunneling  theory. 
At  VGB=3V,  JN  only  increases  about  5%  from  T=-195°C  to  T=100°C.  Therefore  the 
large  temperature  dependence  at  low  voltages  in  MOSC  with  n+silicon  gate  and  p- 
substrate  must  be  from  semiconductor  effects.  In  particular,  the  flatband  voltage  of  the 
MOSC  VFB  =  -(kT/q)loge(NGGPxx/nj2)  shifts  positively  as  the  temperature  increases 
due  to  the  T3/2exp(-E  /2kT)  dependence  of  intrinsic  carrier  concentration  n,.  The 
amount  of  shift  is  0.20V  from  -195°C  to  100°C,  which  changes  V0X(VGB)  in  the  low 
voltage  range  and  gives  the  temperature  dependence.  At  VGB>0.5V,  tunneling  currents 
at  various  temperatures  converge  since  the  oxide  voltage  in  this  strong  inversion  range 
as  described  by  (4.12)  no  longer  depends  on  the  flatband  voltage. 
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Fig.  4.10  Temperature  dependence  of  electron  tunneling  currents  in  MOM 
(W/Si02/W)  and  SOS  (n+Si/Si02/p-Si)  structures.  Oxide  thickness 
XOx=2.0nm.  Temperature  varies  from  100  degree  C  to  -195  degree  C. 
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4.5  Experimental  Correlations 

In  this  section,  correlation  of  the  tunneling  theory  with  experimental  data  on  thin 
oxide  sourced  MOS  capacitors  (SMOSC)  will  be  given.  These  SMOSCs  have  a 
diffusion  region  surrounding  the  gate  area  to  supply  and  collect  minority  carriers  when 
the  surface  layer  of  the  SMOSC  is  inverted.  The  SMOSC  are  better  than  a  regular 
MOSC  in  correlation  with  the  tunneling  theory  because:  (1)  The  measured  terminal 
currents  are  not  limited  by  thermal  generation  of  minority  carriers  in  the  surface 
inversion  layer  of  the  substrate,  (2)  Separation  of  electron  and  hole  currents  are  made 
possible  by  the  additional  source  terminal  and  (3)  It  allows  the  study  of  the  alternative 
tunneling  pathway  through  the  gate-overlapped  source  region.  The  same  overlap  region 
exists  in  the  drain  and  source  diffusion  regions  of  a  MOS  transistor  so  the  results  on 
SMOSC  can  be  readily  applied  to  a  MOST. 

Figure  4.11  shows  the  typical  experimental  current-voltage  characteristics  for 
two  complementary  SMOSCs,  part  (a)  for  the  p+gate  pSMOSC  and  part  (b)  for  the 
n+gate  nSMOSC.  The  DC  bias  configurations  and  the  cross-sections  for  the  two 
SMOSCs  are  shown  on  the  right  side  of  the  experimental  data.  There  are  four  terminals 
in  a  pSMOSC:  gate  (G),  source  (S),  n-well  (B)  and  p-substrate  (C)  and  only  three 
terminals  in  a  nSMOSC  as  both  the  well  and  substrate  are  p-type.  All  terminal  currents 
were  measured  as  a  function  of  the  gate  voltage  between  -IV  and  +1 V.  The  most 
prominent  feature  is  the  complementarity  of  the  gate  currents  in  the  two  SMOSCs:  IG  in 
the  pSMOSC  is  almost  a  mirror  image  of  that  in  the  nSMOSC  around  the  axis  of 
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VGB=0.  There  is  an  inflection  point  in  the  IG-VGB  curve  of  the  pSMOSC  at 
VGB=-0.3V,  which  also  appears  in  that  of  the  nSMOSC  at  VGB=0.3V.  The 
complementarity  strongly  suggests  that  if  CB  electron  tunneling  is  dominant  in  the 
nSMOSC  then  VB  hole  tunneling  must  be  dominant  in  the  pSMOSC.  The  oxide 
voltage  is  too  low  in  this  VGB  range  for  valence  electrons  to  tunnel  elastically  to  the 
conduction  band  on  the  other  side  of  the  Si02  barrier.  However,  valence  electrons  may 
tunnel  to  interface  traps  at  the  other  Si/Si02  interface.  Next,  we  will  correlate  the 
experimental  terminal  currents  to  the  three  energetic  tunneling  pathways  and  the  two 
geometric  tunneling  pathways,  first  for  the  pSMOSCs  and  then  for  the  nSMOSCs. 

In  figure  4.12,  we  breakdown  the  four-terminal  currents  in  the  pSMOSC.  Each 
terminal  current  is  labeled  to  identify  its  energetic  and  geometric  tunneling  pathways  as 
shown  by  the  cross  sections  and  energy  band  diagrams  below  the  experimental  curves. 
Pathways  1  (IG  at  VGB<-0.25V)  and  1'  (Ic  at  VGB>0.3V)  are  hole  tunneling  through  the 
channel  area  (IGBp)  from  the  n-well  and  from  the  p+gate  respectively,  while  pathways  2 
(IG=IS  at  -0.25V<VGB<0V)  and  2'  (IB  at  VGB>0V)  are  hole  tunneling  through  the  G/S 
overlap  (IGSp)  in  the  two  directions  respectively.  IB  measures  electron  tunneling  current 
through  the  channel  area  but  can  not  be  explained  by  any  of  the  three  energetic  tunneling 
pathways  as  will  be  shown  later.  It  corresponds  to  pathways  3  and  4  both  of  which 
involve  Si/Si02  interface  traps.  A  strong  evidence  for  hole  tunneling  is  Ic  at  VGB>0. 
The  p-substrate  collects  holes  that  tunnel  from  the  p+gate  and  then  transverse  the  n-well 
without  recombination.  There  is  a  knee  point  at  0.3V  below  which  Ic  drops  faster.  This 
gate  voltage  corresponds  to  the  onset  of  weak  inversion  at  the  surface  of  n-well.  The 
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inversion  channel  redirects  holes  tunneled  from  p+gate  to  the  source  terminal  and 
therefore  reducing  the  hole  currents  collected  by  the  p-substrate.  Is=Ig  m  tne  entire  VGB 
range.  This  is  because  in  the  strong  inversion  gate  voltage  range,  the  source  terminal 
supplies  the  minority  carriers  (holes)  that  tunnel  through  the  channel  area  and  in  the 
depletion  to  flatband  gate  voltage  range,  the  G/S  overlap  area  is  the  dominant  geometric 
tunneling  pathway.  The  inflection  point  in  IG(Is)-VGB  curve  at  -0.25V  marks  the 
condition  of  equal  tunneling  currents  (IcBp=^GSp)  through  the  two  geometric  areas.  At 
VGB>-0.25V,  IGS  is  larger  than  IGB  despite  its  much  smaller  area  due  to  the  large 
difference  in  the  tunneling  current  densities  (JGBp  and  JGsp)  as  predicted  by  the 
tunneling  theory  in  the  previous  section,  refer  to  figures  4.3  and  4.5.  At  VGB<-0.25V, 
IGB  takes  over  IGS  when  the  difference  between  JGB  and  JGS  becomes  smaller  after 
the  surface  of  n-well  is  inverted.  Next  we  will  quantitatively  verify  the  above  analysis 
using  the  tunneling  theory  developed  in  the  previous  sections. 

Correlation  of  hole  tunneling  theory  with  experimental  data  on  pSMOSC  is 
given  in  figure  4.13.  The  three  sets  of  experimental  data  were  measured  from  three 
wafers  and  the  variation  represents  the  wafer-to-wafer  non-uniformity.  In  part  (a),  IG 
and  Ic  are  normalized  by  the  channel  area  AGB  =  WxL=562u.mx62|im  = 
3.4844xl0_4cm-2.  They  are  matched  by  the  theory  of  hole  tunneling  through  the 
channel  area  (JGBp)  computed  from  (11).  The  best  match  was  obtained  using  impurity 
doping  concentrations  of  Nxx=2xl018cm-3  in  the  n-well  and  PGG=3-5xl019cm"3  in  the 
p+gate  and  an  oxide  thickness  of  1.30nm.  Each  of  the  transistor  parameters  affects  the 
tunneling  J-V  curve  in  a  different  way:  raising  PGG  increases  the  current  density  in  the 
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strong  inversion  range  near  VGB=-1V  but  does  not  change  J  at  VGB=1V;  raising  Nxx 
increases  the  absolute  value  of  the  threshold  voltage  and  shifts  the  knee  point 
negatively;  adjusting  oxide  thickness  gives  a  vertical  shift  of  the  J-V  curve.  These 
features  can  be  used  judiciously  to  extract  transistor  parameters  from  the  tunneling 
current-voltage  characteristics.  The  agreement  between  theory  and  experiment  is 
remarkably  well.  The  theory  computed  from  the  same  formula  matches  both  JG  data  (at 
VGB<-0.3V)  and  Jc  data  (at  VGB>0.5V).  The  deviation  of  theory  and  Jc  data  at 
VGB<-0.5V  was  already  explained  by  the  surface  inversion  channel  in  the  n-well  that 
redirects  hole  currents  to  the  source  terminal.  The  IG  data  at  VGB>-0.3V  is  explained  in 
part  (b)  by  the  theory  of  hole  tunneling  through  the  G/S  overlap  area  (JGsp)-  In  part  (b), 
Is  (=IG)  is  normalized  by  the  G/S  overlap  area  AGS  =  2(W+L)Loverl  =2(562+62)|j.m 
xO.lfxm  =  1.248xl0_6cm-2,  assuming  O.lfxm  for  the  length  of  the  overlap  region.  Two 
theories  were  computed  using  the  same  p+gate  doping  concentration  and  oxide 
thickness  as  in  part  (a)  and  a  p+source  doping  concentration  Pss=PGG=5xl019cirT3.  The 
short-dash  line  is  the  result  of  the  current  tunneling  theory  that  takes  into  account  of  the 
exact  pre-exponential  factor  of  WKB  tunneling  probability  for  a  trapezoidal  barrier, 
refer  to  equation  (4.4),  which  is  in  excellent  agreement  with  the  experimental  data.  The 
long-dash  line  is  the  result  of  a  simplified  tunneling  theory  that  uses  the  traditional 
WKB  tunneling  probability  with  unity  pre-exponential  factor  (T0=l).  The  discrepancy 
between  the  simplified  theory  and  the  experimental  data  becomes  larger  at  lower  |VGB|. 
This  is  expected  from  our  previous  analysis  of  the  unity  pre-exponential  factor 
approximation,  the  largest  error  occurs  at  the  flatband  condition  which  is  at  VGB=0 
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for  the  overlap  area.  At  VGB<-0.3V,  tunneling  current  through  the  channel  area  starts  to 
dominate  and  Is=Iq  is  already  explained  in  part  (a)  by  the  JGB  theory. 

In  figure  4. 14,  we  add  up  two  theoretical  hole  tunneling  currents,  Igbd=Jgbp^gb 
from  the  channel  area  and  Igsd^gSdAgs  from  the  G/S  overlap  area.  The  result 
reproduces  very  closely  the  experimental  IG-VGB  characteristics  of  the  three  pSMOSCs. 
It  clearly  shows  that  the  inflection  point  at  -0.3V  in  the  experimental  data  is  from  the 
crossing  of  IGB  and  IGSp  and  that  the  G/S  overlap  area  is  the  dominant  geometric 
tunneling  pathway  in  the  VGB  range  from  -0.3V  to  1.0V  despite  the  area  ratio  of 
AGS/AGB=l/279!  At  VGB  close  to  1.0V,  the  theoretical  hole  current  can  not  fully 
account  for  the  experimental  IG  data.  This  is  because  the  electron  tunneling  currents  as 
shown  by  IB  in  fig.  4.12  become  comparable  to  the  hole  tunneling  currents.  However, 
the  electron  tunneling  currents  in  the  pSMOSCs  can  not  be  explained  by  any  of  the 
trapless  electron  tunneling  pathways,  as  will  be  demonstrated  in  the  next  figure. 

In  figure  4.15,  we  plot  the  theoretical  electron  tunneling  currents  JGBn  using 
equation  (4.3)  and  valence  electron  tunneling  currents  JGBve  using  equation  (4.10)  for 
two  oxide  thickness  values  of  1.3nm  and  1.6nm.  Experimental  IB  data  are  much  larger 
than  the  theoretical  predictions  for  trapless  electron  tunneling  at  low  voltages.  At 
VGB>0,  due  to  thermionic  barrier  effect,  the  electron  tunneling  (JGBn)  theory  gives  a 
knee  point  near  the  flatband  voltage  ~1V  when  the  conduction  band  edge  Ec  on  the 
anode  (gate)  side  lines  up  with  that  on  the  cathode  (n-well)  side.  Each  of  the 
experimental  IB  curves  shows  a  knee  point  around  0.4V  which  indicates  electron 
tunneling  from  the  conduction  band  of  n-well  into  p+gate/Si02  interface  traps.  This 
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confirms  the  tunneling  pathway  number  4  in  fig.  4.12,  which  is  also  shown  in  fig.  4.15 
below  the  experimental  data.  The  energy  level  of  the  interface  trap  can  be  deduced  from 
the  oxide  voltage  drop  at  the  knee  point,  which  is  about  0.1 8eV  below  Ec.  At  VGB<0, 
JGBn  is  much  smaller  than  the  experimental  IB  data  while  the  trapless  valence  electron 
tunneling  current  JGBve  nas  a  threshold  |VGB|  greater  than  1.7V.  The  theory- 
experimental  discrepancy  suggests  interface  trap  assisted  valence  electron  tunneling 
from  the  p+gate,  which  is  the  tunneling  pathway  number  3  in  fig.  4.12  also  shown  here 
in  fig.  4.15.  IB  is  nearly  constant  at  -0.4V<VGB<0  and  takes  off  at  VGB<-0.4V,  which 
corresponds  to  an  oxide  voltage  of  -0.32V  in  the  pSMOSCs.  This  experimental  feature 
suggests  that  interface  trap  energy  levels  are  more  than  0.32eV  above  the  valence  band 
edge  Ev  at  the  Si02/n-well  interface. 

Following  the  above  approaches  for  the  pSMOSCs,  we  will  next  give  the  theory- 
experimental  correlations  for  the  nSMOSCs.  The  nSMOSC  data  in  Fig.  4.1 1(b)  is 
analyzed  in  detail  in  Fig.  16  to  identify  the  specific  energetic  and  geometric  tunneling 
pathways  in  each  of  the  three  terminal  currents  IG,  IB  and  Is.  In  the  nSMOSC, 
conduction  band  electron  is  the  dominant  tunneling  carrier  species.  Electron  tunneling 
through  the  channel  area  IGBn  is  identified  for  1q=Is  at  VGB>0.25V  (pathway  1)  and  IB  at 
VGB<-0.35V  (pathway  1'),  while  electron  tunneling  through  the  G/S  overlap  area  IGSn 
(pathways  2  and  2')  is  identified  for  IG=IS  at  VGB<0.25V.  Similar  to  the  pSMOSCs,  the 
inflection  point  in  the  IG  curve  at  VGB=0.25  V  is  due  to  the  crossing  of  the  two  geometric 
electron  tunneling  components,  i.e.,  IcBn^GSn'  anc*  tne  ^nee  Pomt  m  tne  ^b  curve  is  due 
to  the  onset  of  (weak)  surface  inversion  in  the  p-well.  All  these  will  be  quantitatively 
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verified  by  the  tunneling  theory.  However,  the  trapless  tunneling  theory  can  not  explain 
IB  in  the  positive  gate  voltage  range,  which  is  likely  due  to  valence  electron  tunneling 
from  the  p-well  to  interface  traps  at  the  n+gate/Si02  interface.  All  the  electron  tunneling 
pathways  in  the  nSMOSC  are  further  illustrated  by  the  cross-sections  and  energy  band 
diagrams  in  figure  4.16. 

Two  sets  of  nSMOSC  data  that  have  different  IG-VGB  lineshapes  will  be  used  to 
correlate  with  the  tunneling  theory.  The  second  set  of  data  are  from  a  previous 
development  technology  with  lower  p-well  doping  concentrations. 

Figure  4.17  gives  the  theory-experimental  correlation  for  three  nSMOSCs  from  a 
recent  technology  which  includes  the  nSMOSC  shown  in  Figs.  4.1 1  and  4.16.  In  part 
(a),  theory  of  electron  tunneling  through  the  channel  area  (JGBn)  *s  compared  to  the 
experimental  data  of  JG=IG/AGB  anc*  Jb=^b^gb-  Trie  best  matches  shown  here  are 
obtained  by  using  impurity  doping  concentrations  of  5xl019cm-3  for  the  n+  gate  and  1- 
2xl018cm-3  for  the  p-well  and  an  oxide  thickness  of  1.60nm.  The  three  theoretical 
curves  correspond  to  p-well  impurity  doping  concentrations  (xl018cm~3)  of  1.0,  1.4  and 
2.0  respectively.  The  first  curve  (1.0)  fits  the  depletion-flatband  range  (VGB<-0.3V) 
well,  the  third  curve  (2.0)  fits  the  inversion  range  (VGB>0.5V)  well,  while  the  second 
curve  (1.4)  gives  the  best  overall  fit.  The  faster  drop  of  experimental  JB  data  towards 
zero  in  the  negative  VGB  range  is  expected  because  of  the  surface  inversion  channel  that 
redirects  electrons  tunneled  from  the  n+gate  to  the  source  terminal,  which  reduces  the 
number  of  electrons  being  injected  into  the  p-well/p-substrate.  IG  at  VGB<0.3V  is 
dominated  by  electron  tunneling  through  the  G/S  overlap  region,  as  proved  by  theory- 
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with  JS=IS/AGS.  The  tunneling  theories  in  the  two  parts  use  the  same 
n+gate  doping  concentration  of  5xl019cm-3  and  oxide  thickness  of  1.60nm. 
Three  theoretical  JgBn  curves  in  part  (a)  correspond  to  p-well  doping 
concentrations  (xlO  cm-3)  of  1.0,  1.4  and  2.0  respectively,  and  the  two 
theoretical  JGsn  curves  in  Part  (b)  corresponds  to  n+gate  doping 
concentrations  (xl019cm~3)  of  3.0  and  5.0  respectively. 
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Fig.  4.18  Two  geometrical  tunneling  pathways  in  n+sourced  n+Si/Si02/p-Si  MOS 
capacitors.  Experimental  Ig_Vgb  curves  are  matched  by  adding  up  two 
theoretical  electron  tunneling  currents,  IGBn  for  electron  tunneling  through 
the  channel  area,  and  IGSn  for  electron  tunneling  through  the  G/S  overlap 
area. 
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experimental  correlation  of  JS=IS/AGS  in  part  (b).  The  same  oxide  thickness  and  n+gate 
doping  concentration  in  part  (a)  together  with  two  n+  source  doping  concentrations 
(xl019cm~3)  of  3.0  and  5.0  were  used  in  the  computation  of  the  two  theoretical  electron 
tunneling  curves  in  part  (b).  Noticeably,  unlike  the  case  of  pSMOSC  in  which  the 
match  is  almost  perfect,  the  experimental  Js  data  in  nSMOSC  exhibits  faster  roll-off  in 
the  low  voltage  range  than  the  theory  even  with  the  refined  WKB  tunneling  theory.  This 
points  to  additional  differences  in  electron  and  hole  tunneling  which  are  not  considered 
by  the  current  theory. 

In  figure  4.18,  the  experimental  gate-terminal  current  IG  is  shown  to  have  two 
geometric  tunneling  pathways,  one  through  the  channel  area  described  by  the 
^GBn^GBn^GB  theory  and  the  other  through  the  G/S  overlap  area  described  by  the 
^GSn^^GSn-^GS  theory.  Similar  to  fig.  14  for  pSMOSCs,  this  theory-experimental 
correlation  verifies  that  the  inflection  point  in  IG  at  0.25V  is  from  the  crossing  of  two 
geometric  tunneling  current  component  and  that  the  small  G/S  overlap  area  is  the 
dominant  tunneling  pathway  in  the  low  voltage  range  from  -1.0V  to  0.25V. 

The  second  set  of  nSMOSC  data  shown  in  Figs.  4.19-20  has  the  following 
features  which  are  different  from  those  of  the  previous  example:  (1)  there  is  no 
inflection  point  in  the  gate  current,  (2)  the  ratio  of  IB/IG  in  the  depletion-flatband 
voltages  is  about  3  times  larger,  and  (3)  the  knee  point  in  IB  shifts  more  than  0.2V  in  the 
negative  direction.  All  these  changes  can  be  explained  by  a  lower  p-well  impurity 
doping  concentration  of  5xl017cm~3  as  compared  to  l-2xl018cm~3  which  is 
demonstrated  in  the  two  figures.  In  figure  4.19,  IG  at  VGB>0  and  IB  at  VGB<-0.55V  are 
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Fig.  4.19  Another  theory-experimental  correlation  for  electron  tunneling  in 
nSMOSCs  with  lower  p-well  doping  concentrations  than  that  of  the 
previous  three  nMOSCs  shown  in  Figs.  17-18.  (a)  Comparison  of  the  JGBn 
theory  with  Jg=Wagb  anc*  Jb=WAgb-  (°)  Comparison  of  the  JGSn  theory 
with  Js-Is^gs-  The  tunneling  theories  in  the  two  parts  use  the  same 
transistor  parameters  including  doping  concentrations  of  5xl017cm-3  in  the 
p-well,  5xl019cm-3  in  the  n+  gate  and  2xl019cm-3  in  the  n+  source,  and  an 
oxide  thickness  of  1.52nm. 
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Fig.  4.20  Two  geometrical  tunneling  pathways  in  the  two  nSMOSCs  with  relatively 
low  p-well  doping  concentrations.  Experimental  Ig-^gb  curves  are 
matched  by  adding  up  two  theoretical  electron  tunneling  currents,  IGBn  for 
electron  tunneling  through  the  channel  area,  and  IGSn  for  electron  tunneling 
through  the  G/S  overlap  area. 
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matched  by  the  theory  of  electron  tunneling  through  the  channel  area  in  part  (a),  and 
IS=IG  at  VGB<0  is  matched  by  the  theory  of  electron  tunneling  through  the  G/S  overlap 
area.  Transistor  parameters  that  give  the  best  match  are:  impurity  doping  concentrations 
of  5xl017cm~3  in  the  p-well,  5xl019cm-3  in  the  n+  gate  and  2xl019cm~3  in  the  n+ 
source,  and  oxide  thickness  of  1.52nm.  The  lower  p-well  doping  concentration 
increases  the  electron  tunneling  current  from  n+gate  to  the  p-well  (IB)  by  reducing  the 
integrated  tunneling  barrier,  refer  to  figure  4.5(a).  The  knee  point  is  determined  by  the 
threshold  voltage  for  surface  (weak)  inversion  in  the  p-well  which  shifts  negatively  with 
lower  p-well  doping  concentrations.  The  lack  of  knee  point  in  IG  can  be  readily 
understood  from  figure  4.20  in  which  the  two  geometric  tunneling  current  components 
^GBn  and  ^GSn  are  separated.  Due  to  the  increase  of  the  channel  component  at  low 
positive  voltages,  IGBn  and  IGSn  no  longer  crosses  at  the  lower  p-well  doping 
concentration.  The  overlap  area  component  still  dominates  the  gate-terminal  current  in 
the  negative  gate  voltage  range. 

4.6  Summary 

Followed  by  an  extensive  review  of  the  literatures,  we  developed  a  theory  that 
applies  to  the  tunneling  of  electron,  hole  and  valence  electrons  through  a  trapezoidal 
barrier  of  the  technologically  important  MOS  structure.  The  new  theory  refines  the 
traditional  WKB  tunneling  probability  by  the  incorporation  of  a  kinetic  energy 
dependent  pre-exponential  factor  which  is  important  for  accurate  modeling  of  MOS 
tunneling  currents  at  low  fields.  The  theory  was  applied  to  study  the  tunneling  currents 
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from  all  terminals  of  thin  oxide  sourced  MOS  capacitors  in  order  to  identify  their 
specific  tunneling  carrier  species  and  geometric  tunneling  pathways.  The  experimental 
data  on  MOS  capacitors  verifies  two  conclusions  predicted  by  the  theory  which  also 
applies  to  cMOS  transistors:  (1)  hole  currents  dominate  over  electron  currents  in  p+gate 
pMOS  transistors,  and  (2)  the  gate  overlapped  source  and  drain  areas  pass  higher 
currents  than  the  channel  area  in  the  depletion-flatband  gate  voltage  range  despite  its 
small  area.  Application  of  the  tunneling  theory  to  extract  transistor  design  parameters 
such  as  doping  concentrations  and  gate  oxide  thickness  was  also  demonstrated. 


CHAPTER  5 

INTERBAND  TUNNELING  NEAR  THE  DRAIN/SUBSTRATE 

JUNCTION  PERIMETERS  IN  MOS  TRANSISTORS 

5.1  Introduction 

The  miniturization  of  MOS  transistors  requires  the  use  of  thinner  gate  oxide  and 
higher  surface  impurity  doping  concentration  in  each  new  technology  generation  in 
order  to  minimize  short  channel  effects  in  scaled  transistors  [102].  As  a  result,  the 
drain/substrate  and  source/substrate  junction  space  charge  regions  become  thinner. 
During  transistor  operation,  they  can  be  further  thinned  by  an  accumulation  gate  voltage 
which  makes  the  depletion  of  substrate  majority  carriers  more  difficult  at  the  junction 
perimeter.  It  is  well  known  that  interband  p/n  junction  tunneling  could  dominate  the 
current-voltage  characteristics  if  the  transition  region  is  thinner  than  about  100A,  as  in 
an  Esaki  diode  [56].  The  interband  p/n  junction  tunneling  at  reverse  drain/substrate  bias 
was  recognized  as  a  major  off-state  leakage  current  path  [103].  DCIV  [5,8]  data  on 
recent  generations  of  MOS  transistors  show  a  rising  baseline  in  the  accumulation  gate 
voltage  range  of  the  Ib-^gb  curve>  which  was  leveraged  by  Sah  as  a  sensitive  indicator 
of  surface  impurity  concentrations  [8].  This  experimental  feature  prompted  us  to  look  at 
the  theoretical  gate  voltage  VGB  dependence  of  the  drain  junction  perimeter  tunneling 
current  IPN,  which  is  presented  in  section  5.2.  We  then  examine  the  effect  of  energy 
quantization  in  the  surface  accumulation  layer  on  IPN  and  the  possibility  of  probing 
discrete  energy  levels  from  IpN"VGB  characteristics  in  section  5.3.  A  similar  idea  was 
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proposed  by  Quinn  et  al.  [104]  to  probe  energy  levels  in  the  surface  inversion  layer  of  a 
p-type  substrate  by  measuring  the  tunneling  current  between  a  p+source  and  the  n- 
inversion  layer.  Detailed  calculation  of  drain/substrate  barrier  height  at  the  surface 
reveals  that  the  direct  interband  tunneling  under  forward  junction  bias  has  a  large 
threshold  |VGB|,  not  compatible  with  experimental  DCIV  data  shown  in  section  5.4.  The 
IB-VGB  data  shows  peaked  structures  superimposed  on  the  rising  IB  baseline  in  the 
accumulation  VGB  range.  In  section  5.5,  the  IB  peak  is  correlated  to  Shockley-Read- 
Hall  (SRH)  recombination  at  interface  traps  (NIT)  located  over  the  gate-overlapped  drain 
region.  The  evidence  of  NIT  in  the  overlap  region  led  us  to  develop  a  theory  for  trap 
assisted  tunneling  which  also  has  a  lower  threshold  |VGB|  to  account  for  the  DCIV  data, 
in  section  5.6.  This  tunneling  current  is  similar  to  the  'excess  current'  in  a  bulk  Esaki 
diode,  where  bulk  traps  are  involved  in  carrier  transport,  such  as  the  gold  centers  in  a 
gold-doped  silicon  p/n  diode  studied  by  Sah  [58].  Section  5.7  concludes  this  chapter. 

5.2  Theory  of  Drain  Junction  Perimeter  Tunneling 

Figure  5.1  depicts  the  energy  band  diagram  of  a  nMOS  transistor  with  its 
drain/substrate  junction  forward-biased  at  VpN,  such  as  in  drain-emitter  DCIV 
measurements.  The  barrier  height  VB  of  the  n+/p  junction  is  increased  by  a  hole 
accumulation  surface  potential  Vsx  induced  by  the  gate  voltage  VGB.  If  VB>E  /q, 
conduction  band  electrons  on  the  n+D  side  can  tunnel  into  valence  band  holes  on  the  p- 
substrate  side,  giving  rise  to  the  junction  perimeter  tunneling  current  in  base  and  drain 
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Fig.  5.1  Cross-section  and  energy  band  diagram  of  a  n+Drain/p-Base  junction  of  a 
MOS  transistor  at  X=0  under  strong  hole  accumulation  gate  bias.  Electrons 
at  the  n+  side  with  kinetic  energy  Ee  tunnel  into  hole  states  with  kinetic 
energy  Eh.  VB  =  n+/p  barrier  height,  VBI=  built-in  potential  of  the  n+/p 
junction,  Vsx  =  surface  band  bending  in  the  hole  accumulation  layer  with  a 
thickness  of  xacc,  and  VPN  =  forward  bias  voltage  applied  to  the  n+/p 
junction. 
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terminals.  Energy  conservation  requires  that  the  electron  kinetic  energy  of  the  initial 
state  Ee  and  hole  kinetic  energy  of  the  final  state  Eh  satisfy  the  following  equation 

Ee    +    Eh    =    qVB    -    Eg  (5.1) 

The  general  approach  used  in  chapter  4  for  the  derivation  of  gate  oxide  tunneling 
theory  also  applies  to  the  interband  p/n  junction  tunneling  at  the  surface  with 
distinctions  in  (1)  tunneling  effective  mass  and  (2)  shape  of  the  tunneling  barrier.  Since 
the  initial  electron  state  is  in  the  silicon  conduction  band  near  k=0.85(27i/a),  where 
a=5.43A  is  the  lattice  constant  of  silicon  crystal,  while  the  final  hole  state  is  in  the 
valence  band  near  k=0,  strict  conservation  of  transverse  (perpendicular  to  the  tunneling 
direction)  momentum  kt  must  involve  phonons.  The  theory  of  phonon-assisted 
tunneling  was  developed  by  Kledysh  [105]  and  Price  and  Radcliffe  [106],  and  was 
applied  to  p/n  tunnel  junction  by  Kane  [57].  The  theory  is  rather  laborious,  and  only 
valid  for  constant  electric  field  in  the  barrier  region.  In  our  case  of  p/n  junction  at  the 
surface,  interface  roughness  scattering  provides  another  means  of  conserving  transverse 
momentum.  To  derive  the  essential  features,  we  shall  relax  the  requirement  of  kt 
conservation.  The  resultant  formula  for  interband  tunneling  with  the  initial  state  in  one 
of  the  six  silicon  conduction  band  valleys  is  [107]: 


47rnqmd 
Jt 


'PN 

h3 


qVB-Eg 

T(E) [f (E)-f (E+qVPN)]EdE  (5.2) 

0 


The  integration  is  over  electron  kinetic  energy  E  that  have  a  direct  tunneling  pathway. 
f(E)  is  the  Fermi  function  that  describes  electron  occupation  of  the  initial  state  at  the  n+ 
side  while  f(E+qVPN)  is  the  electron  occupation  factor  of  the  final  state  at  the  p  side.  r\ 
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and  md  are  respectively  the  degeneracy  factor  and  the  density-of-states  effective  mass  of 
the  subvalley.  Assume  x,  y  and  z  are  all  along  <100>  directions,  which  is  consistent 
with  MOS  transistor  manufacturing  practice,  the  four  transverse  valleys  would  give  the 
most  contribution  to  tunneling  because  of  their  lower  effective  mass  in  the  tunneling 
direction  (y).  Therefore  r)=4  and  md=(mxmz)1/2  =  (mtm,)1/2  =  (0.9 16x0. 191  )1/2m0  = 
0.418m0.  T(E)  is  the  barrier  transmission  probability  which  has  the  following 
expression  from  WKB  approximation: 

y2 


T(E)    =    exp 


K(y)dy 
JYl 


(5.3) 


K  is  the  imaginary  part  of  the  wave  number  in  the  energy  gap  of  silicon,  Y  is  the 
tunneling  direction,  and  yj  and  y2  are  the  boundaries  of  the  classical  forbidden  region. 
There  is  no  potential  discontinuity  at  the  two  boundaries,  thus  the  unity  pre-exponential 
factor  is  a  good  approximation.  Notice  that  in  the  forbidden  region,  in  addition  to  the 
decaying  wave  in  y  direction,  there  are  also  the  propagating  waves  in  the  perpendicular 
x  and  z  directions. 

The  simplest  description  of  K  is  a  parabolic  E(k)  relationship.  A  more  rigorous 
form  of  of  k  can  be  obtained  by  analytical  continuation  of  the  E(k)  relationship  in  the 
energy  band  into  the  energy  gap.  The  result  is  the  a  non-parabolic  formula  which 
reduces  to  parabolic  relationship  at  both  conduction  band  edge  and  valence  band  edge, 
which  was  used  by  Kane  [108]  for  InSb,  and  also  by  Maserjian  [109]  for  Si02.  The  two 
alternative  E(k)  relationships  are: 
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h^-h^-h^k*  =    2m  •  [qV(y)-E]  (parabolic)       (5.4) 

=  2m*. [qV(y)-E] . {l-[qV(y)-E]/Eg} 

(non-parabolic)   (5.5) 

m*  is  the  tunneling  effective  mass,  m*  =2  mcmv/(mc+mv)  [107],  where  mc=mt=0.191m0 
is  the  normal  mass  of  the  four  transverse  subvalleys,  and,  mv,  the  light  hole  mass  near 
k=0  is  nyCA+B)  =  mo/(4.28+0.75)  =  0.199m0  [110].  The  result  is  m*=0.195m0.  V(y) 
is  the  electro-static  potential  in  the  gap  region,  referenced  to  the  conduction  band  edge  at 
the  n+  side.  Eg  is  the  indirect  energy  gap  of  1.12eV  in  Si.  In  the  non-parabolic  formula, 
k=0  defines  the  two  end  points  y{  and  y2  of  the  tunnel  region.  It  can  be  seen  that  for 
non-zero  kx  and  kz,  the  tunneling  distance  and  effective  barrier  height  are  both 
increased,  which  results  in  exponential  drop  of  tunneling  probability  [57,111].  For 
simplicity,  we  only  consider  initial  states  with  nearly  zero  transverse  momentum  that 
has  the  largest  tunneling  probability,  i.e.,  kx=kz=0  in  (5.4)  and  (5.5). 

For  a  parabolic  potential  barrier  (or  linear  field),  only  the  parabolic  formula  (5.4) 
can  be  integrated.  Under  depletion  approximation  and  |E  |»|EX|  at  the  surface  for  the 
one-sided  n+/p  junction,  the  potential  on  the  p-side  is: 

V(y)     =     (qNAA/2esi)     [d2-(d-y)2]  (5.6) 

NAA  is  the  doping  concentration  at  the  surface  of  p-substrate,  and  d  is  the  depletion  layer 
thickness  in  y.  Using  (5.4),  we  derive 

T(E)     =    exp{-l|  (2m*/h2)     (2£si/q2NAA)      [  (qVB-E)  ( 7T/2    - 

siiTH  (qVB-E-Eg)  /  (qVB-E)  ]     -    ll  Eg  (qVB-E-Eg)  ]} 

(linear    F,    parabolic)        (5.7) 

VB  is  the  barrier  height  controlled  by  VGB. 
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For  the  simpler  constant  field  case,  i.e.,  V=qFx,  where  F  denotes  the  constant 
field,  T  is  not  a  function  of  E  since  tunneling  electrons  at  different  energy  levels 
experience  the  same  triangular  barrier.  Analytical  results  for  T  are  obtained  for  both 
E(k)  relationships: 

3/2 


T(E)     =   T   =    exp[-aNm  EgJ/V  (qhF)  ] 


-al 

a   =    4N2/3    =  1.886   (constant  F,  parabolic)       (5.8) 

a   =    7T/(2!|2)  =  1.111   (constant  F,  non-parabolic)   (5.9) 

If  the  same  effective  mass  is  used  in  (5.8)  and  (5.9),  the  non-parabolic  formula  (5.9) 

gives  an  effectively  lower  tunneling  barrier  and  higher  tunneling  probability. 

As  a  first  order  approximation,  barriers  of  any  shape  can  be  approximated  by  a 
linear  barrier  with  a  slope  equal  to  the  average  field  F      .   For  a  parabolic  barrier 
described  by  (5.6),  the  average  field  in  the  tunnel  region  (yj<y<y2)  is: 
Favg   =    [F(yi)+F(y2)]/2 

=    tInaa/2£s1    (IqVe   +    ^lqVB-Eg)  (5.10) 

Figure  5.2  gives  the  theoretical  Jpn'Vgb  dependence  for  a  gated  n+/p  junction 
using  non-parabolic  formula  (5.5)  with  triangular  barrier  approximation  (5.10)  in  part 
(a),  and  parabolic  formula  (5.4)  with  the  more  exact  parabolic  barrier  (5.6)  in  part  (b). 
The  tunneling  effective  mass  in  part  (b)  is  scaled  to  that  of  part  (a)  according  to  eqs. 
(5.8)  and  (5.9)  by  a  factor  of  (1.11 1/1. 886)2=0.347.  The  parabolic  barrier  gives  lower 
tunneling  current  density  than  the  triangular  barrier,  due  to  the  T(E)  dependence  in  (5.7), 
which  drops  sharply  with  increasing  electron  kinetic  energy.  At  an  oxide  thickness  of 
1.5nm,  under  large  negative  VGB,  the  interband  tunneling  current  JPN  is  about  4-5  orders 
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Fig.  5.2  Gate  voltage  dependence  of  surface  p/n  junction  perimeter  tunneling  current 
density,  calculated  from  two  theories:  (a)  average  field  approximation  with 
non-parabolic  E(k)  relationship  (b)  linear  field  (parabolic  potential  barrier) 
with  parabolic  E(k)  relationship.  Junction  forward  bias  changes  from  lmV 
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to  200mV  as  labeled.  NAA  =  5xl018  cm"J  in  p-base,  NDD  =  5xl0iy  cm 
n+Gate  and  gate  oxide  thickness  X0x  =  1.5nm. 
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of  magnitude  smaller  than  the  gate  oxide  tunneling  current  J0x  [see  Fig.  4.8(b)],  for  a  p- 
substrate  doping  concentration  of  5xl018cm-3.  The  threshold  |VGB|  increases  from  1.5V 
to  3.3V  as  the  junction  forward  bias  is  raised  from  lmV  to  200mV.  Due  to  this 
threshold  voltage  change,  at  a  fixed  VGB,  JPN  decreases  as  VPN  increases  except  for  the 
two  dashed  lines  at  VPN=lmV  and  lOmV.  The  drop  of  JPN  at  VPN<10mV  is  due  to 
increased  electron  tunneling  in  the  reverse  direction  from  p-sub  side  to  n+D  side.  Two 
more  theoretical  figures  are  given  which  were  computed  from  the  non-parabolic 
formula. 

Figure  5.3  shows  the  oxide  thickness  (Xox)  dependence  of  JpN-VGB 
characteristics.  The  threshold  |VGB|  increases  almost  linearly  with  Xox,  which  is 
expected  from:  AVGB  =  AVox  in  the  strong  accumulation  gate  voltage  range.  At  all 
oxide  thickness,  JPN  is  insignificant  as  compared  to  Jox  and  thus  its  experimental 
detection  will  be  masked  by  the  gate  oxide  tunneling  current.  JPN,  however,  is  a  strong 
function  of  p-substrate  doping  concentration  (NAA)  while  Jox  is  not.  Figure  5.4  shows 
that  JPN  would  overtake  Jox  at  NAA>1019cm-3,  for  Xox=2.0nm  and  VPN=50mV.  Jox  is 
calculated  using  the  trapezoidal  Si02  tunneling  formula  (4.3),  which  is  nearly 
independent  of  NAA  at  -4V<VGB<-2V. 

5.3  Tunneling  into  Quantized  Surface  Accumulation  Layer 

The  proceeding  JPN  formulas  assumed  that  both  the  initial  state  and  final  state  of 
the  tunneling  transition  can  be  described  by  three-dimensional  (3D)  states  characterized 
by  crystal  momentum  (kx,ky,kz).  It  was  recognized  in  1957  by  Schrieffer  [112]  that 
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Fig.  5.3  Theoretical  surface  p/n  junction  tunneling  current  density  as  a  function  of 
gate  voltage,  for  various  oxide  thickness  from  1.5nm  to  3.0nm.  Forward 
bias  VpN=100mV,  p-Base  and  gate  doping  concentrations  are  the  same  as 
those  in  Fig.  5.2. 
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strong  space-charge  field  due  to  a  large  bias  (VGB)  would  confine  the  motion  of  carriers 
in  a  thin  surface  layer  and  quantize  the  energy  levels  in  the  direction  perpendicular  to  the 
surface.  This  quantization  effect  has  been  studied  for  both  surface  inversion  layers 
[98,113]  and  surface  accumulation  layers  [114]  in  MOS  structures.  In  this  section,  we 
will  develop  a  theory  for  surface  electron  tunneling  from  3D  states  in  n+drain  to  2D 
states  in  the  hole  accumulation  layer  of  p-substrate,  and  evaluate  the  effect  of  quantized 
energy  levels  on  JpN-VGB  characteristics. 

In  the  quantized  hole  accumulation  layer,  hole  kinetic  energy  is  described  by: 
Eh    =    En    +     (fi2/2mv)  (kj    +    kl)  (5.11) 

=    qVB    -    Eg    -    E 
where  En  denotes  the  discrete  energy  level  in  the  x-direction,  and  E  is  the  electron 
kinetic  energy  in  the  initial  3D  state.  The  density-of-states  in  the  2D  hole  accumulation 
layer  is: 


2A          27Tk  dk  A  rt^ 
S(En)    = 


(27T)2  dEh  7Th2 


where   S(Eh)    =0,  0     <    Eh   <    E0 

=    n,  En_x    <    Eh   <    En       (n=l,     2.  .  .  ) 

A  is  the  area  of  the  2D  layer,  E0  is  the  energy  level  of  the  ground  state,  and  En  is  the 

energy  level  of  the  n-th  excited  state.  S(Eh)  is  a  staircase  function  with  increasing  step 
of  unity  at  Eh=En.  The  usual  3D  density  of  hole  states  neglecting  quantization  effect  in 
x-direction  is: 
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2V 


47Tk2dk 


D3(Eh)    = 


V  ihmX 


(5.13 


(2n)2  dEh  ;r2h3 


V  si  the  volume  of  the  3D  system  and  V/A  =  Xacc,  a  characteristic  thickness  of  the 
surface  accumulation  layer. 

Since  the  tunneling  probability  is  directly  proportional  to  the  density  of  states  of 
the  final  state,  we  derive  a  JpN  formula  for  tunneling  into  2D  states  from  (5.2),  (5.11), 
(5.12)  and  (5.13): 


4nrjqimxin2 


'PN 


qvB-Eg 

0 


T(E)  [f  (E)-f  (E+qVPN)  ]E    x 
;m/xacc 
T|2rtv(qVB-E  -E) 


S(Eh) 


dE 


(5.14) 


The  additional  factor  in  [  ]  equals  D2(Eh)/D3(Eh).  The  staircase  function  S(Eh)  in  (5.14) 
predicts  break  points  in  Jpn_^gb  curves>  which  corresponds  to  the  crossing  of  n+  drain 
conduction  band  minimum  and  discrete  hole  energy  levels  on  the  p-side. 

We  will  next  use  the  simple  triangular  well  approximation  to  estimate  the 
discrete  energy  levels  in  the  surface  accumulation  layer.  Stern  showed  that  this  simple 
approximation  gives  good  analytical  results  for  surface  inversion  layers  [98].  However, 
using  the  surface  electric  field  Fs  as  the  slope  of  the  triangular  well  would  give  a  ground 
state  level  E0  greater  than  the  depth  of  the  well  given  by  the  surface  band  bending  in  the 
accumulation  layer,  which  is  not  self-consistent.  This  is  because  the  accumulation  layer 
has  a  much  smaller  surface  band  bending  than  the  inversion  layer  due  to  the  lack  of  a 
depletion  layer.  We  will  use  a  slope  which  is  scaled  down  from  Fs.  A  scaling  factor  of 
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0.03  would  get  the  ground  state  energy  level  in  the  range  of  those  obtained  from  more 
rigorous  numerical  solutions  of  Ando  [1 14]  and  Sune  et  al.  [115].  This  gives  the 
following  expression  [98]  for  the  energy  levels: 

En  =  (h2/2mv)1/3»  [  (3;rq/2)  (0.03FS)  (n+3/4)  ]2/3  (  11=0,1,2...)  (5.15) 
The  accumulation  layer  thickness  can  be  estimated  from  the  x-distance  of  ground  state 
holes  from  the  Si/Si02  interface  (X=0)  as  follows  [98]: 

Xacc    =    2E0/[3q(0.03Fs)]  (5.16) 

Figure  5.5  gives  an  example  of  the  2D  density  of  hole  states  normalized  by 
(mv/7ih2)  in  the  surface  accumulation  layer  at  Fs=106V/cm.  Figure  5.6  shows  gate 
voltage  dependence  of  the  lowest  four  bound  state  energy  levels  (solid  lines)  and  the 
surface  band  bending  qVsx  (dashed  line).  Vsx  was  calculated  using  Fermi-Dirac 
statistics  in  gate  and  drain  layers.  Since  q|Vsx|  indicates  the  depth  of  the  surface 
potential  well,  the  four  energy  levels  are  all  below  q|Vsx|  and  therefore  are  consistent 
with  one-dimensional  bound  states. 

We  will  next  give  a  theoretical  illustration  of  quantization  effects  on  surface  p/n 
junction  tunneling  characteristics.  JPN  and  its  derivative  dJPN/dVGB  are  plotted  in  Fig. 
5.7(a)  and  Fig.  5.7(b)  respectively.  The  solid  lines  were  computed  from  (5.14)  using  the 
2D  DOS  function  of  (5.12),  with  Sj=l,  S2=2  and  S3=3.  The  dashed  lines  are  from  two 
hypothetical  2D  DOS  functions:  S1=S2=S3=1  for  the  long-dash  lines  and  S0=l,  Sj=2, 
S3=2  for  the  short-dash  lines.  Two  breakpoints  at  VGB=  -2.79V  and  -3.73V  in  the  solid 
line  can  be  readily  observed  from  the  derivative  plot  in  Fig.  5.7(b),  also  visible  in  Fig. 
5.7(a)  with  the  help  of  the  dash  lines.  These  break  points  correspond  to  the  crossing 
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DOS  in  hole  accumulation  layer 

Trianglar-well  appr.  Eeff=  0.03  Esx 
Esx  =  1MV/cm,  mv=mLH 
Normalized  to  (Am/7ih2) 
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Fig.  5.5  Normalized  density  of  hole  states  in  a  quantized  suface  accumulation  layer 
with  a  surface  electric  field  Esx  =  lMV/cm.  Quantum  leaps  occur  when 
hole  energy  Eh=  nth  energy  level  of  one-dimensional  hole  bound  states. 
The  light  hole  band  normal  effective  mass  mv  =  0.1991110  was  used. 
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Fig.  5.6  Gate  voltage  dependence  of  the  four  lowest  energy  levels  of  one- 
dimensional  hole  bound  states  in  a  surface  accumulation  layer  with  the  band 
bending  shown  as  the  dash  line.  Same  theoretical  model  as  in  Fog.  5.5  for 
the  discrete  energy  levels.  NAA=1019cm"3  in  p-base  and  Xox=1.5nm. 
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Fig.  5.7  Effects  of  quantized  hole  accumulation  layer  on  surface  p/n  junction 
tunneling  current,  (a)  Tunneling  current-  gate  voltage  characteristcs  and  (b) 
its  first  derivative.  The  solid  lines  were  computed  from  theoretical  stair- 
case 2D  density  of  state  in  (5.12),  while  the  dashed  line  used  two 
hypothetical  DOS  functions.  The  two  breakpoints  due  to  the  discontinuity 
of  DOS  function  at  Eh=Ej  and  Eh=E2  are  more  obvious  in  part  (b). 
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of  the  n+  side  conduction  band  edge  with  the  first  excited  and  second  excited  hole  state 
energy  levels  on  the  p-substrate  side  respectively. 

We  conclude  this  section  by  noting  that  direct  surface  p/n  junction  tunneling  is 
not  significant  in  current  MOS  technologies  with  substrate  doping  concentrations  on  the 
order  of  1017-1018cm~3.  In  samples  with  high  substrate  doping  concentration  so  that  JPN 
dominates  Jox,  it  may  De  possible  to  gauge  quantized  energy  levels  in  the  surface 
accumulation  layer  from  the  breaks  in  the  JpN-VGB  characteristics.  High-order 
derivatives  may  be  needed  to  enhance  the  sensitivity  to  detect  the  quantum  effects  under 
finite  temperatures  since  the  one-dimensional  discrete  energy  levels  would  be  broadened 
to  yield  a  two-dimensional  DOS  with  no  sharp  steps. 

5.4  Base-Terminal  Current  Under  Accumulation  Gate  Voltages 

Recent  experimental  DCIV  data  on  thin  oxide  (<  5nm)  MOS  transistors  [5,116] 
show  base-terminal  current  IB  that  rises  with  accumulation  gate  voltage  |VGB|.  Figure 
5.8  gives  examples  of  three  MOST's  from  two  thin  gate  oxide  technologies:  a  4.5nm 
technology  in  part  (a)  and  a  3.0nm  technology  in  part  (b).  To  focus  on  the  accumulation 
VGB  range,  in  part  (b),  we  only  show  nMOST  data  at  VGB<0  and  pMOST  data  at 
VGB>0.  The  Families  of  DCIV  curves  correspond  to  top-emitter  forward  biases  from 
50mV  to  350mV  on  a  W/L=  20^m/0.25nm  pMOST  in  part  (a)  and  drain-emitter 
forward  biases  from  200m V  to  400m V  on  W/L=10|im/0.32nm  MOSTs  in  part  (b).  All 
IB-VGB  curves  show  peaked  structures  on  top  of  the  baselines  that  rise  with 
accumulation  |VGB|. 
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Fig.  5.8  Experimental  DCIV  data  with  rising  IB  in  the  accumulation  gate  voltage 
range,  (a)  TE-DCIV  data  from  a  4.5nm  oxide  technology, 
W/L=20u70.25u.m,  pMOST.  The  IB-VGB  curves  correspond  to  TE  forward 
bias  of  lOOmV  to  350mV.  (b)  DE-DCIV  data  from  a  3.0nm  oxide 
technology,  W/L=10u.m/0.32nm  for  both  nMOST  (left)  and  pMOST 
(right).  DE  forward  bias  varies  from  200mV  to  400mV  200mV  to  400mV 
with  a  step  of  50mV.  The  arrows  point  to  IB  peaks  in  the  accumulation  gate 
voltage  range. 
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Fig.  5.9  Forward  bias  VPN  dependence  of  the  difference  base-terminal  current  AIB 
for  the  20u.m/0.25nm  pMOST  in  Fig.  5.8(a).  AIB=  |IB(2V)-IB(0V)|,  see 
inset.  The  dash  line  is  a  least-squares-fit  to  AIB  =  I0  exp(qVPN/nkT)  which 
gives  n=2.11±0.02. 


151 

The  rising  baselines  are  mostly  from  currents  flowing  through  the  emitter/base 
p/n  junctions,  which  were  verified  by  comparing  the  emitter  current  with  the  base 
current.  The  following  features  of  the  baselines  are  not  consistent  with  the  theories  of 
surface  p/n  junction  tunneling:  (1)  Experimental  IB.acc  increases  with  VPN  as  well  as 
|VGB|  while  theoretical  JPN  decreases  with  VPN  (see  Fig.  5.2);  (2)  The  rising  baselines 
have  low  threshold  |VGB|  around  IV  which  is  almost  independent  of  VPN  while  the  JPN 
theory  predicts  a  significant  increase  of  the  threshold  at  higher  VPN.  In  particular,  we 
analyzed  the  VPN  dependence  of  IB  at  VGB=2V  for  the  20|im/0.25|im  pMOST  in  Fig. 
5.9.  The  least-squares-fit  (dash  line)  of  AIB  =  |IB(2V)  -  IB(0V)|  to  exp(qVPN/nkT) 
determines  an  ideality  factor  of  n=2.1 1±0.02. 

The  theory-experimental  discrepancy  can  be  understood  from  the  finite  doping 
concentrations  in  drain  and  source  regions.  The  ideal  JPN  theory  neglects  surface  band 
bendings  in  the  gate-overlapped  emitter  region,  i.e.,  VSE=0,  which  is  valid  in  the  limit  of 
infinite  doping  concentrations.  In  real  MOSTs  as  shown  by  Fig.  2.5,  under  channel 
accumulation  gate  voltages,  the  energy  bands  in  the  overlap  region  bend  towards  carrier 
depletion  and  even  inversion  conditions,  in  the  same  direction  of  the  energy  bend 
bending  in  the  channel  region,  thus  decreasing  the  potential  barrier  across  emitter/base 
p/n  junction  at  the  surface  which  prevents  the  direct  p/n  junction  tunneling.  In  Fig. 
5.10(a),  we  compare  the  x-direction  surface  band  bending  in  the  overlap  region  (VSE) 
for  five  n+  drain  doping  concentrations  (NDD)  from  1018cm"3  to  1020cm"3.  The  dash  line 
gives  the  x-direction  surface  band  bending  in  the  p-substrate  (Vsx)  with  a  doping 
concentration  of  1017cm-3.  The  potential  barrier  height  of  the  p/n  junction  at  the  surface 
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Fig.  5.10  Effects  of  finite  n+  drain  doping  concentrations  on  (a)  surface  band  bending 
in  gate-overlapped  drain  region  and  (b)  potential  barrier  height  at  the 
surface  (X=0)  of  the  n+drain/p-base  junction.  The  five  solid  lines 
correspond  to  five  n+  drain  doping  concentrations.  The  dash  line  in  part  (a) 
is  the  surface  band  bending  in  the  p-base  and  another  dash  line  in  part  (b) 
shows  the  potential  barrier  for  infinitely  large  doping  concentrations  in  n+ 
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drain.  Parameters  used  include:  NAA=10   cm    in  p-base,  5x10    cm°  in  n+ 
gate  and  XOx=3.0nm. 
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is  then  computed  from  eq.  (2.18)  and  shown  in  Fig.  5.10(b).  As  the  drain  doping 
concentration  increases  from  1018crrT3  to  1020cm~3,  the  maximum  barrier  height 
increases  from  0.3eV  to  about  0.8eV  under  a  junction  forward  bias  of  0.2V.  In  all  cases, 
the  barrier  height  drops  at  large  |VGB|  due  to  carrier  depletion  at  the  drain  surface,  in 
contrast  to  the  ideal  case  (dash  line)  with  VSE=0  that  keeps  rising  at  large  |VGB|.  Fig. 
5.10(b)  clearly  shows  that  the  potential  barrier  height  at  the  surface  of  the  drain/base 
junction  does  not  exceed  E  =1.12eV  which  prevents  direct  surface  p/n  junction 
tunneling  in  MOS  transistors.  In  general,  interband  tunneling  in  silicon  requires  a 
minimum  potential  barrier  of  E  =1.12eV  to  cross  the  conduction  band  and  valence  band. 
However,  if  there  are  gap  states  such  as  those  from  interface  traps,  the  minimum  barrier 
height  would  be  lowered  for  tunneling  through  the  gap  states  [58]. 

DCIV  curves  in  Figs.  5.8(a)  and  5.8(b)  all  show  peaked  structures  as  pointed  to 
by  arrows  on  top  of  IB.acc,  most  prominent  in  the  20(xm/0.25u,m  pMOST.  The  peaks  are 
around  -1.4V  for  the  nMOST  and  between  -1.0V  to  -1.3V  for  the  pMOSTs.  These 
peak  gate  voltages  and  their  shift  with  VPN  (positive  for  nMOST  and  negative  for 
pMOST)  are  characteristic  of  SRH  recombination  at  interface  traps  in  the  LDER 
overlapped  by  the  gate  (see  Figs.  3.3  and  3.7  in  chapter  2).  Smaller  peaks  not  labeled 
are  visible  in  Fig.  5.8(b)  at  the  weak  inversion  |VGB|  around  0.3V-0.4V,  which  are  from 
SRH  recombination  at  NIT's  located  over  the  MCR. 

The  existence  of  interface  traps  in  the  gate-overlapped  LDER  suggests  interface 
trap  assisted  tunneling  as  the  origin  for  the  rising  IB.acc,  which  will  be  analyzed  in  detail 
in  the  next  section.  In  chapter  3,  we  demonstrated  the  extraction  of  transistor  parameters 
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such  as  channel  doping  concentrations  from  the  DCIV  peak  around  the  weak  inversion 
VGB  range.  We  will  next  show  a  similar  application  to  determine  the  doping 
concentration  in  the  LDER  from  the  DCIV  peak  in  the  accumulation  VGB  range. 

Figure  5.1 1  compares  theoretical  SRH  recombination  rate  RSS1  (dash  lines)  in 
the  LDER  with  the  experimental  DCIV  peak  in  the  20^im/0.25u.m  pMOST  (solid  lines). 
Theoretical  SRH  recombination  rate  per  unit  interface  trap  was  obtained  using  a  mid- 
gap  interface  trap  energy  level  ETI=0  and  equal  electron  and  hole  capture  rate 
coefficients  cns=c  =10_8cm3s_1.  The  nominal  oxide  thickness  of  4.5nm  for  the 
technology  was  used  in  the  theory.  The  best  theory-experiment  match  in  peak  gate 
voltage  as  shown  in  Fig.  5.11,  was  achieved  with  NAA=2.0xlO,8cm-3  in  the  LDER  of 
the  pMOST. 

The  variation  of  peak  gate  voltage  with  LDER  impurity  doping  concentration  is 
shown  in  Fig.  5.12.  Under  a  p+drain/n-base  junction  forward  bias  of  lOOmV,  the  SRH 
recombination  rate  peak  (Rssi-pk)  shifts  from  0.9V  to  2.3V  as  the  LDER  doping 
concentration  increases  from  5xl017cm~3  to  1.5xl019cm-3.  This  demonstrates  the 
sensitivity  of  the  technique  which  is  better  than  10%  for  the  typical  mid-1018cm-3  range 
provided  that  the  experimental  peak  voltage  can  be  determined  within  50mV.  The  IB 
peak  in  the  pMOST  data  of  Fig.  5.8(a)  was  selected  for  the  fit  in  Fig.  5.1 1.  For  the  two 
MOST's  in  Fig.  5.8(b),  the  theory  in  Fig.  5.12  predicts  LDER  doping  concentrations  of 
2.5xl018cm-3  in  the  nMOST  from  VGB.pk=1.40V  and  7xl017cm"3  in  the  pMOST  from 
VGB-Pk=0-95V. 
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Fig.  5.1 1  Theory-experiment  comparison  of  IB  peak  in  the  accumulation  gate  voltage 
range.  Experimental  TE-DCIV  data  are  those  of  the  20|J.m/0.25u.m  pMOST 
presented  in  Fig.  5.8(a).  Experimental  curves  are  scaled  to  match  the 
theoretical  peak  height  in  the  low  VPN  range.  Theoretical  peaks  were 
computed  using  SRH  theory  for  recombination  at  interface  traps  located  in 
the  LDER,  with  the  following  parameters:  mid-gap  trap  level  ETI=0; 
XoX=4.5nm;  Nxx=5xl019cm~3  in  the  gate,  and  NAA=2.0xl018cm~3  in  the 
LDER.  The  LDER  concentration  is  determined  from  the  theory- 
experimental  match  of  peak  voltage  in  the  low  VPN  curves. 
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Fig.  5.12  Theoretical  SRH  recombination  rate  at  interface  traps  in  the  gate  overlapped 
drain  region  at  a  p+drain/n-base  junction  forward  bias  of  lOOmV,  with  the 
impurity  doping  concentrations  in  p+  drain/source  LDER  as  a  parameter. 
Locations  of  the  IB  peak  in  the  20u,m/0.25u.m  pMOST  is  indicated  by  an 
arrow.  NAA=5xl019cm-3  in  p+  gate  and  Xox=4.5nm. 
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The  experimental  IB  curves  in  Fig.  5.11  were  scaled  to  match  the  RSS1  peak  at 
VPN=100mV.  The  scaling  factor  provides  an  estimate  on  interface  trap  density  NIT  in 
the  overlap  region.  The  number  of  interface  traps  in  the  overlap  LDER  area  of  the 
pMOST  is  given  by 

(WxLoverlap)NIT    =    AlB.pk/qRssl_pk  (5.17) 

=    13fA    /  (1.6xlO~19Cx325s1)     =    2.5xl02 
where  Loverl    =0.1u,m  is  the  length  of  overlap  region,  W=20fj.m  is  the  channel  width, 
13fA  is  the  peak  current  at  VPN  =  lOOmV.  The  result  of  (5.17)  normalized  by  the 
overlap  area  gives  N1T=1.3xl010cm~2  in  the  LDER,  which  is  higher  than  the  typical  2- 
8xl09cm-2  measured  in  the  MCR  on  prestress  MOSTs  [5,33]. 

Though  theory-experimental  match  in  Fig.  5.11  is  almost  perfect  for  the  peak 
voltage  and  peak  current  in  the  low  VPN  range  from  lOOmV  to  200mV,  deviations  exist 
at  higher  VPN's.  IB  peaks  at  300mV  and  350mV  appear  much  higher  than  their 
theoretical  predictions.  The  agreement  will  be  significantly  improved  if  the  VGB- 
independent  baseline  is  subtracted  from  each  of  these  two  experimental  curves.  The 
baseline  has  the  exp(qVPN/kT)  dependence  while  the  LDER  peak  has  the 
exp(qVPN/2kT)  from  the  ideal  SRH  theory,  thus  the  baseline  becomes  less  significant  at 
lower  VpN's  to  give  better  theory-experimental  match.  The  experimental  peak  voltage 
shift  lOOmV  more  in  the  negative  VGB  direction  than  the  theoretical  prediction  as  VPN  is 
increased  from  lOOmV  to  350mV.  We  ascribe  this  discrepancy  to  two-dimensional 
effects  such  as  the  fringe  field  in  the  LDER  not  considered  by  the  theory. 
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Similar  experimental  characteristics  on  thicker  oxide  nMOST  were  reported  by 
Acovic  and  co-workers  [117]  who  measured  forward-biased  drain  current  as  a  function 
of  gate  voltage  at  77K.  The  overall  increase  of  drain  current  with  negative  gate  voltage 
was  attributed  to  interface  trap  assisted  tunneling.  Their  sample  has  a  thick  oxide  of 
38nm,  and  there  is  a  primary  peak  in  the  drain  current  that  shifts  from  VGB=-14V  at 
VpN=0.575V  to  VGB=-12V  at  VpN=0.975.  The  direction  of  the  shift  is  consistent  with 
interface  traps  in  the  gate  overlapped  LDER.  Using  the  SRH  recombination  model  just 
described,  we  found  that  a  drain  doping  concentration  of  2.2xl019cm-3  would  perfectly 
match  the  peak  location,  which  is  not  inconsistent  with  the  high  source/drain 
implantation  dose  of  1016cm-2  (120keV,  As)  used  in  their  sample. 

5.5  Interface-Trap  Assisted  Tunneling 

The  rising  IB.acc  can  not  be  explained  by  direct  surface  p/n  junction  tunneling 
due  to  small  potential  barrier  at  the  p/n  junction  surface,  while  the  IB.acc  peaks  provide 
evidence  of  interface  traps  in  the  gate  overlapped  lowly-doped  (drain)  extension  region 
(LDER).  In  this  section,  we  will  consider  interface  trap  assisted  tunneling  (TAT)  in  the 
LDER.  To  enable  efficient  computation  using  one-dimensional  analytical  formulas,  we 
restrict  the  tunneling  in  the  x-direction  perpendicular  to  the  Si/Si02  interface.  In 
general,  tunneling  would  occur  along  the  direction  of  electric  field  at  the  surface,  which 
is  not  limited  to  the  x-direction  near  the  junction  space  charge  region  where  the  field 
distribution  is  two-dimensional. 
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Fig.  5.13  Cross-sectional  view  and  energy  band  diagram  illustrating  interface-trap 
assisted  tunneling  along  the  x-direction  at  the  gate-overlapped  n-i-  drain 
region.  At  large  negative  gate  voltages,  there  is  an  electron  depletion  layer 
(thickness  xd  )  and  even  a  hole  inversion  layer  at  the  surface  of  n-i-  drain, 
while  there  is  a  hole  accumulation  layer  at  the  surface  of  p-base.  The 
energy  band  diagram  shows  three  processes:  (1)  Conduction  band  electron 
at  x=XT  tunnel  into  interface  trap  (with  energy  ET  below  Ec)  at  x=0,  (1') 
Thermal  capture  of  a  conduction  band  electron  by  the  trap  at  x=0,  and  (2) 
Thermal  capture  of  a  valence  band  hole  by  the  trap  at  x=0.  Fp,  FN  are 
electron  and  hole  quasi  Fermi  levels,  VPN  is  the  forward  junction  bias  and 
VSE  denotes  the  surface  band  bending  in  n+Drain. 
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The  trap-assisted  tunneling  pathway  is  illustrated  in  Fig.  5.13  with  a  cross- 
sectional  view  of  the  gate-overlapped  n+D/p-Base  junction  and  an  energy  band  diagram 
along  the  x-direction.  Under  a  large  negative  gate  bias  VGB  and  a  forward  junction  bias 
VPN,  n+D  surface  is  electron  depleted  (n-dep)  and  has  a  hole  inversion  layer  (p-inv). 
The  hole  inversion  layer  is  connected  to  the  hole  accumulation  layer  (p-acc)  at  the 
surface  of  p-Base.  The  triangle  symbolizes  an  interface  trap  in  the  n+D  (LDER)  with  a 
discrete  energy  level  Ej  below  the  conduction  band  edge  Ec.  The  interface  trap  can 
either  thermally  capture  a  conduction  band  (CB)  electron  at  the  surface  x=0  (process  1') 
or  elastically  capture  a  CB  electron  at  x=XT  through  quantum  mechanical  tunneling 
(process  1).  Energy  conservation  requires  that  band  bending  at  XT  equals  VSE  -  ET/q, 
where  VSE  is  the  total  surface  band  bending  in  the  x-direction.  The  relative  probabilities 
of  these  two  processes  depend  on  the  gate  voltage.  Process  2  is  the  thermal  capture  of  a 
valence  band  hole  at  X=0,  which  is  supplied  by  the  base-terminal  through  the  p-Base 
hole  accumulation  layer  and  the  n+Drain  hole  inversion  layer.  The  combination  of 
processes  1'  and  2  is  the  usual  SRH  recombination  process.  At  more  negative  VGB, 
electron  concentration  at  x=0  drops  exponentially  which  decreases  SRH  recombination 
rate,  while  the  probability  of  tunneling  (process  1)  increases  exponentially  due  to  a 
thinner  tunneling  barrier.  The  combination  of  processes  1  and  2  is  the  interface-trap 
assisted  tunneling  process,  which  should  dominate  at  large  negative  VGB. 

We  will  next  develop  a  theory  for  surface  p/n  junction  current  that  takes  into 
account  of  all  three  processes  (1,1 'and  2).  We  start  with  the  steady-state  condition: 
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dnIT/dt    =    0    =    dnIT/8t(SRH)    +    dnIT/dt  (tunnel) 

=      (CnsnsPlT-ensnlT)      "      (CpsPsniT-ensniT)      +    JT^  (5.18) 

nIT  and  pIT  are  the  electron-occupied  and  electron-depleted  interface  trap  areal  densities 
respectively,  and  nIT  +  Pu  =  NITT  is  the  total  density  of  interface  traps.  ns  and  ps  are  the 
surface  electron  and  hole  densities,  and  cns,  c  ,  ens,  e  are  electron  and  hole  capture  and 
emission  rate  coefficients  respectively  at  the  surface.  The  last  term  on  the  R.H.S.  of 
(5.18),  comes  from  electron  tunneling  from  the  conduction  band  at  x=XT  to  interface 
traps  at  x=0  and  JT  denotes  the  tunneling  current  density. 

We  then  derive  a  formula  for  JT  based  on  (5.2)  and  the  DOS  of  the  final  state 
(interface  trap)  in  the  tunneling  process,  which  is  D2(E)=  pIT  8(E+qVSE+ET),  where  E  is 
the  electron  energy  of  the  initial  conduction  band  state  at  X=XT.  The  delta  function 
reflects  the  conservation  of  electron  energy.  Since  the  tunneling  probability  is 
proportional  to  the  density  of  final  states,  we  apply  a  dimensionless  multiplication  factor 
to  the  integrand  of  (5.2)  to  reflect  the  specific  DOS  for  tunneling  into  interface  states. 
The  result  is: 


Jm         — 


47T/jqimym2 


h3 


-qVSE  PIT    6(E    +qVSE+ET) 

T(E)  •  [f  (E)-f  (E  +  qVPN)]E.  dE     (5.19) 

0  m0/(7Th2) 


=    q   pTT   Tn  (5.20) 


T0    -    /?1lmymz    /(m0h)     •    {T  (E)  [  f  (E) -f  (E+qVpN)  ]  E}  |  E=-qVSE-ET  (5.21) 

T0  in  (5.21)  is  the  tunneling  rate  per  second,  which  is  evaluated  at  E=-qVSE-ET. 
It  is  mainly  determined  by  two  factors:  (1)  T(E),  the  barrier  transmission  probability, 
which  increases  with  accumulation  |VGB|  due  to  barrier  thinning  and  (2)  f(E),  the 
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electron  occupation  factor  at  x=XT  which  decreases  with  |VGB|.  To  evaluate  T(E),  we 
use  the  average-field  non-parabolic  formula  of  (5.9)  and  replace  the  barrier  height  E  by 
Ej.  The  average  field  in  the  tunneling  region  is:  F  =  [F(x=0)+F(x=XT)]/2.  The 
electric  fields  can  be  evaluated  using  the  standard  one-dimensional  MOS  equations. 

The  steady-state  electron-occupied  interface  trap  density  can  be  derived  from 
(5.18)  and  (5.20)  in  terms  of  T0,  which  gives: 


Nit  cnsNs   +    eps    +   To 


(5.22) 


Njtt  CnsNs    +    CpsPs     +     ens     +     eps     +    T0 


Equation  (5.22)  reduces  to  the  familiar  SRH  formula  when  the  tunneling  rate  is  much 
lower  than  the  thermal  capture  rates.  In  the  other  limit  when  the  tunneling  rate 
dominates,  the  trap  (electron)  occupation  factor  approaches  unity. 

Combining  (5.18)  and  (5.22),  we  derive  the  following  formula  that  describes 
steady-state  current  follow  due  to  all  three  processes  (1,1'  and  2)  at  the  surface  of  the 
LDER: 

<cnsNs    +    To)     cpsPs    "     enseps 
JLDER     =     3    NITT    ~     QN^rpR-L 

cnsNs   +   cpsPs    +    ens    +    eps   +   T0  (5.23) 

Rj  is  defined  as  a  generalized  recombination  rate  per  interface  trap,  which  reduces  to  the 
SRH  rate  RSS1  if  the  tunneling  process  is  not  considered.  When  the  tunneling  rate  T0 
(process  1)  is  much  higher  than  the  thermal  electron  capture  rate  cnsNs  (process  1'), 
(5.23)  is  simplified  to: 

J-0Cps    s     ~     enseps 
"LDER     =     "5    NITT 

cpsPs   +   ens    +    eps    +   T0  (5.24) 


163 

Equation  (5.24)  is  the  formula  for  trap-assisted  tunneling  (processes  1  and  2).  For  near 
mid-gap  interface  trap  energy  levels,  the  following  equalities  hold:  ens=cnsnj  and 
eps~cpsni-  Therefore  ens  and  e  are  both  negligible  compare  to  c  Ps  and  T0.  In  the  limit 
of  strong  accumulation  such  that  c   Ps  »T0,  (5.24)  can  be  further  simplified  to: 

TOcpsPs 
^lder   =    *3  NITT  ~   qNITTT0 

cpsPs    +   T0  (5.25) 

which  describes  the  TAT-limited  LDER  current. 

We  will  next  examine  the  effects  of  gate  bias,  forward  junction  bias,  n+D  doping 
concentration,  gate  oxide  thickness,  interface  trap  energy  level  and  temperature  on  the 
LDER  current  using  the  most  general  formula  (5.23). 

The  gate  voltage  dependence  of  normalized  LDER  current  (Ri=Jlder^^itt)  ^s 
shown  in  Fig.  5.14,  along  with  that  of  tunneling  rate  T0  [eqn  (5.21)],  the  electron 
thermal  capture  rate  cnsNs  and  the  hole  thermal  capture  rate  c  Ps  for  a  n+D/p-sub 
junction  forward  biased  at  200mV.  Parameters  used  in  the  computation  include 
Xox=4.5nm,  NDD(Drain)=1018cm"3  and  NDD(gate)=5xl019cnT3,  ET=Eg/2  and  T=300K. 
The  gate  voltage  range  (-3.0V  to  OV)  can  be  divided  into  three  regions:  (1)  Tunneling 
current  limited  region  at  VGB<-1.6V,  where  c  Ps  >T0>cnsNs  and  Rj  =  T0;  (2)  SRH 
recombination  current  limited  region  at  VGB>-1.45V,  where  c  Ps,  cnsNs>T0=0,  and 
Ri=R$sp  (3)  A  narrow  transition  region  around  VGB=  -1.5V  in  which  T0  and  cnsNs  are 
compatible.  The  peak  at  VGB=-1.1V  is  from  SRH  recombination  at  interface  traps  in 
the  LDER,  which  coincide  with  the  point  where  cpsPs  =  cnsNs.  In  the  second  region,  T0 
drops  to  zero  since  the  surface  band  bending  is  not  large  enough  to  provide  a  direct 
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Fig.  5.14  Gate  voltage  dependence  of  the  generalized  recombination  rate  Rj  in  the 
gate-overlapped  drain  region  under  a  forward  bias  of  200mV  on  the 
n+drain/p-base  junction.  Three  rates  associated  with  interface  traps  are  also 
shown  which  include  electron  and  hole  thermal  capture  rates  cnsNs  and 
c    Ps,  and  tunneling  rate  T0.  Parameters  used  are:  NDD=1018cm-3  in  n+ 


,19, 


,-3; 


drain,  NDD=5xl0lvcm~J  in  n+  gate,  Xox=4.5nm,  E^O  and  T=298K. 
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band-to-trap  tunneling  pathway.  In  the  first  region  where  the  LDER  current  is  limited 
by  tunneling  rate,  the  LDER  current  increases  with  stronger  accumulation  gate  voltage, 
which  is  the  characteristic  of  IB.acc  data  from  DCIV  measurement. 

Figure  5.15  consists  of  four  parts,  each  part  shows  a  family  of  RrVGB  curves 
with  one  of  the  four  variables  as  a  parameter:  (a)  the  forward  junction  bias  VPN;  (b) 
doping  concentration  in  n+Drain  NDD;  (c)  gate  oxide  thickness  Xox  ar»d  (d)  interface 
trap  energy  level  ETI.  The  solid  line  in  each  part  is  the  same  curve  analyzed  in  Fig.  14 
which  serves  as  a  reference.  To  enable  lineshape  comparison,  the  dashed  curves  in  Figs. 
5.15(b)  and  5.15(c)  are  shifted  along  the  VGB  axis  to  yield  the  same  SRH  peak  at  -1.1  V. 
The  VPN  dependence  in  Fig.  5.15(a)  is  in  qualitative  agreement  with  the  experimental 
DCIV  data  in  the  strong  accumulation  range  (Fig.  5.8):  the  TAT  current  increases  with 
VPN  at  a  fixed  VGB.  However,  the  LDER  current  increases  faster  with  VPN  than 
experimental  IB.acc  data.  In  terms  of  ideality  factor,  the  TAT  theory  gives  n=1.27 
VGB=-3V,  while  the  DCIV  data  gives  n=2  (Fig.  5.9).  Figure  5.15(b)  shows  that  the 
tunneling-limited  LDER  current  at  VGB<-1.5V  increases  at  higher  n+D  doping 
concentrations.  The  peaked  structure  near  the  onset  of  tunneling  region  is  due  to  the 
interplay  between  (1)  the  number  of  electrons  available  for  tunneling  at  x=XT,  which 
decreases  with  |VGB|  and  (2)  the  barrier  transmission  probability  T(E),  which  increases 
with  |VGB|.  The  peak  widening  at  higher  NDD  concentrations  is  due  to  larger  |VGB| 
required  to  invert  the  n+D  surface.  Fig.  5.15(c)  shows  that  thin  gate  oxide  leads  to  high 
LDER  current  in  the  tunneling-limited  range,  which  is  consistent  with  DCIV 
experiments  that  the  rising  IB  in  accumulation  range  is  more  significant  in  thin  oxide 
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Fig.  5.15  Gate  voltage  dependence  of  the  generalized  recombination  rate  Rj  at  (a) 
four  junction  forward  biases,  (b)  four  n+drain  impurity  doping 
concentrations,  (c)  four  gate  oxide  thicknesses,  and  (d)  three  interface  trap 
energy  levels.  The  solid  line  is  the  reference  curve  with  VPN=200mV,  NDD 
=  1018cm~3  Xox=4.5nm  and  ETI=0.  In  each  part,  only  one  parameter  is 
varied  from  the  reference  value.  In  part  a),  the  ideality  factor  n=1.27  at 
VGB=-3.0V  and  n=2.0  for  the  peaks.  In  part  b)  and  c),  the  dashed  lines  are 
shifted  to  match  the  peak  voltage  of  the  solid  line  at  -1.10V. 
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MOSTs.  The  effect  of  interface  trap  energy  level  is  illustrated  in  Fig.  5.15(d).  A  mid- 
gap  trap  has  ETI=0.  For  trap  levels  closer  to  the  silicon  conduction  band  edge,  the 
threshold  for  tunneling  is  lowered  due  to  a  lower  tunneling  barrier  (Ej  =  E  J2  -  Eyj),  a 
peaked  structure  appears  after  the  onset  of  tunneling,  and  the  tunneling  rate  rises  slower 
with  |VGB|  in  the  strong  accumulation  range. 

Theoretical  LDER  currents  at  three  different  temperatures  are  shown  in  Fig. 
5.16.  The  same  parameters  as  those  in  5.14  were  used  in  the  computation.  Figure 
5.16(a)  shows  unsealed  Ri-VGB  curves  which  exhibits  a  significant  temperature 
dependence  even  in  the  tunneling-limited  range,  in  contrast  to  the  notion  that  tunneling 
current  has  only  a  weak  temperature  dependence.  This  is  because  the  tunneling 
probability  is  exponentially  dependent  on  the  tunneling  distance  XT,  which  is  not  fixed 
in  TAT.  XT  has  a  weak  temperature  dependence  through  the  MOS  capacitor  equation. 
In  order  to  compare  the  relative  strength  of  temperature  dependence  in  SRH  range  and 
tunneling  range,  we  scaled  down  the  80°C  curve  and  scale  up  the  -30°C  curve  in  Fig. 
5.16(b)  to  yield  the  same  SRH  peak  as  the  room  temperature  curve.  In  the  tunneling 
range  on  the  negative  VGB  side  of  the  SRH  peak,  the  -30°C  curve  is  above  the  80°C 
curve,  which  means  the  tunneling  rate  has  a  'weaker'  temperature  dependence  than  that 
of  the  SRH  recombination  peak.  On  the  other  side  of  the  SRH  peak,  the  -30°C  curve  is 
below  the  80°C  curve,  since  in  this  range  R{  =  c  SPS  °c  n?  exp(qVPN/kT)  has  a  stronger 
temperature  dependence  than  Rj  <*  rij  exp(qPN/2kT)  at  the  peak. 

In  conclusion,  we  have  developed  a  simple  trap-assisted  tunneling  (TAT)  model 
to  account  for  majority  carrier  tunneling  into  interface  traps  located  in  gate-overlapped 
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Fig.  5.16  Gate  voltage  dependence  of  the  generalized  recombination  rate  in  the 
LDER  at  three  temperatures:  -30°C,  25°C  and  80°C.  (a)  The  three 
unsealed  R,  curves,  (b)  Two  scaled  scaled  R,  curves  along  with  the 
unsealed  room  temperature  curve.  The  scaling  in  part  (b)  matches  the  peak 
levels  at  -1.10V.  Parameters  other  than  temperature  are  the  same  as  those 
in  Fig.  5.14. 
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drain  or  source  regions  (LDER)  of  a  MOS  transistor.  The  model  predicts  a  lower 
tunneling  threshold  |VGB|  than  that  of  the  drain  junction  perimeter  direct  tunneling  (DT) 
theory.  The  TAT  rate  increases  with  both  MOST  accumulation  gate  voltage  and  drain 
junction  forward  bias,  which  are  the  essential  features  of  the  DCIV  base-terminal  current 
(IB)  in  accumulation  VGB  range.  However,  the  VPN  dependence  of  TAT  current  in  the 
LDER  is  not  in  quantitative  agreement  with  DCIV  data,  and  some  DCIV  data  such  as 
those  in  Fig.  3.10  show  no  apparent  IB  peaks  in  the  accumulation  VGB  range.  Two 
sources  may  contribute  to  this  discrepancy:  (1)  the  simple  tunneling  theory  does  not 
consider  the  detailed  interaction  potential  between  the  interface  trap  and  the  tunneling 
electron,  and  (2)  the  tunneling  direction  is  not  strictly  perpendicular  to  the  interface  and 
therefore  the  purely  one-dimensional  theory  fails  to  give  the  quantitative  correlation. 

We  have  also  considered  in  detail  another  possible  mechanism  not  related  to 
interface  traps  which  is  the  injection  current  in  the  LDER  region.  In  this  scenario,  the 
VGB  dependence  of  the  injection  current  comes  from  the  gate  voltage  modulation  of  the 
series  resistance  in  the  surface  (weak)  inversion  layer  of  the  LDER.  The  theory  of 
carrier  injection  over  the  inverted  channel  area  of  a  gated  p/n  junction  was  first 
considered  by  Sah  [3].  However,  the  following  characteristics  of  the  LDER  injection 
current  fails  to  correlate  with  experimental  DCIV  data:  (1)  The  injection  current  would 
quickly  saturate  at  strong  accumulation  gate  voltages  when  the  LDER  surface  is  strongly 
inverted  while  experimental  IB.acc  data  has  no  indication  of  saturation;  (2)  The  LDER 
injection  current  is  independent  of  interface  trap  density  while  the  experimental  IB.acc 
increases  after  channel-hot-carrier  stress  indicating  a  NIT  related  mechanism. 
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5.6  Summary 

Three  original  theories  of  surface  controlled  interband  tunneling  near  the 
drain/substrate  junction  perimeter  in  MOS  transistors  under  forward  junction  bias  have 
been  developed:  (1)  electron  tunneling  from  n+drain  into  classical  three-dimensional 
hole  states  in  p-substrate  surface,  (2)  electron  tunneling  from  n+drain  into  quantized 
two-dimensional  states  of  surface  hole  accumulation  layers,  and  (3)  interface  trap 
assisted  tunneling  (TAT)  in  the  gate-overlapped  drain  extension  region.  The  first  two 
theories  predict  a  high  threshold  gate  voltage  which  increases  as  VPN  increases,  not 
consistent  with  DCIV  data.  The  third  theory  has  lower  threshold  gate  voltage  and  are  in 
qualitative  agreement  with  the  experimental  DCIV  data  in  the  MOST's  accumulation 
gate  voltage  range  which  show  IB  peaks  on  a  rising  baseline.  At  a  given  VGB,  however, 
the  experimental  IB  data  (n=2)  rises  with  VPN  slower  than  that  from  the  TAT  theory 
(n=1.27  at  VGB=-3V).  The  quantitative  correlation  may  be  improved  by  taking  into 
account  of  lateral  voltage  drop  in  the  (weakly)  inverted  LDER  surface,  which  was 
shown  by  Sah  [3]  to  give  a  factor  of  two  increase  in  the  ideality  factor  in  the  limit  of  a 
long  surface  inversion  channel.  The  DCIV  peaks  in  the  accumulation  gate  voltage  were 
correlated  to  SRH  recombination  at  the  interface  traps  located  over  the  LDER  region. 


CHAPTER  6 
SUMMARIES  AND  CONCLUSIONS 

The  success  of  today's  semiconductor  industry  can  be  partially  attributed  to  the 
passivation  of  intrinsic  defects  at  the  interface  between  the  silicon  channel  and  the 
thermal  gate  oxide  to  give  extremely-low  electrically-active  interface  trap  density  on  as- 
manufactured  MOS  transistors.  The  direct-current  current  voltage  (DCIV)  method  has 
the  sensitivity  to  detect  these  low-density  interface  traps  by  measuring  the  gate  bias 
(VGB)  and  the  forward  junction  bias  (VPN)  modulated  electron-hole  recombination 
currents  from  the  base-terminal  (IB).  This  unique  sensitivity  of  DCIV  makes  it  a 
powerful  tool  for  monitoring  the  transistor  reliability  and  for  diagnosis  of  transistor 
design. 

In  previous  chapters,  fundamental  theories  which  relate  DCIV  characteristics  to 
device  and  material  parameters  were  presented.  It  is  shown  that  the  DCIV  peak  current, 
peak  gate  voltage  and  lineshape  can  each  be  used  to  extract  a  group  of  parameters: 
interface  trap  density  and  energy  levels  from  the  VPN  dependence  of  DCIV  peak 
currents;  spatial  profile  of  interface  traps  from  DCIV  lineshape  and  channel  length 
dependence  of  DCIV  peak  current;  and  impurity  doping  concentration  and  gate  oxide 
thickness  from  the  VPN  dependence  of  DCIV  peak  gate  voltages.  These  applications 
were  demonstrated  on  both  the  unstressed  MOS  transistors  as  well  as  hot-carrier- 
stressed  MOS  transistors.  Diffusion  limits  the  minority  carriers  injected  from  the 
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bottom-emitter  to  reach  at  the  Si/Si02  interface  at  low  bottom  emitter  biases  which 
produces  a  wider  and  lower  peak  current  than  the  theoretical  prediction.  The  interface 
trap  energy  levels  are  discrete  rather  than  continuous,  consistent  with  those  of  silicon 
dangling  bonds.  Spatial  density  of  interface  trap  over  the  channel  area  of  unstressed 
MOS  transistors  is  shown  to  rise  sharply  towards  the  source  and  drain  junctions  with  a 
characteristic  length  of  80nm.  Increased  interface  trap  density  was  observed  in  the  drain 
junction  space  charge  region  after  channel-hot-carrier  stress.  Transistor  design 
parameters  (doping  concentration  at  the  interface  and  gate  oxide  thickness)  extracted 
from  DCIV  characteristic  are  in  good  agreement  with  designed  values. 

Transistor  characterization  becomes  more  difficult  on  ultrathin  oxide  MOS 
transistors  in  current  and  future  technology  generations  due  to  the  interference  of  gate 
oxide  tunneling  current  which  appears  as  a  rising  baseline  in  DCIV  base-terminal 
current.  Prompted  by  the  symmetry  of  the  gate  tunneling  current-voltage  (Ig-Vgb) 
characteristics  in  p-channel  and  n-channel  MOS  transistors,  we  developed  a  first- 
principle  gate  oxide  tunneling  theory  that  shows  hole  tunneling  is  dominant  in  p+gate  p- 
channel  MOS  transistors  and  that  the  small  gate/drain  overlap  area  is  not  negligible 
compare  to  the  channel  area  which  gives  the  knee  point  observed  in  the  Iq-^gb  curves 
of  both  nMOST  and  pMOST.  The  theory  improves  the  traditional  WKB  barrier 
transmission  probability  by  taking  into  account  of  quantum  mechanical  wavefunction 
matching  at  the  abrupt  potential  boundaries.  This  improvemnet  is  critical  in  accounting 
for  the  experimental  data  at  low  electric  fields.  We  demonstrated  the  separation  of  two 
geometric  tunneling  pathways  in  gate  currents  and  the  extraction  of  impurity  doping 
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concentrations  in  the  channel  and  drain/source  areas  from  the  theory-experimental 
correlation. 

In  the  gate  accumulation  voltage  range  of  thin  oxide  MOS  transistors,  the  current 
flowing  through  the  drain/base  p/n  junction  increases  with  gate  voltage,  which  is  another 
source  of  rising  baseline  in  DCIV  curves.  We  developed  two  fundamental  theories  of 
interband  tunneling  to  compare  with  the  experimental  data:  (1)  direct  tunneling  at  the 
drain/base  junction  perimeter  and  (2)  interface  trap  assisted  tunneling  in  the  gate- 
overlapped  drain  region.  Energy  quantization  in  the  surface  accumulation  layer  of  the  p- 
base  is  shown  to  give  breakpoints  (not  peaks)  in  the  direct  tunneling  current  at  the  gate 
voltages  when  one  of  the  one-dimensional  discrete  energy  levels  crosses  the  conduction 
band  edge  on  the  n+drain  side.  Direct  tunneling  theory  predicts  a  high  threshold  gate 
voltage  which  increases  as  VPN  increases,  not  consistent  with  DCIV  data.  Interface  trap 
assisted  tunneling  has  lower  threshold  gate  voltage  and  can  qualitatively  explain  the  VPN 
dependence  of  the  experimental  IB  baseline.  In  particular,  the  existence  of  interface 
traps  in  the  overlap  region  is  confirmed  by  experimental  IB  peaks  at  gate  voltages 
consistent  with  the  designed  drain  impurity  doping  concentrations.  Better  theory- 
experimental  agreement  entails  numerical  two-dimensional  field  analysis  near  the 
drain/base  junction  perimeter. 

Tunneling  current  will  eventually  limit  the  further  thinning  of  gate  oxide  below 
1.5nm.  Alternative  high  dielectric  constant  gate  materials  are  being  actively  sought. 
Both  the  DCIV  method  and  tunneling  current  monitor  investigated  in  this  thesis  research 
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will  continue  to  be  useful  for  evaluation  of  interface  reliability  as  well  as  for  diagnosis 
of  transistor  design  in  the  next  generations  of  semiconductor  technologies. 


APPENDIX 
DERIVATION  OF  GATE  OXIDE  TUNNELING  FORMULAE 

The  appendix  presents  a  derivation  of  the  gate  oxide  tunneling  formulae  used  in 

chapter  4,  following  Bardeen's  transition  probability  approach  [59]  and  Harrison's 

independent-particle  tunneling  model  [60].   An  alternative  approach  using  barrier 

transmission  coefficient  (TC)  would  give  the  same  result  if  wavefunction  matching  at 

the  abrupt  potential  boundaries  are  taken  care  of  when  deriving  the  transmission 

coefficient.  However,  the  traditional  TC  method  uses  a  simple  exponential  transmission 

coefficient  with  unity  pre-exponential  factor  which  neglected  the  wavefunction 

matching  [53]. 

Wavefunctions  of  the  Initial  and  Final  States 

Consider  a  trapezoidal  Si02  barrier  and  two  electron  states:  a  right-hand  state 
Wsl  and  a  left-hand  state  4/S2.  *FS]  and  ^9  are  the  quantum  mechanical  wavefunctions 
of  the  two  states.  Because  of  the  one-dimensional  potential  barrier  [V(x)],  the  transverse 
part  of  the  two  wavefunctions  are  trivial  and  can  be  both  represented  by  exp(ikyy+ikzz). 
The  same  wave  numbers  (kv,k7)  in  the  initial  and  final  states  are  a  result  of  transverse 
momentum  conservation.  In  the  following  analysis,  we  will  focus  on  the  normal  part  of 
the  wavefunctions  ^six  anc*  X*/S2X-  For  notation  compactness,  we  will  drop  the 
subscript  X,  so  *FSi,  kj  and  Ej  will  be  used  to  represent  the  x-component  of  the 
wavefunction,  wave  vector  and  kinetic  energy  respectively  in  the  initial  state  (i=l)  and 
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Fig.  A.l  Energy  band  diagram  of  a  trapezoidal  Si02  barrier  and  wavefunctions  of  the 
initial  (right-hand)  state  and  the  final  (left-hand)  state  of  an  elastic  tunneling 
transition  process,  kj  and  k2  (k2>k1)  are  the  real  wavenumbers  outside  of 
the  barrier  region  and  k(x)  is  the  imaginary  wavenumber  inside  the  barrier 
region  (0<x<Xox)- 
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the  final  state  (i=2).  VFS1  and  XFS2  are  sinusoidal  at  x>0  and  at  x<Xox  respectively, 
while  they  decay  exponentially  in  the  opposite  directions  inside  the  potential  barrier,  as 
shown  in  Fig.  A.l. 

From  Bardeen's  transition  probability  point  of  view,  tunneling  of  an  electron 
from  the  right-side  of  the  barrier  to  its  left-side  is  equivalent  to  the  transition  of  an 
electron  in  an  initial  state  Wsl  to  a  final  state  W^.  With  the  free-electron  approximation 
outside  the  barrier  region  and  the  WKB  approximation  inside  the  barrier  region,  the 
wavefunctions  would  be  of  the  form: 


*si    =   Aicos  (k^x+A^ 
B, 


K(X) 


exp 


K(x)dx 


*S2    =   A2cos[k2(x-Xox)+A2] 


ficui 


exp 


^ox 


K(x)dx 


(x<0) 


(0<x<Xox) 


(x>Xox) 


(0<x<Xox) 


[A. la] 


(A. lb) 


[A. 2a] 


(A. 2b) 


kj  and  k2  are  the  wavenumbers  outside  the  barrier  region  in  the  initial  and  final  states 
respectively,  and  k(x)  is  the  imaginary  wavenumber  inside  the  barrier  which  is  given  by 
eqn.  (4.6).  The  various  constants  Ai?  A:  and  Bi  (i=l,2)  can  be  determined  from  (1) 
normalization  of  wavefunctions  and  (2)  wavefunction  matching  at  the  two  abrupt 
potential  boundaries.  In  the  following,  we  will  show  the  procedures  to  determine  Bj  and 
B2,  which  are  the  two  constants  pertinent  to  the  evaluation  of  the  transition  matrix 
element. 
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Continuation  of  both  the  probability  function  and  electron  flux  requires  that  both 
x¥s  and  (l/m)3vFs/3x  be  continuous  across  the  abrupt  boundaries.  Applying  these 
boundary  conditions  at  x=0  for  H^,  gives: 


and 


¥sl(x=0-)    =   AjCOsA!    =   4/sl(x=0+)    =    b1/t|jc(0) 

(l/msl)aTsl/dx(x=0-)    =    -(k1/mslsinA1)A1    = 
(l/mox)aTsl/ax(x=0  +  )    =    -(U(0)/mox)B1 


(A. 3) 


(A. 4) 


Elimination  of  Aj  from  (A. 3)  and  (A. 4)  results  in  the  following  expression  for  the 
magnitude  of  Bj: 


Bil    =    |Ax|x 


k(0, 


msf    k(0) 


mox    kl 


-1/2 


(A. 5) 


Repeating  the  above  procedures  for  T^  at  x=Xox>  we  derive  a  similar  expression  for 
the  magnitude  of  B2: 


A2|x 


ms2    *(Xox^~ 


K(X0X)       mox    kl 


(A. 6) 


Aj  and  A2  can  be  determined  from  the  normalization  of  wavefunctions  in  the  initial  and 
final  states  respectively,  while  assuming  a  small  and  negligible  amount  of  wavefunction 


penetration  inside  the  barrier  region: 
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|*sl|2dx 
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*si I 2dx 
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\Al  \^cos2 (k1x+A1)    dx    =    |A1|2(Lsl/2) 


(A. 7) 
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which  gives 


|AX|     =    (2/L£    vl 


Jsi 


(A. 8a! 
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and  a  similar  normalization  of  ^2  gives 

|A2|     =     (2/Ls2)1/2  (A. 8b) 

where  LS1  and  LS2  are  the  extensions  of  silicon  structure  to  the  right  side  (x<0)  and  to 
left  side  (x>Xox)  of  the  Si02  barrier  respectively.  Thus  the  wavefunctions  of  the  initial 
and  final  states  are  determined  except  for  the  phase  shift  A{  which  does  not  affect  the 
tunneling  transition  probability. 

Transition  Matrix  Element 

Bardeen  [59]  showed  that  the  transition  matrix  element  in  the  tunneling 
probability  expression  of  (4.1)  is  proportional  to  that  of  the  current  density  operator  in 
the  tunneling  direction  (Jx),  which  is  to  be  evaluated  between  the  initial  state  *FS1  and 
the  final  state  WS2  in  the  barrier  region  0<x<Xox-  Using  the  results  of  (A. lb),  (A. 2b), 
(A. 5)  and  (A. 6),  the  matrix  element  is  evaluated  as  follows: 


M12    =    <^S1  |M|^S2> 


■in<vsl|jj'ps2> 


(-ill)  (iB/2mox) 
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where  4^*,  A*  and  B*  (i=l,2)  are  the  complex  conjugate  of  *FSi,  A{  and  Bj  respectively. 
In  deriving  (A. 9c),  we  assumed  that  the  x-dependence  of  4^  and  H^  inside  the  barrier 
region  is  dominated  by  the  exponential  term. 

The  density-of-state  (DOS)  factor  per  unit  energy  interval  on  the  two  sides  of  the 
Si02  barrier  can  be  derived  from  the  DOS  factor  in  the  k-space: 

p1    =    (LS1/7T)  .  (dE1/dkl)~1   =     (LS1/7T)  .  (n^/h2^)  (A. 10a) 

p2    =    (LS2/7T)  .  (dE2/8k2)  _1  =    (LS2/tt)  .  (ms2/h2k2)  (A. 10b) 

The  length  (LS1  and  LS2)  dependence  disappears  in  the  following  dimensionless 

combined  factor  of  tunneling  probability  which  can  be  obtained  from  (A. 8a),  (A. 8b), 

(A.9d),  (A.  10a)  and  (A.  10b)  with  some  simple  algebra: 


p1p2|M12|2   =    l/(27r)2   x    exp 


-2 


'•ox 
K  (x)  dx 
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■moxkl        mslK(0)    -,  r  moxk2 


ms2*r(Xox) 


Lms2lC(Xox) 


moxk2 


Lmsl>c(0)      m^kj 
=    1/  (27T)  2    x    exp(-2u)     x    T0(E1,E2)  (A. 11) 

T0  and  d  defined  in  (A.l  1)  are  the  same  as  those  given  by  (4.4)  and  (4.5)  respectively. 
T0  depends  on  electron  kinetic  energies  in  the  initial  state  Ej  and  in  the  final  state  E2  via 
the  relationship  ki=(2msjEi)1/2  (i=l,2).  The  requirement  of  energy  conservation  gives 
E2=E1-qVox,  where  Vox  is  the  voltage  drop  across  the  Si02  layer. 


Tunneling  Current  Formulae 


In  this  section,  we  will  give  the  details  of  deriving  the  tunneling  current  formulae 


(4.3),  (4.10)  and  (4.1 1).  We  start  from  (4.2)  and  convert  the  summation  over  transverse 
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momentum  (k   and  kz)  to  a  double  integration  using  a  DOS  factor  of  l/(27i)2  in  the  two- 
dimensional  k-space: 


4/rq 
J   =   I 
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(fi-f2: 


|M12|2    pxp2 


(A. 13) 


fj  and  f2  are  electron  occupation  factors  in  the  initial  and  final  states  respectively.  They 
are  described  by  the  Fermi -Dirac  distribution  function: 

fi    =    {exp[Eix+Eit+Ec-EF)  /kTj+1}"1  (i  =  l,2)  (A. 14) 

In  (A.  14),  the  electron  kinetic  energy  above  the  conduction  band  minimum  Ec  is 
decomposed  into  a  normal  part  Eix  and  a  transverse  part  Eit=Eiy+Eiz,  EF  is  the  Fermi- 
level  in  the  silicon  layer. 

The  double  integration  of  fj  over  the  transverse  momentum  can  be  evaluated 
analytically  since  fj  only  depends  on  the  magnitude  of  the  transverse  momentum  for  the 
simplest  spherical  E(k)  relationship.  For  simplicity,  we  drop  the  subscript  Y  in  the 
following  evaluation: 


d^ 


dky   f   =    27T 
_oo 


ktdkt 


2  7ttndkT 


exp [ (Ex+Et+Ec-EF)  /kT] +1 
1  d(Et/kT) 


0      exp[ (Ex+Et  +  Ec-EF) /kT] +1 
=    (27rmdkT/h2)  .loge{l  +  exp[  (EF-EC-EX) /kT]  }  (A.  15) 

Equation  (A.  15)  is  also  valid  for  ellipsoidal  E(k)  relationship  such  as  those  of  silicon 


conduction  band  valleys,  with  the  2-dimensional  DOS  effective  mass  md=(mxmy),/2. 
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The  Fermi-levels  in  the  initial  and  final  state  differ  by  qVGB,  where  VGB  is  the 
applied  gate-to-substrate  voltage  bias.  Thus,  the  [  ]  term  in  (A.  13)  gives  the  following 
results: 

l+exp[ (EF-EC-EX) /kT] 


2  7TmdkT 
[    ]  -   loge 


h2 


l+exp[ (EF-Ec-Ex-qVGB)  /kT] 


(A. 16) 


The  electron  tunneling  current  formula  (4.3)  is  derived  by  using  the  results  of 
(A.  11)  and  (A.  16)  in  (A.  13),  which  gives: 


4;rmdkT 
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exp[-2u(Ex)  ]T0(Ex)log€ 
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,       l+exp[  (EF-EC-EX) /kT] 
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ll+exp[  (EF-Ec-Ex-qVGB)  /kT] 

(A. 17) 


The  integration  in  (A.  17)  is  over  the  allowed  kinetic  energy  (x-component)  range  of  the 
initial  state  from  zero  to  Em=<j)Bn,  where  (})Bn  is  the  Si/Si02  conduction  band  offset. 

The  hole  tunneling  formula  (4.11)  can  be  derived  in  the  same  way  by 
considering  the  silicon  valence  band  hole  states  on  the  two  sides  of  the  Si02  barrier  as 
the  initial  and  final  states  and  taking  into  account  of  all  allowed  hole  kinetic  energies 
from  zero  to  Em=<j)B  ,  where  (j)B  is  the  Si/Si02  valence  band  offset. 

For  valence  electron  tunneling,  the  initial  state  is  a  silicon  valence  band  electron 
state  on  one  side  of  Si02  barrier  and  the  final  state  is  a  silicon  conduction  band  electron 
state  on  the  other  side  of  the  barrier.  Elastic  tunneling  path  opens  when  the  valence 
band  overlaps  the  conduction  band,  which  requires  a  minimum  oxide  voltage  drop  of 
qEg.  Under  such  oxide  voltages,  the  valence  band  states  are  fully  occupied  while  the 
conduction  band  states  on  the  other  side  of  Si02  are  nearly  empty,  which  gives  f j— f2  ~ 
1-0  =  1.   The  amount  of  conduction  band  and  valence  band  overlap  determines  the 
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allowed  kinetic  energy  range  for  elastic  tunneling  0<E1=Elt+Elx<Em=qVOx-Eg-  The 
allowed  range  for  the  transverse  part  of  kinetic  energy  is  thus  0<Elt<Em-Elx  which  sets 
an  upper  limit  for  the  transverse  momentum  kt<km=[2md(Em-Elx)/h2]1/2.  The  result  of 
(A.l  1)  for  the  transition  matrix  element  applies  to  valence  electron  tunneling  as  well, 
where  E^kj)  and  E2(k2)  are  referenced  to  the  valence  band  edge  and  conduction  band 
edge  respectively.  Starting  from  the  general  tunneling  current  formula  (4.2),  we  derive 
the  valence  electron  tunneling  formula  (4.10)  as  follows: 
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Triangular  Barrier  Tunneling 

For  a  triangular  barrier  illustrated  in  Fig.  A. 2,  the  right-hand  initial  state  is 
described  by  the  same  wavef unction  x¥sl  of  the  trapezoidal  barrier.  The  procedures  of 
matching  the  wavefunction  at  x=0  demonstrated  in  (A. 3)  to  (A. 5)  apply  to  the  present 
case  of  triangular  barrier  as  well.  However,  for  the  left-hand  final  state  VFS2,  there  is  no 
discontinuity  at  the  potential  boundary  x=X'.  We  may  write  ¥52  using  the  WKB 
approximation  at  x>X'  and  x<X'  respectively: 
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Fig.  A. 2  Energy  band  diagram  of  a  triangular  Si02  barrier.  El  and  E2  are  the  kinetic 
energies  referenced  to  the  silicon  conduction  band  minimum  in  the  initial 
state  XVSI  and  final  state  4^2  respectively.  Two  boundaries  of  the 
triangular  barrier  are  at  x=0  (abrupt)  and  x=X' (continuous). 
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(x>X')  (A. 19a) 


(0<x<X')  (A. 19b) 


k2(x)  is  the  local  real  wave  number  in  the  WKB  approximation  which  becomes  position- 
independent  at  x>Xox-  Note  that  (A.  19a)  and  (A.  19b)  are  related  to  each  other  by 
analytic  continuation  in  a  complex  x-plane.  Since  [k(x)]1/2°c[V(x)-E]1/4<*:(x-X')1/4,  an 
excursion  into  the  complex  x-plane  around  x=X'  picks  up  a  phase  7i/4,  which  appears  in 
(A.  19a).  A2=  (2k2/LS2)1/2  is  determined  from  normalization  of  W^. 

Both  (A.  19a)  and  (A.  19b)  become  invalid  near  the  classical  turning  point  x=X' 
where  k2(x)  and  k(x)  reduce  to  zero.  This  is  in  contrast  to  4^2  described  by  (A. 2a)  and 
(A. 2b)  for  a  trapezoidal  barrier  which  are  both  valid  near  x=Xox  due  to  the  abrupt 
potential  boundary.  Using  this  approximate  wavefunction  of  the  final  state  inside  the 
triangular  barrier  region,  it  would  be  straight  forward  to  compute  the  transition  matrix 
element  and  the  tunneling  current  formulae  using  the  same  approaches  demonstrated  in 
previous  sections.  The  result  for  the  dimensionless  tunneling  transition  probability  is: 

1                                               4(ms2/mox) 
P1P2  IM12  I  2    =   exp(-2u)    


(27T)2  moxk1/msl>:(0)+msl>f  (0)  /moxk1 

1 
=   exp(-2u)    TqME^  (A. 20) 

(27T)2 

The  only  difference  between  the  trapezoidal  formula  (A.l  1)  and  the  triangular  formula 
(A. 20)  is  the  pre-exponential  factor.  T0  in  (A.  11)  depends  on  kinetic  energies  in  both 
the  initial  state  (Ej)  and  the  final  state  (E2).  This  gives  T0  its  VGB  dependence  through 
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the  relationship  E2=E,-qV0X(VGB).   In  case  of  triangular  barrier,  T0'  in  (A. 20)  only 
depends  on  Ej  and  has  no  VGB  dependence. 

In  conclusion,  we  derived  tunneling  current  formulae  in  the  frame  work  of  time 
dependent  perturbation  theory.  The  result  shows  a  kinetic  energy  dependent  pre- 
exponential  factor  T0,  while  the  simple  barrier  transmission  coefficient  in  traditional 
formulae  gives  T0=T0'=1  [53].  For  trapezoidal  barrier,  T0  was  derived  rigorously  and  it 
gives  a  strong  gate  voltage  dependence  near  the  flatband  condition  which  is 
indispensable  to  account  for  the  experimental  tunneling  current-voltage  characteristics 
presented  in  chapter  4  [Fig.  4.13(b)].  For  triangular  barrier,  T0'  is  less  accurate  due  to 
the  breakdown  of  WKB  approximation  near  the  potential  boundary  x=X'  and  T0'  is 
independent  of  gate  voltage.  Therefore  in  this  high-field  triangular  barrier  range,  the 
traditional  simpler  tunneling  formula  with  T0=T0'=1  can  be  used. 

Summary  of  Detailed  Tunneling  Current  Formulae 

In  the  final  section,  we  list  the  detailed  tunneling  current  formulae  for  three 
tunneling  carrier  species  (electron,  valence  electron  and  hole)  and  two  shapes  of 
tunneling  barrier  (trapezoid  and  triangle),  which  can  be  easily  programmed  to  verify  the 
theoretical  curves  given  in  chapter  4.  These  formulae  give  the  explicit  dependence  on 
voltage  drops  in  the  silicon  and  oxide  layers  by  considering  (1)  direction  of  tunneling  in 
two  gate  voltage  ranges  VGB>VFB  and  VGB<VFB  (2)  the  dependence  of  Fermi-level  at 
the  Si/Si02  interface  on  silicon  band  bendings  and  (3)  parabolic  E(k)  relationship  inside 
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the  Si02  barrier.   Silicon  band  structures  are  taken  into  account  by  effective  mass 
approximations.  The  following  formulae  apply  to  both  barrier  shapes: 
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(A. 21)  is  for  conduction  band  electron  tunneling  from  <100>  silicon.  The 
density-of-state  mass  and  normal  mass  are  mj=0.417m0  and  ms=0.190m0  for  the  4-fold 
degenerate  transverse  valleys  and  mj=0.190m0  and  m^  for  the  2-fold  degenerate 
longitudinal  valleys.  For  valence  band  electron  tunneling  formula  (A. 22)  and  valence 
band  hole  tunneling  formula  (A.23),  the  density-of-state  mass  and  normal  mass  are 
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1  h / 


lh_ 


hh_, 


md  =0.16m0  and  m,!  -0.20m0  for  the  heavy-hole  band,  and  m[jn =0.45m0  and 


171^=0. 28m0  for  the  light  hole  band.  The  effective  mass  inside  Si02  barrier  is 
mox=0.437m0.  EFG  and  ECG  represents  the  Fermi-energy  and  conduction  band 
minimum  respectively  in  the  bulk  of  silicon  gate  material,  while  the  same  symbols  with 
subscript  'X'  instead  of  'G'  are  used  for  those  in  the  bulk  of  silicon  substrate  materials. 
VSG,  Vox  and  Vsx  denote  the  voltage  drops  in  the  silicon  gate,  Si02  insulator  and 
silicon  substrate  respectively.  Electron  and  hole  barrier  height  and  silicon  indirect 
energy  gap  are  (|)B=3.13eV,  (|)B  =4.25eV  and  E  =1.12eV,  which  add  up  to  8.5eV  for  the 
Si02  energy  gap. 

The  specific  forms  of  d(Ex)  and  T0(Ex,ms)  depends  on  the  barrier  shape  which 
are  listed  as  follows: 
Trapezoidal  Barrier 
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Triangular  Barrier 
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B'=     ^Bn-Ex 

(JN' 

=    ^Bp-Ex 

(JP) 

=     ^Bn  +  Ex  +  Eg 

(JVE 

!A.28) 


Ex'=    Ex+q|Vox| 


(JN-     JP) 
=   <5lvoxl-Ex-Eg         (JVE>  (a-29) 

<()B'  is  the  tunneling  barrier  height  for  a  carrier  in  the  right-hand  (x<0)  initial  state  with 
kinetic  energy  Ex,  and  Ex'  is  the  carrier  kinetic  energy  in  the  left-hand  (x>Xox)  final 


state. 
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